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Abstract 



This thesis report deals with the one dimensional Hubbard model. Specifically, 
we describe the quantum objects that diagonalize the normal ordered Hubbard 
hamiltonian, among those the so called pseudofermions. These pseudofermions, 
?7-spin and spin zero objects, are the scatterers and the scattering centers of the 
representation developed here. They have no residual energy interactions. The 
S-matrix of this representation can be written as a simple phase factor, which 
involves the phase shifts of the zero energy forward momentum scattering events. 
The form of the pseudofermion S-matrix constitutes an important new result of 
this thesis report. In contrast to the usual low-energy Luttinger liquid theory, 
the theory reported here allows us to categorize a separation of the charge and 
spin degrees of freedom at a finite energy excitation scale, of quantum objects 
called rotated electrons. The rotated electrons are linked to the electrons by a 
mere unitary transformation. 

Furthermore, we develop a pseudofermion dynamical theory and apply it 
to the evaluation of the spectral function in the one-electron removal and one- 
electron lower Hubbard band addition cases. For any value of the on-site effective 
Coloumb repulsion and electronic density, and in the limit of zero magnetization, 
we derive closed form expressions for these spectral functions, showing explicitly 
the emergence of the characteristic power-law type behavior of correlation func- 
tions of Luttinger liquids. We note however, that our derived expressions are 
valid for the entire elementary excitation energy bandwidth, and not just the lin- 
ear regime. We are able to identify practically all features of the spectral weight 
of the ID Hubbard model, in terms of pseudofermion or pseudofermion hole ex- 
citations. This brings new light to the understanding of the spectral properties 
of the one dimensional Hubbard model. 

The singular behavior of the theoretical spectral weight, as predicted by the 
explicitly calculated values of the relevant exponents and pre-f actors, leads to 
a spectral weight distribution which is detectable by photo emission and photo 
absorption experiments on quasi one dimensional materials. As an important 
contribution to the understanding of these materials, we are able to reproduce 
for the whole energy bandwidth, the experimental spectral distributions recently 
found for the organic compound TTF-TCNQ by high-resolution photo emission 
spectroscopy. This confirms the validity of the pseudofermion dynamical theory, 
and provides a deeper understanding of low dimensional correlated systems. 
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Chapter 1 
Introduction 



1.1 The one dimensional many-body problem 

This thesis deals with a model used to study strongly correlated electrons in 
solids, called the Hubbard model. The model was introduced by J. Hubbard in 
the article series on "Electron correlations in narrow energy bands" p]-|E]- We 
will exclusively focus on the case of one spatial dimension. 

In general, the Hubbard model is viewed as one of the most important models 
of strongly correlated electron systems in solids, and is used frequently in various 
applications both in one, two and three spatial dimensions. Formally, it presup- 
poses the existence of an atomic-like static lattice, upon which valence electrons 
exist in a single energy band. 

The Hubbard model is obtained for the approximative scheme where, due to 
the screening of the long-range Coloumb repulsion, its on-site Coloumb interaction 
term is large as compared to the inter-atomic interactions however not render- 
ing the inter-atomic hopping amplitudes negligible. Thus, the Hubbard model is 
loosely described as a dynamical model of electrons, where the Coloumb interac- 
tion (~ potential energy) competes with the transfer integral between neighboring 
lattice sites (~ kinetic energy), under the influence of the Pauli exclusion princi- 
ple. Despite its conceptual simplicity, the Hubbard model is indeed a very difficult 
model to solve exactly. It has only been solved in one dimension so far, by use 
of the Bethe ansatz technique, originally due to H. Bethe 0, where it was first 
applied to the isotropic Heisenberg chain. 

Even though no "true" one dimensional materials exist, many materials have 
been observed with a quasi one dimensional behavior. For example, some mate- 
rials exhibit very strong anisotropics, for instance in the electrical conductivity, 
where the motion of the electrons is confined to certain specific directions along 
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the Brillouin zone. As an example, "tetrathiafulvalene-tetracyanoquinodimethane" , 
abbreviated TTF-TCNQ ^Hj ^Z] |lH7j - P^7j . is a charge transfer salt consisting 
of linear stacks of planar TTF molecules and planar TCNQ molecules, which in 
the metallic phase exhibit an intra-stack conductivity three orders of magnitude 
larger than the inter-stack conductivity. 

Some other examples of materials that are adequately described by the Hub- 
bard model include V2O3 jH], strontium-copper oxide compounds (also known as 
"chain cuprates") such as Sr2Cu03 jH] 1101 and Sr2Cu02 [H] |12|, the Bechgaard 
salts (organic conductors) [IB] -[IS] and vr-conjugated polymers [IHj [TT?] . 

Historically, one example of the motivation to study such a model refers to the 
unusual metal-insulator transition theoretically predicted by N.F. Mott [201 [21] ■ 
This transition can not be explained by standard band theory j22j. One of the 
first materials observed to undergo a Mott-Hubbard metal-insulator transition, 
was the Cr-doped V2O3 [2S]-[2S]- The transition is characterized by the strong 
mutual interactions between the charge carriers, forcing them to be localized. 
This can intuitively be understood in the case of half filling (i.e. exactly one 
valence electron per lattice site). In this limit, it is not energetically favorable for 
electrons to be delocalized since hopping onto a neighboring lattice site, where 
another electron is present, costs more energy than staying put. In this way, 
a solid with an odd number of valence electrons per lattice site may exhibit 
insulating behavior, contrary to the predictions of standard band theory. 

Quantum systems exhibiting Mott related features continues to be one of the 
main topics in low dimensional strongly correlated electron systems, with some 
recent applications involving ID cuprates 26 , spin frustrated organic conductors 
[2Z], "orbital-selective Mott systems" (modeled by a multi-band hamiltonian) [2H1 
and dispersionless-boson interaction (modeled with a Holstein-Hubbard hamilto- 
nian) [2S]. 

From a theoretical point of view, the strongly correlated one dimensional sys- 
tems have spectral properties not explainable by the usual Fermi liquid theory. 
One of the trademarks of this theory, is the description of the low lying excitations 
in terms of "quasi particles", whose interaction is described by the scattering /- 
functions of Landau [nn[-[SS[- The ground state distribution in a Fermi liquid 
is given by the usual Fermi distribution, becoming a step distribution at zero 
temperature, where the step occurs at the Fermi level. The basic assumption of 
Fermi liquid theory is that there exists a one-to-one correspondence between the 
particle states of the original system and the quasiparticle states of the interact- 
ing system. Thus, there exists a strong 5-function coherent peak in the spectral 
function at the quasi particle energy (at the Fermi level). This picture breaks 
down if for example the interacting system produces bound states, as in a super- 
conductor. Indeed, according to the usual BCS theory, the formation of Cooper 



12 



pairs ultimately destroys the above mentioned one-to-one correspondence |34j . 

However, in the case of one spatial dimension, the breakdown of the Fermi 
liquid picture is more general as ID metals are characterized by the absence of 
Fermi liquid type quasi particles. Instead, the low excitation energy dynamics is 
described by charge and spin elementary excitations propagating independently of 
each other. These properties of interacting ID metals signal the advent of a new 
type of quantum liquid, known as the Luttinger liquid A characteristic 

of a Luttinger liquid is that it exhibits no coherent quasi particle peaks, as all 
low energy contributions to the one-particle spectral weight are of a non coherent 
origin and characterizeable in terms of the separated charge and spin degrees of 
freedom [IHI- Moreover, correlation functions decay algebraically as the Fermi 
level is approached, by an interaction dependent exponent. For example, in the 
case of the ID Hubbard model, the density of states decays with an exponent 
assuming values between and |, where the latter value corresponds to the 
strong coupling limit |41]. 

Using conformal invariance and finite-size scaling, H. Frahm and V. Korepin 
obtained the low lying spectral properties for the Hubbard model, evaluating the 
asymptotics of correlation functions |124j |125j . Since for the low lying elemen- 
tary excitations the Hubbard model belongs to the same universality class as the 
Tomonaga-Luttinger model, the exponents obtained in these references are the 
same as those of a Luttinger hquid. 

Examples of recent experiments related to ID Luttinger liquids 

Recently, there has been a surge of exciting new low dimensional materials, 
which commonly share the traits of Luttinger liquids. For example, in Ref. 
the conductive properties of carbon nanotubes are investigated and compared 
with the theoretical predictions of a Luttinger liquid. The authors find a power- 
law type scaling behavior of the conductance, with a measured value for the 
exponent in good agreement with the theoretical value. 

The material Lai.4_j;Ndo.6Sr2,Cu04 is studied for x = 0.12 in Refs. and 
|47j . and below and above this value in Ref. 148]. The reason for this mate- 
rial attracting attention is due to the quest for a better understanding of the 
copper-oxide superconductors. Unlike conventional metals, the charge carriers 
are confined to one dimensional "domain lines", but where the electronic spins 
in the region between these lines order themselves antiferromagnetically. This 
charge- and spin-ordered state is in these references called a " stripe phase" . This 
exotic charge transport suffers a dimensional crossover ID 2D as the criti- 
cal concentration grows beyond x = 0.12. Numerical calculations employing the 
Hubbard model confirmed the interpretation of "stripe phases" in this material. 

Ref. jini describes a one dimensional optical lattice of ultracold fermions in 
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a harmonical trap in order to study the Mott transition at micro Kelvin tem- 
peratures. In Ref. ^50^, such a lattice is set up by using "thousands" of parallel 
atomic waveguides, to form a nearly perfect atomic quantum wire. This new 
method provides an unprecedented control over the study of strongly correlated 
electron systems, and will be of great interest for anyone aspiring to work with 
low dimensional systems. Related numerical work is presented in Ref. jM], with 
the survival of spin-charge separation far outside the low excitation energy (Lut- 
tinger) regime, as one of the key results. 

Mentioned briefly, Ref. jS21 studies the angle resolved photoelectron spectra 
on metallic nano-wires. Ref. [SHI observes a spin-charge separation in quantum 
wires, which is readily cast into the Luttinger liquid scheme. Ref. jSl] studies the 
electronic transfer between various structures of the DNA double helix. 

The aim of this Thesis Report: 

In light of recent high-resolution photo emission experiments on various mate- 
rials, and in light of significantly improved experimental techniques for the study 
of strongly correlated electron systems in general, the absence of a dynamical 
theory for the theoretical model that is expected to describe the vast majority of 
the properties of these new materials - the Hubbard model - is a serious drawback 
for the complete understanding of many low dimensional quantum systems. 

Indeed, apart from the limit of infinite Coloumb interaction strength, so far 
it is only the properties of the low excitation energy regime that is understood 
to a sufficient degree of satisfaction, by employing theoretical tools only valid for 
that regime. 

However, many recent experiments (see for example the organic conductor 
TTF-TCNQ cited above) suggest the existence of elementary excitations at all 
energy scales, and not just the low lying ones. Thus, the need for a dynamical 
theory for the Hubbard model, capable of capturing the finite energy physics 
of the model hamiltonian in terms of a consistent scattering theory, is urgent. 
The recent high-resolution photo emission studies of the organic conductor TTF- 
TCNQ QZ] is a good example of this. Naturally, with the fine tuning of 
photo emission and photo absorption techniques, as well as with the development 
of new experimental procedures, such as the trapped ultracold atoms technique, 
the need for such a dynamical theory will most likely grow in the future. 

With this thesis report, we aim to derive and apply such a dynamical theory. 

For the derivation part, we will aim at calculating the spectral functions for the 
cases of one electron removal and one electron lower Hubbard band addition. This 
derivation will actually be intertwined with the development of the dynamical 
theory itself, and therefore constitutes the main part of this work, chapters Q 
and The aim of this disposition is to describe how a one electron spectral 
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function is derived, and at the same time gain a significant insight into the general 
physics of the model. However before this, in chapter we characterize the 
symmetries of the model and introduce the quantum objects that diagonalize 
the normal ordered hamiltonian, the so-called pseudofermions, related to the 
likewise introduced pseudoparticles. A scattering theory for theses pseudofermions 
is also developed in this chapter. In chapter (0), we use numerical studies to 
evaluate the one-electron spectral functions from the general expressions of the 
pseudofermion dynamical theory. The results of this chapter are then compared to 
recent photo emission experimental results (using the ARPES technique) for the 
organic conductor TTF-TCNQ, chapter (jH)). A short discussion with a summary 
of the main conclusions is presented in chapter ((Zj). 

1.2 The model hamiltonian 

For the following, and indeed throughout this thesis report, we will assume a 
vanishing magnetization m 0. Consider a solid whose ions are arranged in 
a crystalline structure. Since the ions of the solids are much heavier than the 
electrons, it is reasonable to assume the ions to be static, i.e. not participating 
in the dynamics of the solid. A general hamiltonian governing the dynamics of 
electrons with mass rrie and charge e can be written as 



where Xj denotes the position and pj the momentum of electron j (there are 
a total of N electrons present in the system). V{xj) is a potential with the 
periodicity of the lattice structure. We will assume a maximum of two valence 
electrons per lattice site (with opposite spin projection). This means that we are 
considering atoms where all the other electrons are strongly bound to the ions 
and hence do not contribute to the dynamics. 

The eigenfunctions Xfc(^) of the one-body term of the hamiltonian (jl.lj) are 
just Bloch functions, i.e. Xk{x) = e*'*'^Mfc(a;), where Uk{x) is a function with 
the periodicity of the lattice. These can in turn be expressed in another basis, 
namely the Wannier basis, comprising the Wannier eigenfunctions iplx — Rj) that 
are centered on the lattice site at position Rj. The relation between the functions 
of the two bases is: 






(1.2) 
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where L is the total length of the lattice (many times, units such that the lattice 
constant a is equal to 1 are employed, and hence L = N). 

By using wave functions centered on the lattice sites, we may express the 
second quantization annihilation field operator \I'o-(a;), as a sum over the entire 
lattice: 

N 

^,{x) = Y,,lj{x-R^)c,, (1.3) 

j=i 

where we introduced the second quantization electron annihilation operator Cjo-, 
that annihilates an electron on lattice site j with spin projection a. The creation 
field operator ^J.(a;) is nothing but the hermitian conjugate of the operator given 
in equation ()1.3p . 

In terms of the field operators, the first- quantization hamiltonian above may 
be rewritten as 

H = Y,jdx ^i{x) f m„{x) + J2fdxdy ^l{x)<ijl{y) U ^,{y)^,{x) (1.4) 

where T and If correspond to the first-quantization one-body (kinetic) and two- 
body (Coloumb interaction) operators, respectively. Thus, in terms of Wannier 
states: 

H = Y^ Ujclcj, + ^J2Y1 Uijkiclc]^CkaCia (1.5) 

ij,a ijkl <T(7 

where we have defined the kinetic and the potential matrix elements according 
to: 

(1.6) 

Uijki = j dxdy ip*{x - Ri)ip*{y - Rj) U ijj{y - Rk)'ip{x - Ri) 

Note that in this derivation we have all along assumed a single band model. 
For multi-band solids we would have to use an extended model, except for the 
special case of weak inter-band interactions. 

Let us assume the Wannier states to be "strongly localized", i.e. that tplx — 
Rj) is centered on Rj. This means that the Wannier states of two neighboring 
lattice sites may have a finite overlap with each other, but lattice sites further 
apart may not. Hence, in the first term of Eq. ()1.5|) . the summation over the 
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lattice sites reduces to a summation over nearest neighbor pairs, denoted {ij) 
(where i = j ±1). 

Due to the "strong locahzation" of the Wannier states, we can assume that 
for the two-body term, only matrix elements referring to the same site become 
non-negligable. This makes the Coloumb interaction term become an effective 
potential: the Coloumb repulsion between two electrons on the same lattice site 
dominates completely over the repulsion of two electrons on different lattice sites, 
due to the screening within the actual material sample. 

Introducing these assumptions into Eq. ()1.5|) . we obtain an effective hamilto- 
nian where tij = tij6ij±i and Uijki = Unudijki- 

H = ^ ^ tijcl^Cja + — ^ ] ^ ] Uiiiicl^cl^Ci^Cia (l-^) 
{ij)u i CO" 

Since we are considering systems with identical lattice sites, we assume isotropic 
nearest neighbor transfer amplitudes tij = —t and effective Coloumb interaction 
strengths Una = 2U (where the factor of 2 comes from the spin summation), 
which finally leads us to the Hubbard model hamiltonian: 

H=-tJ2 4aCja + UJ2 (1-8) 
(ij)(r i 

where Ui^ = c\^Ciu is the electron number operator. As we will see in the next 
chapter, the interaction term is written in a different way in order to comprise 
useful sjTumetries, such as the particle-hole symmetry. Note that the i = j term 
in the kinetic part only counts the number of electrons, and can hence be absorbed 
by the chemical potential in a grand canonical ensemble description. 

The hamiltonian derived above allows for some characterization merely by 
"hand waving" arguments. For example, one can see that in the limit {U /t) 0, 
the electrons do not interact and hence are delocalized for almost any value of the 
band filling n = N/Na, where N is the number of electrons and A^^ the number 
of lattice sites (in subsequent chapters, however, we will define n = N/L, where 
L = aNa is the length of the lattice chain, where a is the lattice constant with 
dimension of length). Exceptions are the fully polarized half filled case, and the 
fully occupied case n = 2, respectively. This is the tight-binding model, for which 
the electrons disperse in a cosine-like band 

In the opposite limit, [U/t) —>■ oo, the electrons are localized. At half filling, 
the ground state of the model describes an antiferromagnetic insulator Ac- 
tually, at half filling the Hubbard model is insulating for all values of (U/t) > 
|58j . Moreover, in this limit, any eigenstate of the model can be factorized into 
two eigenstates: one describing spinless fermions and another describing charge- 
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less spins jHH]. It can be shown that the effective hamiltonian for the spin part, 
is nothing but the ID antiferromagnetic Heisenberg spin hamiltonian |13], with 
an effective coupling constant equal to n(4t^/f/)(l — sin[27rn]/[27™]). 

This result can be understood by second order perturbation theory around 
{U/t) = oo. Since the spin configuration is antiferromagnetic, one can think of 
virtual states in which an electron hop to one of its neighbors, with the amplitude 
proportional to t. However, due to the large repulsion U, one of the electrons 
immediately hops back to the originating site, again with an amplitude t. This 
gives rise to an antiferromagnetic exchange of strength ~ t ■ {1/U) ■ t = t'^/U. 
When away from half filling, the value of the effective interaction of the spins will 
then be modified due to the presence of holes. 

Due to the insulating properties of the Hubbard model at half filling, we will 
always stay away from this value of the density n in this thesis report. 
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Chapter 2 



The Hubbard Model 



2.1 Symmetries and exact solution 
2.1.1 SO (4) symmetry 

By rewriting the model hamiltonian of the previous chapter, it is possible to 
arrive at the following general hamiltonian for electrons on a one dimensional 
chain consisting of A^^ equally spaced lattice sites: 

H = f + U=-tJ2 clcja + U Y^icl^c,^ - l){clc,^ - I) (2.1) 

{ij)a i 

where c]^ is the usual electronic creation operator and Qo- is the usual electronic 
annihilation operator. We will assume Na to be an even number. The symbol (ij) 
indicates that the summation is done over neighbouring sites only, i.e. j = i ± 1 
and a stands for the electronic spin projection. We consider periodic boundary 
conditions. The transfer integral t and the effective on-site coloumb interaction 
strength U were introduced in the previous chapter. The hamiltonian conserves 
the total number of electrons, as well as the total number of f-spin electrons 
and the total number of |-spin electrons even though, in the following, we will 
concentrate on somewhat more subtle and interesting symmetries. Any one site 
can either be empty of electrons, singly occupied by a |-spin or a J,-spin electron, 
or "doubly occupied", i.e. occupied by two electrons (with opposing spin projec- 
tion due to the Pauli principle). The number of electrons N, can be written as 
N — N-^ + N^, where A^^. is the total number of cr-spin electrons. The creation 
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and annihilation operators have the usual Fourier representation on a lattice: 



^E^^'^V (2.2) 



The total momentum operator can thus most easily be written: 

^=EE^-(^)^ (2-3) 

where Na{k) = c^ka^ka is the Fourier transformed number operator, and L = 
aNa is the physical length of the lattice chain [a is the lattice constant). This 
hamiltonian is many times written as Hso{a) due to the fact that it is invariant 
under the SO (4) symmetry group [HS] ISEj- The SO (4) group is isometric with the 
SU{2)^SU (2) group and differs in that only half of the irreducible representations 
of SU(2) SU{2) correspond to energy eigenstates, as further explained below. 
Therefore, symmetry properties of the model are usually explained in terms of 
the two related independent SU{2) symmetries, namely the SU{2) spin algebra 
and the SU{2) r^-spin algebra. The generators of the spin algebra (subscript s) 
and of the yy-spin algebra (subscript c) are 



2 \ / ^2 
S\ = Y.AfA S\ = Y.^-inA, (2.4) 

i i 

Ss = Y^ c\^Cii Sc = E(-l)'Q|Ci| 

i i 

where N = J2a iVo- = J2i '^la^^io^ t^e total electronic number operator 

and the a-spin electronic number operator, respectively. They satisfy the usual 
SU (2) commutation relations: 



of A 



2Sl a = c,s (2.5) 



Apart from commuting with the hamiltonian p.lj) . all the generators of the 
spin algebra commute with all the generators of the ?7-spin algebra. Together 
with the square of the total spin and total r^-spin operator, {S°'y (where a = c 
or s), the diagonal generators S^ and the hamiltonian form a set of commuting 
operators. The transformations associated with the SU{2) symmetries are the 
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spin flip (c]^ 4-cr) the //-spin flip (c|^ —>■ (— l)'cjo-), respectively (the latter 
transformation is associated with the so called particle-hole invariance). As is 
trivially understood by physical reasoning, these symmetries are easily broken, 
for example by applying an external magnetic fleld or by introducing a chemical 
potential, since in these cases the full hamiltonian will not commute with the 
off diagonal generators and thus the system will prefer a certain spin projection 
and a certain 77-spin projection, respectively. Mathematically, this comes about 
by adding a term to the hamiltonian equal to fJ^aS^, where in the case of ?7-spin 
[a = c) fic = 2/i where n is the chemical potential and in the case of spin {a = s) 
fis = 2/io/i where /io is the Bohr magneton and h the strength of the applied 
external magnetic fleld. 

Applying any of the ladder operators, 5*^ or 5*0, (where a = c or s), to an 
eigenstate of the model ()2.1|) . leads either to zero or to a new state for which the 
total spin projection (a = s) or the total ?7-spin projection (a = c) has changed 
by one. One can then successively apply this ladder operator until we reach either 
the "top" of the ladder (by repeatedly applying SjJ or the "bottom" of the ladder 
(by repeatedly applying Sa)- The state reached by successive application of the 
operator Sa is usually called a lowest weight state (LWS), whilst the state reached 
by successive application of the operator S"^ is usually called a highest weight state 
(HWS), respectively. Trying to climb lower than the bottom of the ladder, as 
well as trying to climb higher than the top of it, yields zero: S'a|LWS) = and 
S^IjIHWS) = 0. The eigenvalue of 5"^, denoted by S^, ranges from —Sa, —Sa + 1 , 
. . . , Sa — 1, Sa, where {8°")^ has eigenvalue Sa{Sa + 1). 

Because of the above found symmetries, it will be sufficient to consider elec- 
tronic densities n = ^ and magnetization m = ^^1-J1l^ where N = N-^ + is 
the number of electrons, such that < raa < 1 and < ma < na. For example, 
considering electronic densities larger than one is equivalent to considering hole 
concentrations smaller than one, and thus we are back at the original mathemat- 
ical formulation. In the remainder of this thesis report, na and ma will always 
be assumed to obey these inequalities. There is one final remark to be made, 
originally due to the discoveries of Ref [HZj '■ the hamiltonian would not have com- 
muted with the generators of the //-spin algebra had we chosen the number of 
lattice sites Na to be an odd number (this is easily verified by explicit calculation 
of the commutators), which is why Na is assumed even. 

When forming S^ -\- S^ = iV^ — we see that S^ + SI only takes integer 
values. This implies some restrictions in the types of multiplets that are allowed 
by the model: states for which both 5*^ and S^ are integers, or for which both are 
half-odd integers, are allowed, whilst states for which one is integer and the other 
is a half-odd integer, are prohibited. This is the reason for why the hamiltonian 
()2.H) does not possess a full SU{2) ® SU{2) symmetry. 
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2.1.2 The Bethe ansatz solution 



By using an ansatz wave- function, Lieb and Wu managed, in their famous paper of 
1968 to reduce the problem of diagonahzing the hamiltonian ()2.H) into solving 
a set of coupled non-linear equations. The variables in these equations are two sets 
of numbers usually referred to as charge momenta {kj} and spin rapidities {A;}. 
These numbers can be finite or infinite and are in general complex. However, it is 
these two sets of numbers that characterize the eigenfunctions of the model. As a 
side note, it is worth to mention that this eigenf unction, the "Bethe ansatz wave 
function" was explicitly presented by F. Woynarovich in 1982, for (U/t) ^ 1 [HH] . 
Unfortunately, for finite values of (U/t), it is not suitable for direct calculation 
of correlation functions, due to its complexity. The coupled non-linear equations, 
also called the "Lieb-Wu equations", are: 



h T -i-r Ai — sm k^a — 

e^'^ = n W J = l,...,N 

fj^ Xi - sinkja + ij 



nXi — sm kja — t-t Ai — Am — i^ , , 

J 1 4_TTJ 2 / = i^...^M (2.6) 



where M is the number of |-spin electrons. Here and in the following, we use the 
notation of T. Deguchi et al [00] (with the exception of n and m in that reference, 
which we here denote by n and m, on order to separate them from the density 
and the magnetization). Taking logarithms on both sides of the two equations 
fl2.6j) . introduces the quantum numbers, which are integers or half-odd integers. 
Even though the solution provides us with the ground state energy as well as 
the energy spectrum, it does not provide us with an association of the quantum 
numbers with the electrons, nor to the configurations of electrons on a lattice. 

Moving on, in 1972, Takahashi reformulated the original solution using the 
so called string hypothesis P|. In this paper, the hypothesis is used to classify 
the finite quantum numbers of the problems into "strings", which are valid as 
the system size Na becomes very large. Basically, the numbers kj and A; are 
grouped into strings according to the value of their real parts, and are distributed 
symmetrically with respect to the real axis. Like this, Takahashi arrives at three 
different groupings of numbers; one that only involves real charge momenta kj 
(type I), one that only involves (complex) spin rapidities A/ (type II), and finally 
one that involves both complex charge momenta kj and complex spin rapidities 
A; (type HI). By using these relationships between the numbers of equation ()2.6p . 
we arrive after some algebra to the following set of equations, one for each type. 
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usually called the " discrete thermodynamic Takahashi equations" : 



, ^ ^ T ^ v-^ ^ / . sin k^a — A" \ ^ ^ / , sin kM — A'" \ 

kjL = 2nij ^^^^^^ (4 — ) - 2 E E (4 — mr^j 

n=l a=l ^ ' n=l a=l ^ ' 

N-2M' / ■ 7 An\ °° -^"^ / A" \fh\ 

-2 E arctan 45!ii|^ ^ + E E 9" (*^) 

j=l ^ ^ m=l 13=1 ^ ^ 

^ (^arcsin(A'^ + i^) + arcsin(A'^ - i^)^ = (2.7) 



Here a enumerates the string of length m (where "length" translates into 
" amount of numbers on a given string" ) . A^ and A'^ are the purely real midpoints 
of strings of type II and type III, respectively. and are the total numbers 
of strings of type II (with length n) and of type III (with length m), respectively. 
Finally, M' is the total number of charge momenta numbers involved in a string 
of type III, M' = YlnLi ^^n- The number of cr-spin electrons are related to the 
numbers of strings through: 

00 

n=l 

N^ = N-N^ (2.8) 



and the function Qnfh{x) is given by: 



©nm('^) ^ 



2 arctan + 4 arctan (^z^i) + ■ ■ ■ 

• • • + 4 arctan (^^) + 2 arctan (^) , n ^ m 
4 arctan (|) + 4 arctan (f ) + ■ ■ ■ 
. . . + 4 arctan (2^:^) + 2 arctan (^) , n — m 



(2.9) 



It is important to note that these equations introduce the quantum numbers 
of the Bethe ansatz. These are a set of purely real numbers: Ij, J" and J'". Ij is 
an integer if ^^(M^ + M'^) is even and a half-odd integer if odd, J" is an integer 
UN — Mfi is odd and a half-odd integer if even, J'^ is an integer if A^^ — A?" -|- 
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is odd and a half-odd integer if even. They obey the following inequalities: 

|2/,|<A^a 

oo 

\2J2\<N -2M' -l-J^t^nM^ (2.10) 

m=l 

oo 

\2JZ\ <N,-N + 2M'-1-J2 tmnMi, 

m=l 

where tmn = 2 min(m, n) — Smn- The quantum numbers are equidistant from each 
other, for example Jj+i — Ij = 1. With these numbers specified, one can solve 
Eqs. ()2.7p for the numbers kj, and A'^. The specification of the occupancy 
configurations of the above quantum numbers allows for the construction of the 
ground state and any excited state of the original hamiltonian |BT] . 

Some very important insights into the Bethe-ansatz solution was given by 
F.H.L. Efiler, V.E, Korepin and K. Schoutens in 1992 [02] - [HII. In Ref. ^ and 
[nni it was shown that the Bethe ansatz solution only accounts for either lowest 
weight or highest weight states of the one dimensional Hubbard model. However, 
in Ref. it was shown that, after taking into account all the states reached by 
the off-diagonal generators of the spin- and the ?7-spin algebras, the Bethe ansatz 
solution is indeed complete, in that the total number of states present in the 
solution, gives the accurate dimension of the Hilbert space of the original model. 
In this notation, the total energy and momentum (modulo 27r) are expressed as: 

E = -2t J2 cos%a + 4t^^ReWl - ( A'^ + 2^ j +UN (2.11) 

j=l n=l a=l V V / 

N-2M' oo A/4 

j=l n=l a=l 

where N = (iV^ — 2iV)/4 and Re denotes the real part. 



Re { arcsin ( A'" -|- 



(n+1)- 
a 



(2.12) 
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2.2 Pseudoparticles and rotated electrons 



2.2.1 Pseudoparticles - historical overview 

It was mentioned earlier that the Bethe ansatz does not provide information 
about the connection of its quantum numbers to the original electrons of the 
problem. However, these numbers have been associated with various quantum 
objects, different from the electrons, in various descriptions valid for some strict 
subspace of the Hilbert space of the model. Some examples include: charge 
pseudoparticles and spin pseudoparticles, introduced in Refs. - [ZEI, in the 
study of low-lying excitations. 

In 1990, M. Ogata and H. Shiba used the Bethe ansatz wave-function factoriza- 
tion to calculate physical quantities, such as the momentum distribution function 
and the spin correlation function, for the case of strong coupling, f//t — > cx3 [77j . 
They used the fact that, in this limit, the Bethe ansatz ground state wave func- 
tion factorizes into a charge part and a spin part. The charge degrees of freedom 
are then calculated via a Slater determinant of "spinless fermions" and the spin 
degrees of freedom are described by the one dimensional 5" = | Heisenberg model. 
This framework, with spinless fermions and a "squeezed" spin wave taken from 
the 5* = I Heisenberg model, was subsequently used by K. Penc et al in various 
publications [TH] - in the study of the exact one electron spectral function 
for {U/t) oo. K. Penc and B.S. Shastry then adopted a similar technique for 
the Schultz-Shastry model [HI]. A related representation valid in the (U/t) oo 
limit was presented by R.G. Dias and J.M.B. Lopes dos Santos |9Uj . 

In Ref. jHE] the SO (A) symmetry of the model is used in case of exact half 
filling na = 1, to describe the excitation spectrum and an S-matrix in terms 
of spinon and holon excitations. In this work, the elementary scatterers carry 
either spin but no charge (thus, they are dubbed spin -(-| and spin — | spinous, 
respectively) or charge but no spin (consequently dubbed r^-spin +i and r^-spin 
— I holons, respectively). Note that these "holons" and "spinous" are not the 
same quantum objects as those which will be introduced in this thesis report. 

The common ground between these works is the fact that they are only valid 
for some strict subspace of the model, and that the introduced quantum object 
description, either derived from the numbers introduced by the string hypothesis 
or by considering the symmetries of the model, accounts for the famous separation 
of the electronic degrees of freedom: the spin and the charge degrees of freedom of 
the electron are described by different quantum objects that propagate through 
the system with different velocities. 

In 1997, J.M.P. Carmelo and N.M.R. Peres managed to generalize the pre- 
vious pseudoparticle picture, to be valid for the entire Hilbert space j^I- The 
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energy bands and the residual interactions of these pseudoparticles were exphcitly 
presented. Furthermore, a normal ordered (relative to the ground state) formu- 
lation of the problem allowed the hamiltonian to be rewritten solely in terms 
of momentum configuration distribution operators. However, the connection be- 
tween these pseudoparticles and the original electrons remained unknown until 
2004, when the electrons were "re-connected" to the problem by relating them 
to quantum objects baptized rotated electrons, obtainable from the electrons by 
a mere unitary transformation for all values of U/t jSB]- The pseudoparticles are 
then " constructed" in terms of the decoupled spin- and charge degrees of freedom 
of these rotated electrons. 

To give a flavour of how the pseudoparticles of Ref. jH2! were born out of 
the quantum numbers of Takahashi, and also to conform with the notation that 
we will use from now on, we let a — > j and n —>■ u in the Takahashi string 
hypothesis, Eq. ()2.7j) . Furthermore, we will denote the numbers Ij, and 
J'^ by J|°, Ij'^ and Ij'^ respectively. Correspondingly, let kj, and A'2 be re- 
baptized into kj (unchanged), Agu and Aa, respectively. Since these quantities 
depend on each other, according to j Ij^ kj (and similarly for the other two 
quantum numbers Ij'^ and Ij'^), we actually have by Eq. ()2.7p . that kj = k{Ij^), 
Agty = Asui^Ij") and Acu = Acv^Ij"). Associating the quantum numbers with 
pseudoparticle momenta qj, according to qj = {27i/L)I°"^ (where 1°"^ is shorthand 
for any of the three types of quantum numbers), we can finally write k = k{qj), 
Asu = Asu{qj) and Acu = Acv{qj) respectively. One should not confuse the different 
gj's with each other: they are different in that they are equal to (27r/L)/; 
{271/1)1^" and (27r/L)J|^ respectively. In short: 

Ij — > kj kj 



cO 
j ' 



J^^Jf A^^A,, (2.13) 

j/n jcu j.in . 

and 

k = k{qj) 

A,u = AUqj) Qj = ^ir (2-14) 

Acu = Acu{qj) 

In the following, we will write au as a collective symbol standing for all of 
the cO,cz/ and su branches (unless otherwise specified). This change of notation 
helps us to identify the sets of quantum numbers as occupational configurations 
of pseudoparticles, whose discrete momentum spacing is the usual 2ti/L. Fur- 
thermore, it identifies some functions of these momenta, collectively referred to 



(Ij = 


2vr 


(Ij = 


L ^ 
2vr 


Ij = 
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as rapidities, namely k, Asu and Acu, that are related to each other via Eq. ()2.7|) . 
The inequalities of Eq. ()2.1()|1 define the minimum and maximum values of the 
quantum numbers, and thus defines the occupancies inside the effective Brillouin 
zone of the pseudoparticles. For the ground state, the functions k{qj) and Aay{qj) 
are odd functions of their arguments |94j . 

2.2.2 Rotated electrons - historical overview 

The mapping that transforms the electrons into their rotated counterparts, the 
rotated electrons, follows the previous work by A.B. Harris and R.V. Lange jHl] 
and the work of A.H. MacDonald et al |^. In these publications, a unitary 
transformation from now on denoted V{U /t), is introduced for large values of the 
on-site coloumb repulsion U, that cancels all terms in the original hamiltonian 
that change the number of doubly occupied sites. Later, the same transformation 
is successfully used in Ref. [SI?, in order to calculate the one electron spectral 
function in the limit U /t oo. Because of the "large U history" behind the 
introduction of this transformation, it will be our starting point as well. In this 
limit, V{U/t) can be written as an expansion in powers of (t/U), but we should 
note that formally, this expansion is not the definition of V{U/t), which we will 
give later. 

The basic consideration behind this transformation is the fact that double 
occupancy is a good quantum number in the limit U/t ^ oo. This is easily 
seen by investigating the Coulomb interaction term of the hamiltonian ^2.1\ : a 
doubly occupied site gives a contribution to the total energy of the system equal 
to f/, whilst other types of occupancies does not. Thus, states with j and j ± 1 
number of doubly occupied sites have an energy difference equal to [/, so that 
this difference goes to infinity as U . 

The double occupancy quantum number D, is nothing but the expectation 
value of the double occupancy operator, D = (D), where 

Na 

D = ^ni^nii, ni„ = c\^Cia (2.15) 

i=l 

In the following, we will adopt the following notation for an arbitrary operator 
X, transformed by V{U/t): 

X = V\U/t)XV{U/t) (2.16) 

where X is the corresponding "rotated operator". Later we will see that V{U/t) 
is in fact unitary, which means that X and X share the same set of eigenvalues 
and preserves the norm of the eigenstates, which will turn out to be very useful. 



27 



The kinetic term of (j2]Tj), responsible for the dynamics of the model, can be 
written as a sum of three terms, T = Tq + Tu + T^u, according to the energy dif- 
ference that the hopping results in ("before" the hopping as compared to "after" 
the hopping). For example, a lonely electron on site i that hops to a neighbour- 
ing site i ± 1, where one electron already is present, will have increased the total 
energy of the system by U. Thus, the job that V{U/t) has to accomplish is to 
cancel Tu and T^u- It is then possible to rewrite the hamiltonian according to 

H = H^^^ + H^^^ + H^^^ + ... (2.17) 

where H^^'' allows hopping with a total number of j doubly occupied sites. One 
should note that Eq. ()2.17|) is actually nothing but a large - U expansion of the 
original hamiltonian, where H^^^ is the U = oo term, and the following terms 
are corrections of order {t/Uy . The first eight terms in the series ()2.17|1 were 
explicitly calculated and presented in Ref. |H5] . 

The question arises, then, how one obtains the operator expressions that con- 
stitute the terms of ()2.17|) . To answer this, we remember the fact that any unitary 
operator can be written as the exponential of an anti-hermitian operator, which 
here will be called Y{U/t): 



H =V{U/t)HV^{U/t)) 
V(U/t) = e^(^/*) J 



^ = (2.18) 



Now, using the Baker-Hausdorff Lemma, we see that this can be rewritten as 



H = H 



Y{U/t),H +- Y{U/t), Y{U/t),H +... (2.19) 



which immediately produces very many terms as we need to introduce H = 
Tq + Tu + T^u + U (note that U = UD). Requiring {D) to be a good quantum 
number, we need to have 

D,H =0 (2.20) 
Now, assuming that the following series expansion exists 

Y = + + + _ _ _ (2.21) 

where Y^^^ ~ {t/Uy, we can express Eq. ()2.19p by using Eq. 1)2.211) together with 
H = Tq + Tu + T-u + U . The assumption of the existence of Eq. 1)2.211) is actually 
not such a gamble as it may seem: just like an ansatz solution to a differential 
equation "if it works it works, if it does not it does not". In this case "works" 
translates into "deriving a closed form expression for y(^)" . By retaining only the 
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first term of tfie expansion of Y in the above mentioned sclieme, one can deduce 
tliat in order for Eq. ()2.2()j) to be valid, we need to have: 

= ^{Tu- f^u) (2.22) 

Inserting into Eq. ()2.19p with an added unknown Y^'^\ we find after some algebra 
that in order for Eq. ()2.2()j) to hold, we must have 



y(2) ^ ±_ 
f/2 



(2.23) 



Tc/,T_ 



u 



(2.24) 



and so forth. In this way, the first terms of Eq. ()2.17p becomes: 

This expansion can be continued to higher orders in {t/U), finding closed form 
expressions of V{U/t) that are successively valid for a larger range of (t/U) values, 
but it is actually not needed in order to create a valid theory for arbitrary values 
of this parameter. What is enough is the overall definition of the transformation, 
which is a combination of Eqs. fTT^i . ^TK)^ and fT^ : 



ViU/t) 
D,H 

H = H + 



.y(u/t) 



(2.25) 



Y{U/t),H 



YiU/t), Y{U/t),H 



+ 



Even though this definition seems quite abstract, there are some things that can 
be said about V{U/t) by pure physical reasoning. Since electronic double occu- 
pancy is a good quantum number for [U /t) = oo, electronic |-spin and electronic 
J,-spin single occupancies, as well as electronic no occupancy are also good quan- 
tum numbers in this limit (we will come back to this in section 1^. 2. 3|) . For rotated 
electrons, however, these numbers are always good quantum numbers. This means 
that when a finite-(f//t) energy eigenstate is acted upon by the operator V(U/t), 
the resulting state must bear some similarities with the corresponding [U /t) = oo 
eigenstate. This also implies that V{U/t) 1 when (U/t) oo due to the sim- 
ple fact that in this limit, the rotated electrons are identical to the electrons, or 
in other words, y = in the expansion given by Eq. 1)2.211) . The formal proof of 
the existence, uniqueness and unitariness of V{U/t) is given in Ref. [HSI- Lastly, 
since V{U/t) does not change the lattice (neither the amount of lattice sites nor 
the lattice constant), we must have that 



V{U/t),P 







(2.26) 



considering that the momentum operator is the generator of lattice translations 
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2.2.3 Connecting rotated electrons to pseudoparticles 

By the definition ()2.25|) . we have actually defined a new quantum object, dubbed 
the rotated electron, that has the same spin and charge and that exists in the 
same lattice as the original electron, but for which double occupancy is a good 
quantum number for all values of U / 1. Let us look a little bit deeper into this 
claim. 

Define the number operators for electronic double occupancy, no occupancy, 
single t-spin and single j-spin occupancy, respectively: 

Rc- = c|^Ct|c[| Q| 
i 

Rc+ = ^ Cj|c|^Q|c|| 
i 

4- = $^c;^Qtc;^Qi (2.27) 

i 

Rs+ = ^ C-^Ci|cJ^Cj| 
i 

These operators are not independent. In fact, they can all be expressed in 
terms of the electronic double occupancy operator D: 



Rc- ■ 
Rc+ ■ 
Rs- ■ 
R.«-t- ■ 



D 

-■Na-N + D 
-■N -N^-D 



(2.28) 



Note that this description of the different types of electronic occupations will be 
very important for the framework that we are going to use. Indeed, it accounts 
for all the electrons in the system, 2Rc- + Rs+ + Rs- = N. Now, let us define the 
rotated electronic creation and annihilation operators: 



V{U/t)clV\U/t) 
V{U/t)c,,V\U/t) 



(2.29) 



The operators and Cjo- create and annihilate some quantum objects whose 
double occupancy number, single occupancy f-spin and single occupancy |-spin 
number, as well as its no occupancy number, are all good quantum numbers for 
any value of {U /t). This is seen by inserting Eq. ()2.29p into Eq. ()2.27|) . and thus 
forming the corresponding rotated number operators, Rau where a = c or s and 
/ = — , +. Like this we obtain by definition ()2.25j) and by Eq. ()2.28|1 that: 



H, Rc- 




H, Rc+ 




H, Rs- 




H, Rs+ 







(2.30) 
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According to the studies of Ref. |H3|, it is the separated charge and spin degrees 
of freedom of the rotated electrons that constitute some exotic quantum objects 
called pseudopartides. The occupancy configurations of these pseudoparticles are 
described by the quantum numbers given by the Takahashi string hypothesis [6T] . 
The connection to the original electrons was not reached because of the fact that 
the relevant electronic occupational numbers were not good quantum numbers, 
so that any separation of the original electronic degrees of freedom is bound to 
lead to quantum objects consisting of several very involved superpositions of the 
original electronic occupational configurations. The picture becomes much more 
elegant if, when describing the various excitations of the model, the starting point 
consists of quantum objects whose occupancy number operators commute with 
the original hamiltonian, like those of the rotated electron. Like this, it will 
be much more simple to describe the properties of the new quantum objects, 
here called the pseudoparticles, in terms of the occupancy configurations of the 
previous ones, the rotated electrons. 

To proceed any further, it is necessary to explain what the pseudoparticles are 
constituted of. Since V{U/t) is unitary, all eigenvalues of the electronic number 
operators of Eq. ()2.27p are equal to the eigenvalues of the corresponding rotated 
operators, Rai- Now, there are N electrons in the system together with A^^ = 
2Na — N electronic holes. Let the number of sites that are singly occupied by 
rotated electrons be equal to N^. This means that there are N — Nc rotated 
electrons on doubly occupied sites, and hence a number of (A^ — Nc)/2 doubly 
occupied sites. By the same reasoning, we have {N^ — Nc)/2 sites doubly occupied 
by rotated holes ("empty sites"). Next we baptize the |-spin |-spin pair, on the 
sites doubly occupied by rotated electrons, — | holons. There is a total number 
of = {N — Nc)/2 of such quantum objects. Whilst the — ^ holons are spin 

zero objects, the value of its 77-spin projection is — |, which justifies the choice of 
name. Equivalently, we form a total number M^_,_i = (A^'' — Ac)/2 of +| holons, 
from the sites doubly occupied by rotated holes. The Nc rotated electrons on the 
singly occupied sites decouple into Nc chargeons (with the same charge as the 
rotated electron but with no spin degrees of freedom) and Nc spinons (with the 
same spin as the rotated electron but with no charge degrees of freedom) . Now, 
the total number of j-spin spinons is ^1 = N-^ — {N — Nc)/2, which is nothing 
but the total number of j-spin rotated electrons of the system, minus the f-spins 
of the doubly occupied sites. Equivalently, = N^ — {N — Nc)/2 is the total 

number of |-spin spinons. Below, we will make use of the number A"^ defined 
as the number of lattice sites not singly occupied, i.e. either doubly occupied 
or empty: A"^ = Na — Nc- Lastly, we note that the Nc sites singly occupied by 
rotated electrons also carry A'c rotated electronic holes. The charge part of these 
rotated holes, living on singly occupied sites, will be called antichargeons. To 
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summarize, we have the following amount of quantum objects: 



i //-spin holons M^_i={N - N^) /2 



t r^-spin holons + 1 = (A^^ - N^) /2 

i-spin spinous _ i = A^^ - (A^ - N^) /2 (2.31) 

t-spin spinous +i=N^ - {N - N^) /2 



The following should be observed: 



1 , 1 =N^ 



Yl = + = Na (2.32) 



So far we have only regrouped the rotated electrons, giving them new names 
according to the occupancies of the lattice sites. Following the interpretation of 
Ref. [Sni; let now the Nc chargeons and the Nc antichargeons recombine into 
= Nc cO-pseudoparticles. Let further a total number of u +| spinous and 
u — I spinons form one su pseudoparticle, where the number z/ = 1,2,... . 
Equivalently, let a total number of u +^ holons and i> — | holons form one 
cu pseudoparticle, where the number u = 1,2, . . . . 

The claim of Ref. [S2] is that the quantum numbers describing the momentum 
occupancies of these created pseudoparticles, are nothing but the Takahashi num- 
bers, given by the string hypothesis. Thus an au pseudoparticle (where a = c 
or s and u = 1,2, . . .), that contains a number u of holon (a = c) or spinon 
(a = s) pairs, corresponds to a string with length z/ of type II (a = s) or of type 
III (a = c). There will be a total number of Ng^, sv pseudoparticles, and a total 
number of Ncu cv pseudoparticles. 

These here introduced quantum objects cannot, however, account for the en- 
tire Hilbert space, as all of them are ?7-spin zero and spin zero objects (remember 
that we always combine an equal number of " + |" and " — |" to form one av 
pseudoparticle). Since we have electronic densities na and magnetization ma in 
the ranges < na < 1 and < ma < na, this means, for example in the spinon 
case, that some spins would be "left out" in the pairing process of constructing 
the su pseudoparticles (because in general N-^ > N^). The physical properties 
of the holons and spinons that were "not able to make it" into any of the au 
pairs are quite different from those of the pseudoparticles, as further discussed in 
section ^TT^ . 

The description of pseudoparticles in terms of holons and spinons, rotated 
electrons and electrons, respectively, is given in Fig. !^2.1\i . 
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Figure 2.1: Flow chart describing how the electrons and the electron holes, due 
to the rotation transformation given by V{U /t), are described in terms of rotated 
electrons and rotated electron holes. These, in turn, are then forming the holons 
and spinons. These holons and spinons combine in a way described in section 
()2.2.3|) . to form the pseudoparticles. 
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2.2.4 Pseudoparticles, Yang holons and HL spinons 

Starting from the total number of electrons = Yla ^o-' where is the to- 
tal number of cr-spin electrons, we can derive expressions for the corresponding 
numbers of the pseudoparticles Nco, Ng^ and No, {Nau gives the number of av 
pseudoparticles: Nc2 = 3 means that we have 3 c2-pseudoparticles consisting of 4 
holons each, 2 with r^-spin projection +| and 2 with r^-spin projection — |). We 
will see that, in order to account for all the electrons in the system, we need to de- 
fine some new objects that are inherently linked to the pseudoparticle description, 
but that do not contribute to the dynamics of the Hubbard model. 

Our starting point is the symmetry consideration that led us to conclude that 
it is sufficient to study the region of the parameter space where < na < 1 and 
< ma < na. This means that: 



+ 1 - _ 1 =Na-N >Q 



M^^^.-M^^_.=N^-N^>0 (2.33) 



Since the total number of ±| holons (a = c) and ±| spinons (a = s) that take 
part in the au pseudoparticles can be written as Yl'^=i ^^au , we can write the 
difference between the total number of holons and spinons and those paired up 
in au pseudoparticles, as: 

oo 

^a,±i=M„,±i -^z/iV,, (2.34) 
v=l 

These remaining particles, there are a number of of them, have been 

given the names ±| Yang holons (for a = c) and ±| HL spinons (for a = s) [5^ . 
[9T] . "Yang" stands for C.N. Yang who authored Ref. 55 , whilst "HL" stands for 
Heilmann and Lieb, who authored Ref |93J. It is not difficult to hint a relationship 
between the total 77-spin and spin on the one hand, and the total number of Yang 
holons and HL spinons on the other. In this case of a LWS, we have L^^f = 

Na-N = 2Sc and L^^f = in the holon case, and Lj^^f = A^^ - A^^ = 2^, and 
l^Lws = Q in f^i^Q spinon case. Now, by acting onto the LWS by S}^, we increase 

S, 2 

the number of _i by one, on the expense of which decreases by one (the 

z-component of the ?7-spin and/or spin changes accordingly). Like this we can 
continue until we reach the HWS. Hence L^ _^i = Sa — 2(±|)S'^ = 0, 1, ... , 25'^ 
and furthermore La = L^_^_i + L^ _ 1 = 25*0,, where Sa is the total ?7-spin (a = c) 
or spin (a = s) of the system. 

These HL spinons and Yang holons behave quite differently inside the many- 
body system, than the cO-pseudoparticles and the au pseudoparticles. Firstly, the 
creation and annihilation operators of the former objects do not commute with 
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the generators of the ?7-spin and the spin algebras (as aheady noted above), whilst 
the corresponding operators of the latter objects indeed do. This is easily seen 
by understanding that the cO-pseudoparticles have no r^-spin or spin degrees of 
freedom, and that the pseudoparticles are singlet ?7-spin (a=c) and spin (a = 
s) quantum objects, and thus yield zero when acted upon by any of the ?7-spin and 
spin generators, respectively (this is easily confirmed mathematically by forming 
a singlet state of rotated electrons corresponding to a certain pseudoparticle, and 
then letting any of the ?7-spin or spin generators act upon this state). 

The physical implication is that all {28^ + 1) energy eigenstates inside any 
given ladder have the same occupancies of pseudoparticles: "climbing up" the 
ladder, from a LWS towards a HWS, will only change the numbers L^^^i. Due to 
the fact that the occupancies of pseudoparticles are insensitive to the application 
of the generators of the two SU{2) algebras, the discrete momentum values and 
hence their corresponding rapidity numbers stay unaltered as well. If we define 
the number operators Nau = J2q-^au{(l), with eigenvalues Nau, we have that 



a' = c,s (2.35) 



valid for all au branches. Next, we claim that all the six generators of the two 
SU{2) algebras commute with the unitary operator V{U/t), whilst the operators 
Nauil) do not. Actually, the last claim is trivial: if they would commute, then the 
pseudoparticles as described by the rotated electrons would be the same objects 
as described by the original electrons. However, since double occupancy is not a 
good quantum number for the original electrons, we know that the pseudoparticles 
can not be described as simple occupancy configurations of these objects. For the 
generators of the two SU (2) algebras, we have the opposite: these operators create 
the same quantum objects in the unrotated as well as in the rotated frame. Unlike 
the operators for the pseudoparticles, the generators of the two SU{2) algebras 
can be easily expressed both in terms of electronic as well as rotated electronic 
creation and annihilation operators: 

Si = E(-l)^^l^^T = V{U/tyslV{U/t) = $^(-l)^gliQT (2.36) 

i i 

where 5"^ serves as an example of a typical SU{2) generator. The point is that 
Sl has the same expression in terms of electrons and rotated electrons. This was 
studied in detail in Ref. [HI], where the electronic double occupancy expectation 
value Del (q) was obtained as a function of the occupancies of various pseudopar- 
ticles. It was found, for example, that when creating a cl-pseudoparticle, the 
electronic average double occupancy did not in general increase by one (which 
would be the naive guess), but rather, it was found to depend on the momenta q 
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and on the value of (U/t), according to 



DcM < D,i{q) < D,i{n - 2kF) = 1 



In Fig. 4 of that reference, we see that -Dci(O) ^ 1 as {U/t) —>■ oo. This 
can actually serve as a " measure" for how the electron - rotated electron unitary 
transformation depends on (U/t); the closer -Dci(O) gets to 1, the closer V{U/t) 
gets to unity. The fact that the generators of the two SU{2) algebras commute 
with V{U /t) implies that the Yang holons and HL spinous are the same quantum 
objects in terms of electrons as they are in terms of rotated electrons, due to the 
fact that the creation and annihilation operators for these objects are nothing 
but the off-diagonal generators themselves, multiplied by a normalizing constant. 
Indeed, it was found in Ref. jHl], that creating a — | Yang holon always yields 
a double occupancy increase by one, D ^ D + 1, independently of {U/t). Fur- 
thermore it was found that creation of a Yang holon or a HL spinon does not 
change the expectation value of the kinetic operator, thus deeming these objects 
"localized" in terms of (lack of) charge and spin transport. Summarized: 



v{u/t),si = V{U/t),S^ = V{U/t),S, 



which implies that 







a 



c, s 



(2.37) 



V{U/t),L 



V{U/t),L 







(2.38) 



where L^^_^i is the ±| Yang holon {a = c) and ±| HL spinon {a = s) number 
operator, respectively. 

Finally, the proof of the consistency of this pseudoparticle, Yang holon and HL 
spinon picture, with the Bethe-Ansatz solution, in terms of the counting of the 
states, was given in Ref. [221 • ^^^^ reference, it is shown that this representation 
accounts for all the 4^" eigenstates of the model, as well as showing that the 
number of states in a subspace consisting of a fixed number of Mq, = ^a,a 
holons or spinous with a fixed r^-spin or spin value S^, equals the number of states 
obtained by au pseudoparticle occupancy configurations (where au = cO, cz/, su), 
according to the previously described scheme of the recombination of rotated 
electrons into pseudoparticles. 
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2.2.5 Occupational configurations of the pseudoparticles 

Since the pseudoparticles are derived from the Takahashi string hypothesis, with 
the occupational configurations given by the quantum numbers, we can define 
the functions Nco{qj), Nsu{qj) and Ncy{qj) respectively, where Nay{qj) = 1 means 
that the discrete pseudoparticle momentum qj is occupied and Nauiqj) = means 
that the discrete momentum qj is unoccupied. To each set of allowed quantum 
numbers correspond a unique set of rapidity numbers, which in turn correspond to 
a unique energy eigenstate [00] [HI] ■ The pseudoparticles obey a generalized Pauli 
principle, known as Haldane statistics |lU8j . In short, the Haldane particles affect 
the number of states available to any other particle in the many-body system 
(whilst in the case of exact fermions, only the state occupied by one fermion is 
forbidden to the next fermion). 

Let us define pseudoparticle creation and annihilation operators: h\^^^ and 
hq au respectively, where q is the momentum and av = cO, cz/, su. The pseudopar- 
ticle number operator can be written as: 

Nauiq) = (2.39) 

which when summing over all momenta gives the total number of au pseudopar- 
ticles 

Nau = J2^-u{q) (2.40) 

q 

The statistics obeyed by the pseudoparticles is 



ir) I ,p«""^(9'-'?)/2 



iV*,sin(^-(pO)]- Ql^ = ±n 

(2.41) 

where N^^ is defined below and Q^,^ measures the quantum shake-up effect, and 
is introduced in section ()2.3.2|1 . Q^^^ is zero for the ground state (by construction) 
and nonzero if the actual excited state is described by quantum numbers shifted 
as compared to those of the ground state (from integers to half-odd integers 
or vice- versa). The momentum dependent creation and annihilation operators 
are formally defined locally on an effective av lattice, with lattice constant aav, 
defined so that the length of such a lattice is av independent and equal to L: 

= «^ (2-42) 

au 

where A*"*^ is the number of av lattice sites. That in general N^^ ^ Na stems 
originally from the upper and lower bounds on the quantum numbers from the 
Takahashi string hypothesis ()2.1U|) , which controls the value of the largest possible 
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occupied momentum (the positive boundary of the effective Brillouin zone, the 
largest quantum number) and the smallest possible occupied momentum (the 
negative boundary of the effective Brillouin zone, the smallest quantum number). 
Only in one case (for the cO-pseudoparticles) , does the total number of allowed 
discrete momenta equal the number of "real" lattice sites Na- The sum of the 
number of occupied momentum values and the number of unoccupied momentum 
values, must always equal the number of au effective lattice sites. 

= iV„, + (2.43) 

This gives us instantly that A^*q = Na due to Eq. ()2.32|) . Ref. |^ provides us 
with an expression for the total number of au pseudoparticle holes: 

oo 

N^, = L^ + 2 J2 (2.44) 
which also can be expressed as: 

oo 

Nl = iV, - iV,o - 5^ (z/ + z/' - |z/ - z/'l) iV,,, v>l (2.45) 

oo 

Nt = N,o - ^ (z. + z/' - |z/ - z/'l) Ns., V > 1 



u'=l 



The number of allowed momentum values in the ground state reads: 





= Na 








= Na 


-N 






= N^ 






su 




-N, 


v>2 



(2.46) 



(the corresponding numbers for an arbitrary excited energy eigenstate will be 
given in section ()3.1.6j) ). With these numbers well defined, we can then relate 
the local av pseudoparticle operators with the Fourier transformed momentum 
av pseudoparticle operators: 



TV' 

1 ge-^«^-^6,,„. (2.47) 
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where fej^,^ creates a au pseudoparticle on lattice site position j (with space 
coordinate xj = jaau) and bj^au annihilates a au pseudoparticle on lattice site 
position j. From the equalities of Eq. ()2.30|) . we have that Nau commutes with 
the original hamiltonian. 

As can be deduced from Eqs. fl2.44|l and ()2.4(ij) (together with the total number 
of Yang holons and HL spinous), we find that the ground state is completely void 
of siy pseudoparticles for z/ = 2, 3, . . . and of cz/ pseudoparticles for all = 1, 2, . . . 
|82j . Moreover, the two branches that have finite occupancies in the ground state, 
as given by Nco{q) = N^oiQ) and Nsi{q) = N^i{q), are both densely packed around 
a minimum energy - zero momentum point with well defined left and right Fermi 
points. In the finite ground state system, this means that the occupied quantum 
numbers are symmetrically distributed around zero, with the exception of (l/L) 
corrections (see below). We can define the smallest possible quantum number 
for the au branch as /"'^ and the largest as J"*^, defining the limiting momenta 
for the effective Brillouin zone, and similarily for the occupied momenta, the 
negative (left) Fermi point and the positive (right) Fermi point Ip'^. The 
"true" limiting momenta for the effective Brillouin zone and Fermi momenta will 
be shown to be (l/L) corrections to the momenta given by these numbers. In our 
new notation, we can reformulate the conditions on when the quantum numbers of 
the Takahashi string hypothesis must be integers and when they must be half-odd 
integers. For all au ^ cO branches, we have that 

A^'*^, even ^> J"'^ half-odd integer 
N*^ odd =^ If" integer 

whilst for the cO branch, 

— ^ — Nay even ^> if* half-odd integer 

au=cv^sv 

-Y~ ^""^ ^ integer (2.48) 

ai'=cu,su 

For all au branches, let us define the quantum numbers introduced above accord- 
ing to 

N* -1 N* -1 

- ~ 2 +~"~2 

- 1 A^ - 1 

= --^^^ /r+ = = '-^^^ (2.49) 

To obtain the limiting momenta for the effective Brillouin zone as well as the 
Fermi momenta, let us define some temporary variables Kp^u '^Fa^ 
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respectively, according to: 



(2.50) 



Next, let us consider the ground state configuration. To make matters simple, 
we will only consider LWS ground states, such that L^_i =0 which translates 
into Nsi = Ni since there are no other j-spins in the system. Furthermore, we 
have that Ncu = N^^, = for all z/ = 1, 2, ... in the case of the cu branches and 
u = 2, 3, . . . in the case of the su branches. The only relevant k"^^ and k^q^, 
becomes 
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(2.51) 



Fsl 



where we have used Eqs. ()3.28|1 and (j2.49|l together with A^^i = Ni. The differ- 
ent cases of limiting momenta for the effective Brillouin zone q^^, and of Fermi 
momenta q^^^ are, according to Ref. jHS], in the au = cO case: 



aiy = cO Limiting momenta for the effective Brillouin zone: 

^« \^ AT ^ ± ± 

— - 2^ N^u even =^ Qco = ^^co 



2 



2 



^^'^ ^ ^c^o = ± 7 



2 

ai>=cu,su 



au = cO Fermi momenta (LWS): 

Nau and both odd or both even ^> qp^Q = Kp^Q 



av=cv,sv 



~ X] odd and N even =^ qp^^ = - ^ 

au=cu,si' 

^ - ^ Nau even and N odd q$^Q = Kp^Q + ^ 

au=cu,su 



40 



and in the au = si case: 



aiy = si Limiting momenta for the effective Brillouin zone: 

N*^ = N-f even or odd (i.e. always) =^ qf^ = 



si 



av = si Fermi momenta (LWS): 

A^*^ = N-^ and A^^i = both even or both odd ^> qp^^ = 



Fsl 

N*^ = even and A^",! = odd =^ q$^^ = k$^^ + y 



iV*^ = A^-f odd and Nsi = Ni even =^ q^^^ = Kp^^ 

where the k^j^'s and the K^au^ defined in Eq. ()2.51|) . 

So far everything has been described for the finite system, however, later we 
will frequently use the same quantities in the thermodynamic limit. In this limit, 
the notation becomes less heavy due to the fact that we neglect {1/L) terms. 
Let us define = lim^^oo Q'^^ and qpaiy = ^it^l^oo qpau- The ground state 
occupancy configurations of the pseudoparticles become: 



NM = SiqFcO-\q\) \q\<q% 

NM = OiqFsi-\q\) \q\<q'si (2.52) 
Nl{q) = \q\<ql 

where here ai' stands for all other pseudoparticles 7^ cO, si and the Fermi mo- 
menta qpco and qpsi and the limiting momentum values of the effective Brillouin 
zone g°Q, q'^^ and respectively, become: 

qFsi = kFi = ^FT (2.53) 

and 

q^su = kF^-kFi z/ = 2,3,... 

g°^ = --2A;i. = 7r(--n) z/=l,2,... (2.54) 
a a 

where kp = {kp^^ + kp\^) /2 = {Tm-^ + 7mi)/2 = Tcn/2 is the usual Fermi momentum. 
Note that we cannot define a corresponding Fermi momentum for the cv and the 
su (z/ = 2, 3, . . .) bands due to the absence of these pseudoparticles in the ground 
state. However, the well defined Fermi points and effective Brillouin zones of the 
cO and si momenta, allows us to define some typical ranges and values for the 
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ground state rapidities: 



|g| < 2kF =^ \k°{q)\ < k\2kF) = Q 

\q\ < kp^ =^ \Al,iq)\ < AO,(A;pJ = B (2.55) 

where actually k°{-2kF) = ~k\2kF) = -Q and kl^{-kFi) = -A°i(A;^i) = -B, 
defining the quantities Q and B. The occupancy configuration functions in mo- 
mentum space will be important when deriving expressions for the pseudoparticle 
energy bands and phase shifts, amongst other quantities. Some well known lim- 
iting values of these entities include: 



(z/ > 2) (2.56) 

where the two last limits are particularly interesting: the actual bands vanish and 
thus the entire dynamics of the system is described by the cO and the si bands. 
The limit na 1 is thus the limit where the Fourier momentum space of the cu 
{v > 1) pseudoparticles disappears, just as the limit ma — > is where the Fourier 
momentum space of the sv {v > 2) pseudoparticles disappears. Later, we will see 
that due to this effect, the Fermi points of the cO and the si bands, are in the cv 
and sv cases simulated by the limiting values of the effective Brillouin zone, equal 
to being the only momenta points to survive in these limits. Half filling is a 
limit that we will avoid in the dynamical theory of chapter however the zero 
magnetization limit will indeed interest us. 



k{q) —>■ q U/t oo 

B -^oo ma — » 

g° ^ na ^ 1 
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2.2.6 Energy and momentum deviations (Introduction) 

In order to formulate a dynamical theory, which is necessary in order to obtain in- 
formation about the spectral properties of the model, we need to find expressions 
for the deviations (from the ground state) in energy, associated with pseudoparti- 
cle excitations in the many-body system, in terms of pseudoparticle energy bands 
and pseudoparticle number deviations. Even though the proper representation 
for the dynamical theory, introduced in section ()2.3.1|) . is different from that of 
the pseudoparticle representation, it is similar enough to allow for the subsequent 
study. We will see that the quantities derived here will lead us naturally to the 
new representation of section ()2.3.1|) . 

In the following, when we talk about "transitions" to an excited state, we 
mean that the occupancies of the Takahashi quantum numbers go from the ground 
state distribution (which is a densely packed distribution of numbers around a 
minimum energy point) to some other distribution. A particle-hole excitation 
of a pseudoparticle means that the occupied quantum number in the ground 
state, becomes unoccupied at the expense of some other quantum number that 
in turn becomes occupied. On the other hand, when adding pseudoparticles, or 
when removing pseudoparticles, there will be a net increase, or decrease, in the 
number of occupied numbers, leading to a new excited state configuration. We will 
only consider excited states that differ from the ground state in the occupancy 
of a small number of pseudoparticles, even though the formal requirement is 
much more general: the number of excited electrons must remain finite [LSTlj . 
This means that, when evaluating correlation functions, we only allow operators 
whose expressions involve a finite number of electronic creation and annihilation 
operators. 

The pseudoparticle number deviations, depending on the momentum, will be 
a key quantity is this analysis, since all other quantities will ultimately depend 
on the occupation and non-occupation of pseudoparticles for different momenta. 
To obtain the lowest order corrections in energy, we should consider "small" 
deviations from the ground state pseudoparticle number configurations. This will 
be enough since, as studies of following sections will confirm, an overwhelmingly 
large portion of the total spectral weight of the one-electron addition and removal 
processes are generated by excitation of only a few pseudoparticles. Since the 
Takahashi string hypothesis is only valid when A^^ ^ 1, we often consider the 
continuous momentum limit such that g^+i — qj = (27r/L) 0. This means that 
we can stop talking about discrete quantum numbers altogether, and replace the 
sums by integrals in the Eqs ()2.7|) . We thus arrive to the "continuous momentum 
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Takahashi equations": 



q y / (ig' A^sjy(g') arctan 

^ J-qO 



j dq' Ncy{q') arctan 



yco 



sinA:(g)a - Kyj^q') 
/ sin k{q)a - KM') 



1 / 

kcvio) = Q'H — / c^?' Nco{q') arctan 



vr ./_„o 

9c0 



Acu{q) — sin k{q')a 



1 ~ /-g" , 

l" dq' N,,iq') 



VU 

KM)-Ky,{q') 
u 



+ 



(2.57) 



= g- - 



^ oo 



c/g' Nco{q') arctan 



As^iq) - sin k{q'] 



dq' Nsu'{q) ^uu' 



uu 

AsM - ^su'jq'Y 
u 



+ 



where kcu{q) = ^Re {arcsin(Aci/(g) — ivu)} and u = U/At. The function 9yjy/(x) 
can be found in Eq. ()2.9p . 

The energy and momentum can according to the Takahashi string hypothesis, 
Eq. ()2.11|1 . easily be re-expressed in the continuous hmit: 



E = U— 

271 



u 
27r 



27r 



1 /•^"o 
2 



-1% 



I dq Nco{q) cos k{q)a + > / dq Ncu{q)Re iyl - {Acu{q) - i^u) \ 
dq Nco{q) -u dq uNcu{q) 

J-pO 

dq N,o{q)k{q) + > / dq N,,{q) - - k,,{q) 



a ' 2 



/ " iV,o(g)g + / iV,,(g) (^^ _ + ^ rfg iV,,(g)g 



a 



The two equivalent expressions for the momenta are obtained by using Eqs. 
()2.57|) . The constant term ^M^_i is obtained by using Eq. ()2.8p together with 
Eq. ()2.12j) . and shows the constant momentum value of the — | holons.The occu- 
pancy functions Nau{q) have a well defined value for each energy eigenstate. 

What we will do next is to allow a small deviation ANau{q) to perturb the 
ground state occupancy configurations, 



Nau (q) = A^l (g) + A AT,, (g) au = cO, cu, su 



(2.59) 
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where the N^i,{q) are given in Eq. ()2.52jl . Eq. ()2.59jl then describes an excited 
energy eigenstate. 



2.2.7 Energy deviations and the ^ai^^a'i^' functions 

Since it is the deviations from the ground state that interests us, we would hke 
to shape our theory so that all quantities are expressed relative to this ground 
state. This is because all dynamical quantities in the following will depend on 
ANau and not on Nau- Therefore, in order to capture the relevant dynamics of 
the problem, we will formulate all quantities in a " normal ordered relative to the 
ground state" fashion. 

Our starting point is to express the energy of the excited state of the system 
as a ground state energy plus higher order corrections due to the introduction of 
the pseudoparticle deviations Eq. ()2.59j) : 

oo 

E = ^E(J) (2.60) 

j=0 

The energy deviation E^^^ will be expressed as proportional to some pseu- 
doparticle energy band relative to the ground state, multiplied by the corre- 
sponding pseudoparticle occupancy first order deviation. Therefore, we define 
the pseudoparticle energy bands e°j,(g) such that 

^^'^ = ^ E E I'^dq elM^N^q) (2.61) 

a=C,S u=Oa,B 

The relative to the ground state pseudoparticle energy bands e'^j^^q) are de- 
fined as the functional derivative of the energy with respect to the occupancy 
configuration deviation: 

The second term in the expansion of Eq. ()2.6()j) would contain bilinear com- 
binations of the ANau^s and would correspond to residual energy interactions 
between the different au pseudoparticles. 

Now, a nonzero deviation in the au occupancy configuration yields, as can be 
seen in Eq. ()2.57|) . a corresponding deviation in the rapidity functions: 

= k^{q) + Ak{q) 

NaM = Nliq) + AN^M) =^ { KM) = Aliq) + AA^q) (2.63) 

-Al{q) + AAUq) 
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where au = cO, cu, su. The ground state rapidity functions come from the solution 
of the Takahashi equations for the particular case of having the occupied quantum 
numbers in their ground state configurations, i.e. such that the momenta obeys 
Eq. (I2.52|l . They are obtainable by solving ()2.57|1 with Nay{q) = N^^{q). Before 
deriving the energy bands, however, we will investigate the continuous momentum 
Takahashi equations a little bit further, in order to derive relationships between 
a new set of functions (denoted ^au,a'u'{'l, q') below) that will be used to express 
the energy bands in an elegant way. 

By introducing the pseudoparticle occupational deviations together with the 
rapidity deviations (both of them are given in ()2.(jH|l ). together with Eq. ()2.52j] . 
into the Takahashi equations ()2.57p . we obtain equations separable order by order. 
Focusing on the zeroth order contributions (the ground state) and the first order 
deviations, we note that as the algebra turns out, we can simplify matters a lot by 
applying d/dq to the zeroth order equations and inserting them into the first order 
equations. To simplify matters even further, we define a new quantity AQau{q) 
according to: 

Ak{q) = ^^AQM 
(iA° (a) 

AKu{q) = — ^AQ,,(g) a = c,s (2.64) 

which allows the first order contributions of the Takahashi equations be written 
as: 

AQ.o(g) = — r dA ^^^^^^^ + E E r'^' ^co,a'AQ, q')^N^'u'{q') 



AQM = -^l'dk ^Q^oik) cos ka ^ ^ ^ ^ q')ANMl') 



I I ( {(l) -sin ka \ ^ 
y uu J 



TlVUj^n , I At.(a)-sinka\ ^. — ' J _„0 



AQ,,{k)coska 1 ^.AA^.^^^/ (Kuiq)-^ 



AQsuiq) = / dk ^^^^ -2 dA AQ,i(A)e',i + 

TiyuJ_Q ^ ^ f AUQ)-sinka Y 2tXU J _b \ U 



0° 



+ E E / „ rfg'2;,,,,v(g,g')AiV„v(g') (2.65) 



a'=C,S u'=6^i ^ '^a'v 



where we have introduced the function AQay^X) = AQau{X{q)) = AQau^q) and 
the functions Zau,a'u'{q,q') which can be obtained explicitly and are given below. 
0'(x) = dQ{x)/dx from Eq. ()2.9|) and Q and B are defined by Eq. ()2.55p . The 
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equations ()2.65|1 express relationships between the rapidity deviations, since 



{dk''(q)/dq] 



by definition. Examining the mathematical form of the relationships between the 
different AQ„,^(g)'s, we see that there is a possibility of expressing them as linear 
combinations of the pseudoparticle deviations, 

^Qauiq) = J2 '^"-."'-'(^' 9')AiV,v'(g') (2.67) 

or equivalently, 

if we take the functional derivative with respect to these deviations, i.e. applying 

6/6ANa'u' to the equations ()2.65j) . By doing this, and by changing variables 
according to 



k{q)a ka = r or r' 

A„i.(g) Aau ^ = r or r' 



$a^,aV' (g, q') ^a^,a'^' (r, v') (2.69) 



where the indices of $ always indicate what au branch the variable belongs to 
(the first variable is always unprimed and the second always primed), we find 
after some algebra: 

1 /--s/" $ ^ , Jj," 
$co,aV' (r, r') = - / dr" — + ^^CaV (r, r') 



TT 



1 r"^"/" . $.n.',/fr". 



$...V'(r, r') = - / dr" 7;';';:,,;. ^ " (2-70) 

rfr" <l,i,„v'(r",r')e:,i(r - r") + z..,„v'(r,r') 



1 

'2^ 



where the auxiliary functions Zav^a'v' are defined by: 



47 



ZcoM^y = Zco,cu{r, r') = -(j)y{r - r') Zco^suij, r') = -(pyij- 

Zcu,co{r, r') = (j)^{r - r') Zcu,cu'{r, r') = Quu'{r - r')/2n Zcu,su'ir, r') = 
Zs,y,coi'^^ = -(P^r - r') Zcu,cu'ir, r') = Zs,,,s,y'ir, r') = Qyy'{r 

Here (pui.^) = arctan(a;/z/)/7r and Zay^a'u'{.f^f') = Zau,a'u'{(l,(l')- 

We note that the functions defined in ()2.7()|1 obey the following symmetry: 
^au,a'u'{f^, f') = —^au,a'v'{.~f, ~'^') i which together with the oddness of the ground 
state rapidity functions implies that ^au^a'u'il^ l') = —^au,a'u'{.—<l, —q')- The same 
line of thought can be applied to the energy, by use of the energy expression ()2.58j) . 
introducing the pseudoparticle deviations and separating contributions order by 
order. The newly derived ^ai',a'u'{r,r') enters the calculation via Ak{q), and we 
can by comparing the different resulting first order terms with Eq. (j2.6Hl . obtain 
the pseudoparticle energy bands: 

/■Q _ U 

€%{q) = -2t cos k°{q) a + 2ta / dk sin A; $co,co(fc, (<?))- tt 

J~Q ' 2 

e°^(g) = 4tRe | ^1 - (A°^(g) - ii^u^^ + 2ta J dk sin ka $co,c^(fc, A°^(g)) - uU 

eliq) = 2ta / dk sinfca ^,oAk,KM) (2-71) 
J-Q 

where we have used an alternative to ^au,a'u'iq,q') = ^au,a'u'{k^{q), -^auil))- 
a bi-product of this calculation, we have that the zero order energy term, i.e. the 
ground state energy, can be expressed as: 

E^^^ = -2t— dq cos k\q)a (2.72) 

27r J^2kp 

where the ground state rapidity function k^{q) satisfies the first equality of Eq. 
()2.57|) . with Nau{q') = N^uil') (^^^ ground state configuration). 

We want, according to convention, to fix the reference levels of these bands so 
that the cO and the si bands gives zero at their respective Fermi points. The cu 
(z/ > 1) and the su {v > 2) bands will then have their reference levels adjusted 
according to their u = {a = c) and u = 1 {a = s) counterparts. This adjustment 
is a consequence of breaking one or both of the SU{2) symmetries of the model. 
When this happens, the energies will depend on the chemical potential and the 
magnetic field strength, respectively. For cu {v > 1) this energy difference is 
proportional to the number of doubly occupied rotated electron sites belonging 
to the cu pseudoparticle, whilst for su {v > 2) it is proportional to the number 
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of |-spin singly occupied rotated electron sites belonging to su. Hence both of 
these contributions are equal to u. We define thus the following energy bands for 
the pseudoparticles: 



eco(g) = e°o(g)-e°o(2M 



(2.73) 



where /ic = 2/i and /i^ = 2fiQh were defined in section fl2.1.H) . The energy bands 
eco(q'), esi(g) and e^^{q) (for au ^ cO, si) are even functions of q, and are such 
that: 



where z/ > 1 and u' > 2. The mathematical exercise of deriving the ^au,a'u'{q, q') 
functions payed off, judging by the beauty of the derived energy band expressions. 
We use Eq. (12.7111 to numerically obtain the dispersion relations shown in Figs. 
()2.2|) , ()2.3|) and ()2.4|) , for a " very large" , " intermediate" and " very small" value of 
{U/t). The filling dependence on these relations is discussed separately in section 
()5.4j) (for a further discussion on these energy bands, see Ref. ^|). In order to 
obtain these dispersions, we need (for example) the functions $co,cO; "^co.si, ^si,co 
and These are plotted in Figs. ()2.5|) . ()2.6|) and ()2.7|) . respectively, for three 

different values of {U/t), and are further discussed in Ref. [HHl- 

This analysis can be carried to higher orders, where the terms of order j in- 
cludes i factors of different AiV„,^(g)'s. By keeping 2nd order terms, including 
products of type [/S.Nau{,q)^Naiu'{,q')], we can derive the residual energy inter- 
action term E^'^'^ between the pseudoparticles. That this term is finite shows 
that the pseudoparticles have residual energy interactions. By following the same 
general scheme as for the first order (presented here above), we arrive to 



after a considerable amount of algebra. Here we have introduced the pseudopar- 
ticle group velocity Vau{q) and the pseudoparticle Fermi velocity Vau, defined by: 



eeo(2M = esi(fcn) = - 2M 



eL'(^FT - kpi) = 



(2.74) 




(2.75) 




V, 



(2.76) 
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Figure 2.2: Energy dispersion relation for the au = cO pseudoparticle, for (U/t) = 
0.3, 4.9, and 100 and for n = 0.59 and m 0, in units of t. The (U/t) = 4.9 
curve is visible in between the other two curves. Note that the dispersion for 
—2kF < q < 2kF becomes successively deeper for (U/t) oo, however always 
keeping the bandwidth constant at 4t. 

That -E*-^-* only contains terms of order [ANau{Q)^Na'u'{Q')] can be seen by in- 
serting Eq. ()2.67|) into Eq. ()2.75p . We note that the last term of E'^'^^ is actually 
of order (l/L) due to the square of AQauiJlPau)- 
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Figure 2.3: Energy dispersion relation for the au = si pseudoparticle, for (U/t) = 
0.3, 4.9, and 100 and for n = 0.59 and m — > 0, in units of t. Note that the energy 
bandwidth is a decreasing function of (U/t) and that si pseudoparticle becomes 
dispersionless in the (U/t) — oo limit. 



t"ci(q) 
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\u=0.3t 


/ 1 














U=100t\\ 



-(:7r-2kF) (;r-2kF) 

Figure 2.4: Energy dispersion relation for the av = cl pseudoparticle, for (U/t) = 
0.3, 4.9, and 100 and for n = 0.59 and m ^ 0, in units of t. The energy 
bandwidth is a decreasing function of (U/t) and the cl pseudoparticle becomes 
dispersionless in the (U/t) oo limit, and is a horizontal line along the zero 
energy level in the figure. 
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Figure 2.5: The functions ^au,a'v'{'i,<i') for {U/t) = 100, n = 0.59 and m — > 
0, arranged according to: ^co,co{<l,<l') (upper left), $co,si(?, (upper right), 
^siMl^l') (fower left) and $si,si(g,g') (lower right). 
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Figure 2.6: The functions ^au,a'u'{<l,(l') for (U/t) = 4.9, n = 0.59 and m — 
0, arranged according to: ^co,co{<l,(l') (upper left), ^co,sii.1yl') (upper right), 
^si,co{l^l') (lower left) and $si,si(g, g') (lower right). 
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Figure 2.7: The functions ^au,a'v'{,<i,<i') for {U/t) = 0.3, n = 0.59 and m 
0, arranged according to: ^cq,co{<1,<i') (upper left), ^co,si{q,<i') (upper right), 
^siMl^l') (fower left) and $si,si(g,g') (lower right). 
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2.2.8 Momentum deviations and inverse rapidities 



According to the two equivalent expressions for the value of the total momenta 
P given in Eq. ()2.58|1 . we can in the same way as with the energies introduce 
the pseudoparticle number deviations into these two expressions. Since we have 
one expression for P involving only the momenta {Pmom) and another involving 
only the rapidities (Prap), we hope that by equalling the two expressions, it will 
be possible to deduce some interesting relationships between the pseudoparticle 
momentum and the ground state rapidities. In the same "normal ordered relative 
to the ground state" spirit as before, let now 



P = J]pO) where = P^^ = p(^) (2.77) 

j=0 

i.e. we will equal the terms order by order in the deviations. By reasons that 
will become obvious in section ()2.3.1|1 . we will focus on the j = and the j = 1 
terms. By using Eqs. ()2.fifi|l and ()2.fi7j) . we arrive to the following relationship: 

dq Ak{q) = / dq AN,o{q) / dk $cO,co(fc, k{q)) + (2.78) 

-2kF J^qO, J-Q 

rQ 



+ V V / dqAN^^iq) rffc l',o,a^(fc,A„^(g)) 
^-aS„ J-Q 



We find that by comparing the first order deviations with each other, the fol- 
lowing relationships between the pseudoparticle momenta and their corresponding 
rapidities can be derived: 



q = k\q)^ r dk^,o,co{k,k'iq)) 
J-Q 



q=klM- I dk^,oAk,^liq)) (2.79) 

-Q 
-Q 



q= dk^^oAK^liq)) (2.80) 
J-Q 



From these relationships it is possible, at least in principle, to obtain the explicit 
dependencies of the ground state rapidities on the momenta. However, there is 
more to be done than just to make an abstract mathematical claim. By defining 
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the inverse of the ground state rapidities as functions of momenta according to 

rQ 

qco{k) = k+ / dk' ^co,co{k',k) 

J-Q 

rQ 

qsu{A)= dk'^^oAk'A) (2.81) 

J-Q 

2 ^ 
qcu{^) = -3^ {arcsin (A - iuu)} - / dk' $co cu{k', A) 
« J-Q 

we can by taking derivatives define the following density functions: 
dqco{k) 



2Tcpcoik) 



dk 



1+ / dk' ^^^,^^,{k',k) 

J —Q 



27ra,,(A) = = ^ dk' ^^co^k', A) (2.82) 

2vr^c.(A) = %^ = -Re 1 ^ i - ^ dk' ^T^coAk', A) 

^ [ y/f - (A - zuuf J -^-Q 



Note that these functions can be written as 

dk\q) _ 1 
dq 27ipco{q) 
dALiq) 1 



dq 2Tiaay{q) 



au = cu, sjy (2.83) 



by inverting the equalities of Eq. ()2.82|) . where 2Trpco{q) and 2TTcrau{q) are the cor- 
responding pseudoparticle momentum dependent density functions. By inserting 
the former into the first derivative of the zeroth order "deviations" of the contin- 
uous Takahashi equations, we can obtain coupled integral equations in terms of 
these density functions. 
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2.3 Pseudofermions 



2.3.1 Introduction 

In this section we will introduce some new quantum objects, specify some of their 
properties and relate them to the pseudoparticles. The physical interpretation of 
some derived quantities in terms of relevant Hilbert spaces for our ground state 
— > excited state transitions, and in terms of phase shifts and scatterers, will be 
presented in section ()2.3.2|) . 

Eq. ()2.75|1 shows that the pseudoparticle residual energy interaction contains 
the quantities AQau{q), for au = cO, cz/, su. These quantities were, in turn, intro- 
duced to facilitate the mathematics, but they also carry an important physical 
meaning. In the following, we will incorporate these quantities in the definition of 
a quantum object related to the pseudoparticle, namely the pseudofermion. The 
pseudofermion will not have any residual energy interaction terms, and will allow 
us to construct a dynamical theory for the ID Hubbard model. In the following, 
we will illustrate the "birth" of the pseudofermion picture by using the example 
of the cO pseudoparticle, but is obviously valid for any pseudoparticle branch (by 
letting k^{q) Aai,{q) whenever au ^ cO). 

We have 

k{q) = k\q) + Afc(g) = k\q) + ^^Ag,o(g) (2.84) 

dq 

and equivalently for the other av branches. However, a normal Taylor expansion 
of k^{q) yields 

k%q + 5{q)) = k%q) + —^S{q) + ... (2.85) 

dq 

where \S{q)\ is a small number. 

By defining Q%{q) = LAQco{q), we see that the two expansions become equal 
if we define 6{q) as 

6{q) = AQUq) = (2.86) 

That S{q) is indeed of order (l/L) should be clear by observing that AQco{q) 
carries a factor of (27r/L) in the discrete system. The summation over the effec- 
tive Brillouin zone only contributes when the pseudoparticle deviation is nonzero, 
which only happens a finite number of times (for the cases that will become rel- 
evant when studying the one-electron spectral functions, we will see that this 
happens typically no more than two or three times, i.e. for two or three momen- 
tum values), which should leave no doubts concerning the order of magnitude of 
AQM- 
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Up to first order, we can thus write the following relationship between the 
excited state rapidity and the ground state rapidity 



L 



A^^{q) = A^^lq + 9lM\ a = c,s i/=l,2,... (2.87) 



This is quite remarkable because it states that the excited state rapidity can 
be expressed by the ground state rapidity, if we shift the momenta by an amount 
Qauil)/^- W^e see that all excited states that we are interested in can thus be 
expressed in terms of the ground state rapidities, if we are using a slightly shifted 
value for the discrete momenta. Also, all other properties of the pseudoparticles 
(for example their constitution in terms of rotated electrons, the fact that the cu 
pseudoparticles are ?7-spin zero objects and that the si> pseudoparticles are spin 
zero objects) remain intact. The "cruncher" of this new formulation is, however, 
that if we use these shifted discrete momentum values, the energy deviation ex- 
pansion corresponding to Eq. ()2.60|1 will yield zero or non physical expressions 
for all terms other than the first two: 

E = f2E^^^ = E^'>^ +E^'^ = Egs + AE for q ^ q + 9llM (2.88) 

j=0 

which is most easily understood by investigating the expressions for E^'^\ Eq. 
()2.75|) . since this quantity is proportional to Qauil) itself. Since explicit calcula- 
tion of the second order case is very lengthy with no contributions to the physical 
understanding, it seems more fruitful to present some simple reasons as to why 
this is true. 

We remind ourselves that Q^uil) ^ measure of the discrete momentum 
shift, due to the ground state — *• excited state transition. By letting the original 
momentum values include this shift already "from the start", we have that there 
is no extra shift in the momenta to use in the rapidity expansions, since this 
shift is already recorded by the momentum values q + Q%j^{q)/L. Hence, we 
should put Qauio) = in the expression for E'^'^\ Eq. ()2.75|) . which renders 
E^"^^ = exactly. Hence, the conclusion is that if we define a new set of quantum 
objects, with momentum values equal to the momenta plus this deviation, we 
find that these new objects undergo scattering events associated with the ground 
state — > excited state transition with no energy exchange. We will call these 
new objects pseudofermions. The pseudofermions have momenta q, which we 
will call canonical momenta, due to the canonical pseudoparticle-pseudofermion 
transformation which will be defined in section ()2.3.3|) . This canonical momentum 
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is defined as: 

(a) 

q = q+ ^'^^y^' OLV = cO, cu, su (2.89) 

Ij 

Note that for the ground state Q^j,{q) = since all the deviations in Eq. 
()2.(j7j) are zero. This means that for the ground state we have that q = q- The 
absence of residual energy interactions between the pseudofermions will simplify 
our calculations of the one-electron spectral function tremendously. In fact, with- 
out this property of the pseudofermions, it would be pointless to introduce these 
new quantum objects. 

Similarily to the pseudoparticle, we will define creation and annihilation op- 
erators for the pseudofermions, which lead to a formal definition of the pseud- 
ofermions number operator Afauiqj), section ()2.3.3j) . However, by physical rea- 
soning, there are some things that can be claimed without further due (see also 
Ref. [HII)- First off, since in the ground state the pseudoparticles and the pseud- 
ofermions are exactly the same objects, we have 

KM) = Nliqj) (2.90) 

and moreover, when a pseudoparticle with momenta q is found in a configura- 
tion belonging to an excited energy eigenstate, we have that the corresponding 
pseudofermion has a canonical momentum value of qj = qjiqj) according to Eq. 
()2.89|) . Since the ground state configuration of the two representations are equal, 
this means that 



AJ\fauiqj) = ANauiqj) 

^c..%)=^c..%) -KuiQj) (2.91) 

which implies that J^auiqj) = N^uiqj)- However, this does not imply that in 
the continuous system Maviq) = Nau{q), in fact this is in general not true. We 
remember that we reach the continuous system by letting {2tt/L) —>■ but that 
our pseudofermion theory carries physically relevant terms of order {1/L). 
Let us define the inverse of Eq. ()2.89|) in the discrete system, 

N", , 

27r 

Qj = (iMj) = ~ ^ $Z 5Z '^"^."''^'fe' ^■)^"V'(gj) (2.92) 

a'v' j'=l 

and then investigate the jacobian of the q ^ q coordinate transformation: 

J^FcAq) = ^ [dq F^,{q) = ^ [ dq TaM^ (2-93) 
„ 27r J 2tt J. dq 
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where Fau{q) and J^av{<i) are some functions of the momenta and canonical mo- 
menta, respectively. The jacobian becomes 



dq{,q) 
dq 




dq 




dq 



AAfa'.'iq') 



(2.94) 



which yields one only if J-'au{q) is proportional to AAfa'u'{q'), since we do not 
include second order terms in our theory. The canonical momenta spacing will 
be further discussed in section (I2.H.H|1 . 

By similar reasoning, we see that the energy bands of Eq. ()2.7H) transform 
according to 



but since the energy bands always multiply the corresponding pseudoparticle 
or pseudofermion deviation, we find that the second term of the pseudofermion 
energy band is of second order in the pseudofermion number deviations, and hence 
falls outside the realm of our pseudofermion theory. It is therefore safe to use the 
same energy bands as previously derived. 

Since we do not change the total number of pseudoparticles, we only shift the 
momenta of them, the number operators as well as the corresponding number 
deviation operators, will have the same eigenvalues in the pseudoparticle basis as 
in the pseudofermion basis '^4^. This is due to the fact that the pseudoparticles, 
whose number operators commute with the hamiltonian, have the same compo- 
sition in terms of rotated electrons as the pseudofermions. Therefore, the pseud- 
ofermion number operators should also commute with the hamiltonian. Thus, 
due to the equahty of eigenvalues and the common eigenstates, it would be ex- 
pected that by using a formal operator language, we find that the pseudoparticles 
and the pseudofermions are related to each other by a unitary transformation. 
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2.3.2 The PS subspace and quantum shake-up effects 

As already mentioned before f section |2.2.6|) . we are interested in formulating a 
normal ordered theory, relative to the ground state. This implies that all quan- 
tities needed to describe the dynamics of the model will be expressed in terms 
of deviations from the ground state configuration. Like this, we obtain different 
theories for different ground states. According to the pseudofermion picture, we 
will only retain terms of order {1/L). The fact that the pseudofermions do not 
have any physical properties of order {l/Ly for j > 2, makes it possible to for- 
mulate the normal ordered theory with only first order terms. Hence, we consider 
ground state excited state transitions such that the number of pseudofermions 
change according to 

AN^,{q) = N^M - Nliq) - qL < q < qL (2.96) 

where we use the same symbol for pseudoparticle number deviation and the pseud- 
ofermion number deviation since they are always equal to each other. The devi- 
ations of au pseudofermions, can equally be expressed in terms of deviations in 
the electronic numbers, according to: 

M,,_i=^(Ar-Ar,o) 

oo 

M,__.=V.+^i.Ar,, (2.97) 
i/=i 

which by taking deviations leads to 

oo 

AN = AN^o + 2AL^^_i + 2 ^ z/AA^,^ (2.98) 

i/=i 

Secondly, 

oo 

M,^_i=L,^_i+5^z/Ar,, (2.99) 

implies that 

oo 

A(A^I - Ni) = AN^o - 2AL^ _i - 2AN,i - 2 ^ uAN^^ (2.100) 

!/=2 

Thirdly, due to Eq. ()2.H2j] . we have that 

A(M,,+ i+M,,_i)=-AAr,o 

a(m,^^^+M,A=AN,o (2.101) 
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which, by taking deviations of Eq. ()2.34j) and equalhng the resulting expressions 
with the ones above, ehminating the deviations for the +^ ho Ions and spinons, 
respectively, leads to 

oo 

i/=i 

oo 

AA^,o - 2AA^,i - 2 ^ uNsu - L^^_i (2.102) 

i>=2 

which implies that the numbers AL^ and AL^ are not independent and thus 
it suffices to specify only L^_i and L^_i when dealing with deviations of Yang 
holons and HL spinons. We note that for the LWS ground state, Ncu = Ngu = 
L,^_i = = ^ AN,, = N,, (for z/ > 1), AA^,, = A^,, (for z/ > 2), 

AL^ _i = _i and AL„ _i = _i . 

(,,2 ^'2 '2 '2 

These relationships limit the number of pseudofermions created or annihilated 
whenever a finite number of electrons are created or annihilated. This is hence our 
first restriction: to only allow processes that create or annihilate a finite number 
of electrons. Ultimately, this will be seen in the number of electronic creation and 
annihilation operators present in the operators of any correlation function that 
we wish to calculate. Later, we will focus on the one-electron spectral problem, 
for which this issue is trivial (being only one electronic creation or annihilation 
operator). The finite electron creation and annihilation operator limitation, im- 
plies that the collection of excited states reachable by these operators, span a 
strict subspace of the entire Hilbert space of the model, a subspace that we will 
call the pseudofermion subspace, abbreviated PS. 

Having limited the number of electrons created or annihilated from the system, 
we can investigate how this change in the total number of electrons affect the 
number of pseudofermions and their lattice configurations. Within our theory, 
each au pseudofermion is only existing inside the many-body quantum system, 
and doing so on a specific au dependent lattice. This means that not only does 
the lattice constant and the total number of available canonical momentum values 
differ for each au branch, but the latter number also changes whenever creating 
or annihilating electrons. This implies that the quantum numbers describing the 
occupancies of the specific au branch under consideration, change from being 
integers or half-odd integers to being half-odd integers or integers, respectively. 
This effect is usually called the quantum shake-up effect. This effect takes place 
for all av branches, even for the cO branch in spite of the fact that N*q = Na 
is constant, because for this branch it is the number apparent in Eq. ()2.48j) . 
and not N*q, that decides whether or not the occupancies of the cO branch are 
described by integers or half-odd integers. During a transition to an excited state. 



62 



if the changes of the following numbers are odd, the quantum numbers describing 
the occupancies of the corresponding branch, change according to the quantum 
shake-up effect (derived using Eq. ()2.45|) ): 



av=cv,sv 
oo 

^{v' + v -\v' ~v\)ANsu' (2.103) 

oo 
u'=l 

The expressions of these number deviations will be simplified in section ()2.3.4j] . 
Note that the deviational numbers are purely expressed in terms of occupational 
numbers of pseudofermions. This is necessary since for the same electronic cre- 
ation or annihilation process, the resulting quantum mechanical state may be 
a linear combination of several states, with different set-ups of pseudofermions. 
Thus, in order to properly account for the shake-up effect, we need to count the 
deviations of pseudofermions in the particular state that we are investigating. 
This is a consequence of the occupational numbers of electrons not being good 
quantum numbers: the same number of electrons, as well as the same numbers of 
|-spin and j-spin electrons, can be fitted with many different quantum mechani- 
cal states, whilst a certain specified set of pseudofermionic occupational deviation 
numbers specify one and only one quantum mechanical state. 

Since these shake-up effects are measured relative to the ground state, we 
should study how they alter the momentum values in the ground state, which 
has g = g as already mentioned in section 1)2.3.111 . The "shifts" in the quantum 
numbers, which can be written as 

/-^/- + Jo Jo = -^ , 0, ^ (2.104) 

implies that the ground state momenta changes according to 

so that we can write the shift in the ground state momenta as 

+ % QL = -^,0,7r (2.106) 
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su : AiV;, = A (iV,, + iV,o) - 

cu : ANl = A (iV,, - iV,o) - 



(2.105) 



which defines the new quantity Q^u- This is not in contradiction with the defini- 
tion of the canonical momentum, Eq. ()2.89j) . What the shake-up effect entails, 
is that the momentum values, the g's in Eq. ()2.89|) . are shifted, but with the 
definition of the canonical momentum intact. One should note that if IQ^ul = ^) 
we have two resulting states that differ from each other in terms of the positions 
of the occupied quantum numbers, as shown in Fig. ()2.8|) . 

One could, and we certainly will, view the shake-up effect as producing a 
virtual excited state, which is a first step of any ground state excited state 
transition, where the quantum objects have momenta q + Q^^^/L. The true final 
state would then be the one which is described by the canonical momenta, where 
the scattering events between the pseudofermions are governed by Qauil)- 
section ()2.3.4p . where the scattering theory will be developed, we shall see that 
the usual quantum mechanical picture, "a shift in the momenta of the quantum 
particles produces a shift in the phase of the wave-function" , will also apply here. 



2.3.3 Pseudofermion operators and momentum spacing 

By using the fact that pseudofermions, with canonical momenta q, live on the 
same lattice as the pseudoparticles, with momenta q, and that furthermore there 
are as many pseudofermions as pseudoparticles for every au branch, we can with- 
out further due define pseudofermion creation and pseudofermion annihilation 
operators: 

N* 

ft _ piqja-au ft 

y ^"^av j=i 

4- = ^E^^^''"^/^-- (2.107) 

V^^au j=i 

where fj^^^ creates a au pseudofermion on a au effective lattice site position j 
(with space coordinate xj = jaau) and fj^au annihilates a au pseudofermion on 
a au effective lattice site position j. Note that there is only one q for each q. 
In other words, the function q = q{q) is unique for each ground state — > excited 
state transition, due to the uniqueness of the solution of Eq. ()2.70|) . 

The pseudofermion operators and the pseudoparticle operators are related 
through a unitary transformation Vau'- 

fq,au ^}iti/^\,au^OLV fq,au ^au^q,au'^au (2.108) 

where 

Vau = exp I ^ 6j [65,,, - } (2.109) 
q 
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Figure 2.8: Schematical figure illustrating the shake-up effect on a toy lat- 
tice (where the filled circles depict pseudofermions and the empty circles pseud- 
ofermion holes). Let US, clS Si pedagogical example, suppose that we create a 
pseudofermion au hole onto a ground state with a symmetrical distribution of 
pseudofermions around the zero momentum point (here indicated by a vertical 
dashed line), in such a way that we have a non zero shake up effect, according to 
Eq. ()2.1()H|1 . We note that, in this example, the top configuration is energetically 
more favorable than the bottom configuration. The constant spacing between the 
lattice points is 2tt/L and the size of the shake-up shift is tt/L. 
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These relationships, as well as the proof of Vau being unitary, are given in Ref. 
|94j . The pseudofermion picture thus allows us to write the energy as 

E = Egs + AE = Egs{{nI}) + AE{{AN^,}) (2.110) 

without any higher order terms in the deviations and 

= (j2^lauKau) = {Y.fio.uhau) (2.111) 



Before we move on, there are some important properties that need to be 
clarified. The name pseudofermion stems from the fact that the operators f^,^^ 
and fq^au^ just like the pseudoparticles, satisfy the fermionic anticommutation 
relations almost. To evaluate the anticommutator between these two operators, 
we use the fact that their local counterparts satisfy the fermionic anticommutation 
relations exactly. 

{fL.. = %^ E E e--(--"'^'H4.' (2.112) 

By using the following equalities 

{fj,au^ fj',au} = Sjj' {fj,au, fj',a'u'} = 

E^" = ^^^73r ^2.113) 

we arrive to the following anticommutators 

sin ^ 

/ft f, -1 ia^Aq-q')/2 i{Q^Aci)-Qo.Aq'))/2 ^ ^ 



au sm ( ^ ' 

{fq,auJq',a'u'} = (2.114) 

which will play a key role in the development of the theory. That the local pseud- 
ofermions satisfy the fermionic anticommutation relations exactly is a property 
due to the rotated electrons, whose number operators commute with the hamil- 
tonian pUnj . 

The second property of the pseudofermion description that we will mention 
here regards the discrete canonical momenta spacing Aq. The discrete momenta 
spacing is, just like for free fermions, constant and equal to (2,71 /L): 

Aqj = Qj+i -qj = -j- (2.115) 
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The discrete canonical momenta, however, satisfies (for au = cO, cu, sv) 

Aqj = - = Ag, + ^"^^^^'^ (2.116) 
where the difference (Q*(gj+i) — Q'^{qj)), due to Eq. ()2.67|1 . is 

7 , / dq' [^au,a'u'iqj+l, q ) - ^aiy,a'u'{qj, ? )] ANa'u'iq') = 

AN^'u'iq') (2.117) 



„0 



which means that, since Aqj = 27t/L, 

Aqj = ^ + 0(1/ L') (2.118) 

or, in other words, since {l/L"^) terms do not contribute to the physics in the 
pseudofermion picture, we are erroneously led to conclude that Ag^ = Aqj, or 
that the discrete momenta and the discrete canonical momenta are equal to first 
order in {1/L). These are erroneous conclusions based on the fact that for larger 
deviations, we have qj/ — qj = {2tx / L){j' — j) for the momenta, whilst in general 
qji — qj 7^ (2,71 /L){j' — j) for the canonical momenta. This is most easily seen 
when j' and j are very far apart, say j' — j ~ {Na/2). This means that we have 



qr -qj- Y 



^ + 0(l/L2) 



~- + 0(l/L) (2.119) 



by averaging the values of the derivatives of $Q!/,a v {q, q') which yields a factor 
proportional to j' — j. This difference between Aqj and Aqj is of the order of 
(1/L), which is a difference inside the realm of the pseudofermion physics, i.e. 
a non-negligible difference. This difference implies that whenever we want to 
replace a sum by an integral, i.e. when A^^ ^ oo we are not allowed to use the 
"standard" replacement Aq (27r/L) but instead we will need to use a state 
dependent jacobian (due to the state dependence of [(5*(gj+i) — Q'^{qj)\). 

Normal ordered operators are sometimes written as :X: with the definition 
:X: = X-X^, for X^ = \ {GS\X\GS)\ (where \GS) is the ground state). A typ- 
ical eigenstate of such a normal ordered operator is thus the deviation between 
the eigenvalues of the excited state and the ground state respectively. Since our 
pseudofermionic operators are normal ordered by construction, it seems super- 
fiuous to use this notation for them, whilst the hamiltonian and the momentum 
operator becomes :H: and :P: respectively. 
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In the pseudofermionic basis, these operators are 




j=l u=l j=l 



(2.120) 



2.3.4 Virtual states and pseudofermionic subspaces 

To describe the scatterers and the scattering centers of the theory, we need to 
describe the excited eigenstates in terms of pseudofermions. We will follow the 
standard non-relativistic description of a quantum scattering theory ^7\, in which 
the scattering S-matrix, which maps the "incoming" quantum state into the scat- 
tered "outgoing" quantum state, will play a central role. Since all our pseud- 
ofermions are either ?7-spin zero or spin zero objects, the scattering matrix will be 
of dimension 1x1, i.e. just a complex number, in contrast to the representation of 
Ref. |SS1, where the scatterers are ?7-spin | and spin | objects. In this reference, 
the S-matrix has a larger dimension due to the the coupling of the 77-spin and the 
spin channels. In the following, we will use the usual definition of a phase shift 
such that a shift in the momentum Si{q) of a quantum object with momentum q, 
produces a shift in its wave-function equal to e'^^'-^'^\ where / stands for a collection 
of quantum numbers used to fully describe the original (unscattered) incoming 
wave (also known as the in asymptote) . The job of the S-matrix is then to supply 
the incoming wave with this phase shift, and transform it into the outgoing wave 
(also known as the out asymptote) [HZj. Since the incoming and outgoing waves, 
according to the general quantum scattering theory, preserve the total momenta 
and the total energy, we see that we must introduce another scheme when de- 
scribing the transitions. This is because the quantities that we have associated 
with a general ground state —>■ excited state transition, namely the energy and 
momentum deviations according to 



do not preserve neither the total energy nor the total momentum. We therefore 
divide the entire transition into two steps: one scatteringless step and one in which 
all the scattering events occur. The scatteringless step yields a virtual state, or 



E = + E(^) = Egs + 
p = p(o) + p(i) = + AP 
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an intermediate state, which brings the system from having energy Eos and mo- 
mentum Pgsi to having energy Egs + ^E and momentum Pgs + ^P respectively 
|98j . Under this ground state — > virtual state transition, the quantum numbers 
describing the occupancies of the pseudofermions may "shake-up" as described 
in section ()2.3.2j) . Also, the energies and the momenta change according to Eq. 
()2.121|) . i.e. the ground state virtual state transition is a scatteringless finite 
energy and finite momentum transition. The virtual state is the in asymptote in 
the scattering theory. Thus, this virtual state is occupied by pseudofermions with 
momenta q + Q^^^/ L. Then, the virtual state will undergo scattering events gov- 
erned by the quantities Q%y{q), for av = cO, cu, su, that preserve total energy and 
total momentum. In other words, the additional state dependent shift Q%j^{q)/ L 
implies no extra energy nor momentum terms in the deviation expansions. 

Before we move on, we need to specify the pseudofermion deviations that 
characterize a typical virtual state. Since electrons are the only quantum objects 
that can be created or annihilated, we have to classify the types of subspaces we 
obtain by fixing the deviations AA^ and A(A^| — Ni). This has actually already 
been done in section ()2.3.2|) . where we saw that the electronic deviations can 
uniquely be expressed as 

oo 

AN = AN CO + 2AL^_i + 2 ^ uANcu 

oo 

A[N^ - N^)=ANco - 2AL,__i - 2AA',i - 2^ z/AA^,^ (2.122) 



u=2 



which means that for each fixed set of numbers {ANau} au=co,cu,su and {AL^ _i }a=c, 
we span one strict subspace of the entire Hilbert space of the model, which cor- 
respond to the actual electronic deviations at hand. One should then collect all 
possible sets of these pseudofermionic deviation numbers, to arrive to the com- 
plete set of virtual states that emerge due to nonzero deviations in Eq. ()2.122j) . 

The total energy and total momentum acquired during a ground state 
virtual state transition, is easily obtained from Eq. ()2.58|) : 

oo „gO^ 

AE = —\2 V l^dqeau{q)AN^u{q)+uJo 

2'''" X J-q° 

"—C,S U=dct,s ""^ 

' dq qAN,o{q) + dq qAN,,{q)+ (2.123) 



27r 



5:/,?'(^')^'^<■'W 



a 
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where we have defined the minimum energy excitation ooq as 

uo = _i + 2noh{M^ _i - Nsi) (2.124) 

where = M^ _i — Ngi = for the initial LWS ground state. As expected, 

this term is nonzero if we have broken the 5*0 (4) symmetry of the model. Hence, 
ujq serves as a gap parameter that tells us whether or not our excitations live in a 
gapped or in a gapless system. 

Note that for these expressions, we have qj+i—qj = 271 /L due to the scattering- 
less property of the ground state — > virtual state transition. However, the actual 
occupancy positions may shift globally according to the shake-up effect. The con- 
ditions on whether or not the quantum numbers for a particular au branch are 
shaken up can be simplified as compared to the expression given in Eq. ()2.103p . 
Since (z/ + z/' — |z/' — z/|) is always an even number, we can exclude the u summa- 
tion in Eq. ()2.1()H|1 altogether since we are only interested in terms that have a 
possibility to be odd. This is also the reason for why we neglect — 2AiVco in the 
cu case, after having added and subtracted AA^^co on the right hand side of that 
equation. We thus arrive to the following statement: if the following deviations 
are odd, in connection with our ground state — > virtual state transition, then 
the quantum numbers for the actual virtual state change from being integers (or 
half-odd integers) to being half-odd integers (or integers): 

cO: 

ai'=cu,su 

au : ANo^u + AA^co (2.125) 

2.3.5 The S-matrix 

According to the standard quantum mechanical scattering theory, the S-matrix is 
a unitary operator that maps the pre-scattered state |0i„) into the post-scattered 
state |0oMt)- As we have seen above, these states have many names, according 
to the rich history of scattering theory in general. An example of frequently 
used names are "incoming" waves and "outgoing" waves, due to the classical 
analogue of colliding billiard balls (thus incoming balls and outgoing balls). A 
more mathematical nomenclature includes "in asymptote" and "out asymptote", 
due to the mathematical formulation of modern scattering theory. In this case 
the pre-scattered state is regarded as the "untouched" state that existed at a 
time t — oo (i.e. as the time t approaches this limit, the state approaches some 
asymptotic idealized form) and the post-scattered state is similarily regarded as 
the asymptotic state at time t +oo. From now on, we will choose this latter 
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nomenclature for the pre- and post-scattered states. Thus in general, 



5 1 0m) = \(t)out) 



Sw = e 



2i5, 



hi' 



(2.126) 



where Sui is the total phase shift for the / — > V scattering process and 5*//' is 
the corresponding matrix element of the S-matrix (/ and are sets of quantum 
numbers that fully describe the scattering quantum objects). In our case, since 
the different av branches do not mix with each other, we can define a Say{q) 
for each branch at momentum q. This quantity describes the scattering events 
of a av pseudofermion of momentum g, with pseudofermions of all other a'v' 
branches created by the transition, as it travels around the lattice. We shall see 
that this operator can be described by products of quantities called Sau,a'u' (<?, q') , 
that gives the form for the individual scattering events between the au and the 
aV pseudofermion. This should not be confused with the usual scattering notion 
of "mixing between different scattering channels". In our theory, there is no such 
mixing: we describe scattering events between r^-spin and spin zero objects, which 
preserve the individual au branches. Moreover, since the scattering itself does 
not change the energies nor the momenta of the scatterers nor of the scattering 
centers (upon which the scatterers scatter), the scattering events are of a trivial 
zero energy forward scattering type. We will clarify this claim in more detail later 
in this section (however, we refer to Ref. |l()4j for the main results and to Ref. 
jHH for a detailed analysis). 

Formally, in finding our expression for the total phase shift 6au{q), we note 
that the phase of the in asymptote changes as our pseudofermion scatters with 
all the scattering centers of the system, i.e. as the scatterer travels around the 
lattice ring once, to arrive to its original starting position (remember that we 
adopted periodic boundary conditions for the original hamiltonian ()2.Hl ). There 
are different choices of coordinates available for this picture, giving different defin- 
ing expressions for 6au{q), however always yielding the same S-matrix jHE]- We 
should note that all of the au branches live on au effective lattices with the same 
lattice length L. We will choose the following coordinates: let our pseudofermion 
depart from lattice position x = —L/2 and arrive at lattice position x = L/2, 
with a phase difference equal to Sau{q)- We have 




(2.127) 



which yields 



qL 
2 



qL 
2 



q + 



QL + Qtiq) \ L qL Ql^ + Qtjq) 
L J 2 2 2 
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by comparing the momenta shifts between the ground state and the final state. 
Hence we can write the total phase shift of an au pseudofermion as, Eq. ()2.(j7j) . 




= % + ^ I] X5 "^-.-'-'fe' qj')ANa'u'iqj') (2.129) 



and we obtain, by the formal definition of the S-matrix ^7\, that 




av j'=l 



From this equation, we see that 7f^au,a'u'{qj,qj') measures the phase shift of the 
au pseudofermion at momentum qj due to the individual scattering event with 
the aV pseudofermion at momentum qj/. Note that it is the latter (primed) 
pseudofermions, that were not present in the original ground state, that make 
the S-matrix to differ from unity and hence it is these pseudofermions that are 
the scattering centers of the theory. The au (unprimed) pseudofermions, on 
the other hand, are the scatterers of the theory. Note here that would we have 
chosen our pseudofermion to originate from x = 0, and travel around the lattice 
ring until x = L, the resulting total phase shift would have been the same as 
above multiplied by 2. However, the S-matrix would have remained the same by 
letting the phase shift of the individual scattering event be twice the expression 
given above, 2TT^au,a'u'{qj,qj')- Moreover, there is nothing in this picture that 
distinguishes pseudofermions from pseudofermion holes, which means that an au 
pseudofermion hole is also a scatterer on equal footing with the pseudofermion. 
In other words, whenever ANa'u'iqj') < 0, we have that it is pseudofermion 
holes that act as scattering centers. Thus, the number of pseudofermions and 
pseudofermion holes for which the S-matrix has the form of Eq. fl2.1H()j) . equals 
the number of lattice sites for every non-empty branch (i.e. for every branch not 
consisting entirely out of holes). 

Thus, the pseudofermion or pseudofermion hole S-matrix for the one dimen- 
sional Hubbard model is just a phase factor, given by Eq. ()2.130|) . This statement 
is consistent with the previous claims that the scattering events do not mix dif- 
ferent ?7-spin or spin channels, and that the scatterers as well as the scattering 
centers are ?7-spin and spin zero objects. Moreover, it is not only the total energy 
and the total momentum that is conserved during these scattering events, but 
also the individual au energies and canonical momenta components. This can 
most easily be seen on the form of the energy and momentum deviations, Eq. 
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()2.123|1 . that do not mix pseudofermion deviations from different av branches 
after the substitution q ^ q. 

Finally, we should note that the pseudofermion anticommutation relations of 
Eq. (j2.114|l . can be solely expressed in terms of the pseudofermion S-matrix, 
according to 



Sau{q) 



nl/2 



Im 



Sau{q) 



1/2 



sin[aai.(g - q')/2] 



(2.131^ 



with the important implication that the S-matrix introduced here fully controls 
the one electron spectral properties of the normal ordered ID Hubbard hamilto- 
nian, as will later become apparent by use of the anticommutation relations in 
the evaluation of matrix overlaps in chapter Q . 



2.3.6 Properties of the pseudofermion scattering 

Before closing the section on pseudofermion scattering theory, there are some 
properties that the theory implies which is worth mentioning, here numbered 
from (i) to (v). 

(i) As noted by the explicit form of the S-matrix, we have reduced the many- 
body scattering events into two-body scattering events, as shown by the definition 
of SaiJ,a'u'{qj, qj') in Eq. ()2.13()j) . This means that the relative ordering between 
any pair of two-body scattering events is independent of the final expression for 
the S-matrix (mathematically due to the commutativity of complex numbers) 
|98j . This is a stronger result that what an S-matrix satisfying the Yang-Baxter 
Equation |lUUj - piI^ could claim [22] |lU3j . In these references, a representation 
different from the one of our scattering theory regarding the active scattering cen- 
ters is made. Indeed, the scatterers and scattering centers of that representation 
are "spinous" and "holons" with ?7-spin and spin projection equal to whilst 
in our theory the scattering centers are //-spinless and spinless. This explains 
the differing dimension of the S-matrix in those references, as compared to our 
representation. 

(ii) One should note that the ±| Yang holons and the ±| HL spinous have 
not played any role in the pseudofermion scattering theory. On the contrary, they 
have constant momentum values during the ground state —>■ virtual state — > final 
state transitions, only one of them nonzero (the — | Yang holon with momentum 
equal to tt). This means that the S-matrix for these quantum objects equals unity, 
due to the absence of phase shifts. We remind ourselves that these objects are 
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exactly the same in the original electronic frame as in the rotated electronic frame, 
and thus have no quantity corresponding to the AQ{q), Eq. ()2.(i7|l . previously 
defined for the quantum objects described by the Takahashi string hypothesis, 
Eq. ()2.57|1 . We thus conclude that the ±| Yang holons and the ±| HL spinous 
are neither scatterers nor scattering centers. 

(iii) There is an elegant theorem, called the Levinson's theorem [105J, which 
states that in the infinite wavelength limit, i.e. when the momenta of the scatterer 
tends to zero in the reference frame of the scattering center, the phase shift 
becomes an integer multiple of vr, where this integer is nothing but the number of 
bound states Nf,. In our notation, this means that the momenta of the scattering 
center q' should be replaced by and that the momenta of the scatterer q should 
be replaced hj q — q' and we should thus have 

lim 9MS^LS1 = n, (2.132) 

according to the theorem, since we in our case have that the scatterer feels the 
effect of the scattering centers during the virtual state final state transition 
only. For our theory to comply with this theorem, we should have that the 
above limit is equal to zero, since we by construction have no bound states in 
the theory (this can also be seen mathematically: our S-matrix has no poles). In 
section ()2.2.7|1 it was found that ^au,a'u'i(l,<l') = —^au,a'u'{—(l, — <?') which in the 
alternative reference frame translates into 

^a,y,a',.'iq- q',0) = -^au,a'u'i-iq- q'),0) ^ hm " ^'^ 0) = 

(2.133) 

which means that Eq. ()2.132|) is fulfilled for our scattering theory. 

(iv) There is another exact result for which we can check our derived results, 
known as the Fumi theorem |106j . |lU7j . This theorem, originally formulated for 
electrons in a metal, states that the total energy Ei due to the existence of an 
impurity upon which otherwise free electrons scatter, can be written as an integral 
over all the phase shifts caused by this impurity: 







where E{k) ~ k"^ for free electrons, and / represents a set of relevant quantum 
numbers (for three dimensional scattering events of electrons, / usually denotes 
the angular momentum components). 

In addition to the pseudoparticle and the pseudofermion representation, one 
can introduce a third related description in terms of quantum objects that carry 
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rapidity momentum kj. It will be shown elsewhere jOH] that the phase shift of 
such quantum objects obey the Fumi theorem. The energy Ei above will then 
correspond to that part of the energy deviation E^^\ which is of scattering origin. 

The occupied rapidity momentum kj = k^{qj) obeys for the ground state 
—Q < kj < Q, i.e. the Fermi points of these quantum objects are defined by 
the value Q. Moreover, their dispersion relation goes as E{k) ~ cosfc, and their 
phase shift is Sco{k) = Qco{k)/2, where Qco{k) is the quantity equivalent to Qcoil) 
in this representation. 

(v) The following final properties will only be briefly mentioned here, since 
they correspond to quite exotic cases of the pseudofermion theory, and are only 
valid for the cu {v > 1) and the su {u > 2) branches and thus will not be 
considered in the dynamical theory. The interested reader should go to Ref. 
^98ij for a complete analysis. In section ()2.2.7|1 we saw that the energy bands 
in the case of au = cu, sv equal zero for momenta equal to the limiting value 
of the av effective Brillouin zone. Thus at these momenta points, the energies 
of the cv and the sv pseudofermions becomes the sum of the energies of the 
individual quantum objects of which they are constituted. One can understand 
this by the "handwaving" analogue that the "binding energies" between the ±| 
ho Ions (a = c) and ±^ spinous [a = s) vanish, so that there is nothing to hold 
the pseudofermions together. Interestingly enough, one can show that for these 
momentum values, and for < na < 1 and < ma < n, we have the following 
equality 

±Aq = ±Aq + ^""^ ^""^ = (2.135) 
Ij 

which implies that the canonical momentum spacing vanishes, as opposed to 
canonical momenta spacings at other points in the au canonical momentum 
Fourier space. We should recall that the corresponding pseudoparticle shift is 
nonzero whenever AN*^ ^ 0, i.e. whenever we have a finite shake-up. This 
means that not only does the cu and su pseudofermions break up at the bound- 
aries of the effective Brillouin zone, but their canonical momentum values at these 
boundaries are the same for the ground state as for the final state, as there are no 
nonzero momenta spacings that allows a shift in the canonical momentum values, 
thus becoming non dynamical. 

This phenomena is ultimately demonstrated by the fact that at these canonical 
momentum values, the cu and the su pseudofermions become invariant to the 
unitary operator transforming electrons into rotated electrons. There is however, 
in the many-body system, some "memory" of these pseudofermions left, since it 
is possible to express ^ai',a'u'{qFau, i^Qa'u') (where a'u' = cu' or su' and l = ±) 
solely in terms of ^au,co{lFau,±qFco) and ^au,si{lFau,±qFsi), [98 . This means 
that even though the cu and su pseudofermions fall apart into their constituents. 
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their momentum values are "carried over" in the system, in such a way that 
the cO and the si pseudofermions feel the usual two-body scattering events with 
these pseudofermions as if they were cO and si pseudofermion scattering centers, 
respectively, at their corresponding Fermi points. 

By letting na ^ 1 and/or ma 0, we have that the cu and/or the su 
bands shrink until they finally disappear in the half filled case (ci/) or in the case 
with zero magnetization (su). This is easily seen by the fact that the limiting 
momentum values for the two bands become zero, i.e. for cu we have that (7r/a) — 
2A;i? — > whilst for su we have that kpi — kpi ^ 0- 
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Chapter 3 

Pseudofermion Dynamical 
Theory 

3.1 Excited energy eigenstates 
3.1.1 Introduction (spectral function) 

The goal of this chapter is to derive a pseudofermion dynamical theory which 
will enable us to calculate the one electron spectral function for the Hubbard 
hamiltonian ()2.H1 . By "one electron spectral function" we mean the spectral 
function for one electron removal and one electron addition, respectively. The 
removal function gives rise to spectral weight in the so called Removal Hubbard 
Band (RHB) whilst the addition function gives rise to the Lower Hubbard Band 
(LHB) and the Upper Hubbard Band (UHB), respectively. For (U/t) — >• oo, these 
latter two bands are separated by an energy proportional to the effective Coloumb 
interaction strength U, since the LHB gives the spectral weight for an electron 
added at an empty site whilst the UHB gives the spectral weight for an electron 
added at a singly occupied site. There is only one band in the removal case due 
to the fact that our ground state, upon which we act with suitable electronic 
creation or annihilation operators, is a LWS that is void of doubly occupied sites. 
In this thesis report, the spectral function for the RHB and the LHB will be 
calculated, even though the method used here is general and can perfectly well 
be applied to the UHB as well as to correlation functions involving creation or 
annihilation of several electrons [021 • We will use arbitrary values for the 
parameters (U/t) and n, but we will keep a small magnetization ma > and 
later let ma — ^ (to confine ma to zero at the start of the calculations has been 
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seen to be quite pathological, for reasons that we will give later). Formally at 
zero temperature, a spectral function is defined as the imaginary part of the time 
ordered Green's function G{k, uj) at electron momentum k and electron energy u, 
multiplied by a constant for which there is no conventional fixed value, but that 
is uniquely defined by applying suitable sumrules. These sumrules stem from the 
fact that the spectral function is interpreted as a probability function, and that 
thus the integral over the domain of this function must equal a certain positive 
value. The definition of the one electron removal spectral function B^{k,uj) and 
the one electron addition spectral function B^{k,uj) is 

B-ik,co)= 5^ 5^|(/_|cfc.|G5)|'5(cu- AE_) (RHB) 

i?+(fc,.;)=^ 5^|(/+|cL|G5)|'5(^-Ai?+) (LHB) (3.1) 

<^=T,i /+ 

where the summation over a will yield nothing but a factor of 2 in the zero mag- 
netization limit since creating or annihilating a |-spin electron will give exactly 
the same spectral function as creating or annihilating a |-spin electron. These 
spectral functions are then directly proportional to the probability of finding the 
added electron or the added electron hole at momentum k and energy uj, respec- 
tively. Many times we will summarize these two functions by using I = ± (which 
we will treat equivalently to / = ±1): 

B^{k,uj) = J2\{fi\cL\GS)\'6{u; - AE,) I = ± (3.2) 

fl 

where c^a = and = c^^. \fi) denotes a final state, i.e. the energy eigen- 
state of the + / electron system, where N is the number of electrons in the 
ground state, here denoted by \GS). Due to the nonzero phase shift of the pseud- 
ofermions, the states |//) and \GS) have different boundary conditions for each 
final state and au pseudofermion branch. This implies that the evaluation of the 
matrix overlaps leads to the orthogonal catastrophe, originally due to the canoni- 
cal momentum shifts of order (l/L) of the theory |132j |133j . Since the scattering 
phase shift is state dependent, we would expect a different contribution due to the 
orthogonal catastrophe for each ground state final state transition. AEi is the 
energy difference between the ground state and the final state, defined according 
to 

AEi = I {Ej^ - Egs) (3.3) 

This definition of AEi measures the energies relative to the chemical potential, 
that hence never enters the calculations explicitly. The relation to the zero tem- 
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perature Green's functions is 



B-{k,u) 



1 



lm{G{k,uj)} 



TT 



1 



lm{G{k,uj)} 



(3.4) 



to>0 
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where Im is the imaginary part. The Kramer- Kronig relations give us the inver- 
sion of these relationships, expressing G{k,uj) in terms of B\k,uj): 



where 6o is a positive infinitesimal quantity. This Green's function can, at least 
formally, be used to obtain the expectation value of any one electron correlation 
function The momentum distribution function is just the u integral 

over B^{k, uj). Integrating this function over k gives then the density of electrons 
na. The remaining spectral weight, from the LHB and the UHB spectral functions 
respectively, must then by construction have a k and u integrated value of 2 — na. 
These integral values constitute the sum rules for the spectral functions 



where the approximative sign in the last equality stands for the (very) weak {U /t) 
dependence in the sum rule This dependence is due to the allocated weight 

in the UHB, which varies slightly as {U /t) varies, and is approximatively equal 
to na. We see that at half filling, na 1, the LHB weight vanishes, as all the 
one electron addition spectral weight is transferred to the UHB. 

One could argue that since the exact wave function of the model is known 
■ the calculation of i?'(fc, a;) is just a matter of explicit brute force calculation. 
Unfortunately, however, the complex form of this wave function and the fact that 
it has remained unknown how to express the generators of the excited energy 
eigenstates in terms of electronic creation and annihilation operators, it has so 
far been practically impossible to calculate B''{k,uj) by brute force. This means 
that, from a theoretical standpoint, we have to choose an approach between 
either finding alternative methods, using the exact solution in some limit that 
simplifies the expressions, or discovering the missing link between the electrons 
and the quantum objects that diagonalize the normal ordered hamiltonian. The 
pseudofermion theory allows to do the latter, but the main contributions in the 
literature on the subject has primarily been focused on the former. 




(3.5) 




(3.6) 
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In Ref. |112j a lattice dividing technique is used to calculate the spectral 
weight for small system sizes, with (U/t) = 4, raa = 1 and na = |. The pur- 
pose of dividing the lattice into small "clusters" is that the Green's function can 
be obtained by exact diagonalization, when the system size is very small. The 
" intercluster" hopping integral is then treated perturbatively to obtain the full 
Green's function. In the perturbation theory, the exactly solvable hamiltonian is 
taken to be the entire Hubbard hamiltonian for one cluster (using cluster size of 
12 lattice sites). Like this, a spectral weight for the cases of RHB and LHB with 
somewhat distinguishable spin and charge dispersions is obtained, even though 
the shape of the two dimensional surface in the {k, u) plane could be improved. 

Another method was used in Ref. using parameter values suitable for 

comparison with experimental results on the charge transfer salt TTF-TCNQ, 
namely filling na = 0.59 ~ 0.6 and effective Coloumb repulsion {U/t) = 4.9. For 
this the "dynamical density matrix renormalization group method" (DDMRG) 
was employed. With open boundary conditions and system sizes up to 90 lattice 
sites, the spectral weight was calculated by using the eigenstates of the particle- 
in-a-box problem in the DDMRG routine. A more reliable association of the 
spectral weight with different quantum object dispersions is made, as well as 
some estimates for the exponents with which the spectral function diverges along 
the dispersive lines. Unfortunately, the DDMRG routine becomes non applicable 
as the system size approaches the thermodynamic limit. 

In Refs. |114j and |115j . the one electron spectral function is investigated 
using Green's functions, conformal field theory (briefly described below) and the 
Bethe ansatz solution in the half filled Mott-Hubbard insulating phase, for finite 
values of {U/t). Using a holon and spinon picture, some lines along which these 
holons and spinous disperse are identified to display singular features of the spec- 
tral function, with momentum line shape dependent exponents (however, these 
holons and spinous are different than the quantum objects dubbed "holons" and 
"spinous" in this thesis report). 

In Refs. jZS|-[Hn!, the RHB and the LHB spectral functions were calculated 
in the {U/t) oo limit, by using the exact solution and exploiting the wave 
function factorization in that limit. The quantum objects that describe the oc- 
cupancies of the two resulting parts of the factorized wave function are called 
"spinless fermions" and "spinous", respectively, and account for the spin-charge 
separation in the {U/t) — > oo limit. The momentum of the spin wave, obtained 
by mapping the spin part of the Hubbard hamiltonian to the ID Heisenberg 
spin hamiltonian, imposes a twisted boundary condition on the otherwise peri- 
odic lattice of the spinless fermions. This is however, the only remnant of the 
coupled Takahashi equations ()2.7p . Moreover, the spin spectrum collapses as 
demonstrated by the fact that the group velocity of the propagation of the spin 
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wave goes to zero in this limit. The "spinless fermions" introduced in these ref- 
erences are nothing but the cO pseudofermions as {U/t) —>■ oo. Thus they are 
not "real" fermions, but rather quantum objects obeying the Haldane statistics 
|in8j . introduced in section ()2.3.3|1 . However, these quantum objects account for 
all the excited states available since at infinite repulsion there can exist no cu 
pseudofermions at finite energy, and at zero magnetization there can exist no su 
pseudofermions, due to the fact that the su band is non existent for the ground 
state at zero magnetization. The technique of calculating the matrix overlap of 
the charge part, expressing the full matrix overlap in the spectral function as a 
determinant of the anticommutators of the spinless fermion operators, will in this 
work be generalized to finite values of (U/t). 

Some of the main difficulties of this model lay in the fact that the Hubbard 
hamiltonian cannot be treated by perturbative methods due to the non pertuba- 
tive character of the electronic interactions. Indeed, in ID the Coloumb interac- 
tion parameter U, however weak, qualitatively changes the correlations between 
the electrons, as compared to the free system. By restricting the Hilbert space 
to low energy eigenstates only, we can apply various methods that ultimately 
depend on the linearization of the elementary excitation energy bands. The as- 
sumption one makes is that all relevant low lying excitations can be constructed 
by taking into account states with momenta close to i/cj? only. The technique 
of bosonization separates the hamiltonian into two bosonic hamiltonians, one de- 
scribing the charge part and another the spin part. In this way, the problem 
is reformulated into two massless bosonic theories describing the charge (c) and 
spin (s) degrees of freedom, respectively, with dispersions u!c,s{k) = Vc,s (k ± /ec.s), 
where f c,s and fcc,s are the Fermi velocity and Fermi momentum of the c (charge) 
or s (spin) branches, respectively. The charge-spin separation of the electronic 
degrees of freedom is thus explicit for the bosonization technique |116j - PTH] . In 
the case of weak coupling, it is then possible to compute the critical exponents of 
the correlation functions for the Hubbard model |119j - P^^ . 

Other powerful techniques used to calculate correlation functions in the same 
low energy excitation regime are conformal field theories and finite size scaling. 
The basic ideas of these techniques are simple scaling arguments, due to the fact 
that at large distances the behavior of correlation functions does not depend on 
the microscopic hamiltonian. Moreover, correlation functions for systems with- 
out any internal scale have to decay algebraically, for example as simple power 
laws, due to the universality class of the Hubbard hamiltonian. The exponents 
of these power laws of the conformal theory are then used to obtain finite size 
corrections of the energy and momentum |123j - [T^ . In this way, the low lying 
excitations can be obtained as "towers" of states by adding {±2ttVc^s/ L) to the 
energy and {±2tt/L) to the momentum |134j . The spectral function and its ex- 
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ponents have been described near the Fermi points of the elementary excitations 
(usually referred to as charge and spin excitations respectively) |127j - P^n] . In 
these references, the spectral function of the related Tomonaga-Luttinger model 
was examined, and yielded the characteristic Luttinger-type power-law behavior 
in the vicinity of the elementary excitation Fermi energies: [uc^s — Vc,s {k ± kc^s)]^^ , 
where the exponent 71 is given in these references. 

In the following, we will attempt to calculate this spectral function for all en- 
ergy scales, i.e. for the entire {k, u) domain. Obviously, we will exploit some of the 
ideas briefly discussed here, for example, we will see that for each pseudofermion 
branch, the spectral weight close to the dispersive lines in the {k, uj) plane obeys 
simple power law behaviors whose exponents are related to the tower of states 
close to these lines. In section fl3.1.3|) we will classify the processes leading to 
the final states of the model, as well as the various subspaces that these final 
states span. However, we will simplify many expressions of the general theory (as 
presented in Refs. |13Uj and |131p . These simplifications are heavily dependent 
on the findings of Ref. fSH], in which the partial sum rules, i.e. contributions to 
the total sum rule from different excited state subspaces, are measured. Thus the 
use of the theory developed here will involve approximations in terms of compli- 
ance with the sum rules, however all relevant features of the one electron spectral 
function will be accounted for. The theory presented in the following is formally 
developed in Refs. |13Uj and |131j . 

3.1.2 Fourier transform and rotated electrons 

Since the final states are energy eigenstates of the hamiltonian, which in turn 
is diagonalized in the pseudofermionic basis, it would be suitable to describe all 
quantities in terms of pseudofermions. Thus, apart from the description of the 
final states in terms of occupancies of pseudofermions, we need to express the 
electronic creation and annihilation operators in terms of their pseudofermionic 
counterparts. The goal is to allow for a unique description of the generators 
of all relevant eigenstates in terms of pseudofermionic creation operators acting 
onto the vacuum. Our first step however, is a little bit more modest. Since 
the unitary transformation that maps electrons onto rotated electrons is defined 
with local operators, we need to Fourier transform the c|j,^ operators appearing 
in the defining expressions for the spectral functions, into operators creating or 
annihilating local electrons. Using Eq. ()2.2|) . we obtain a sum over the lattice 
sites j inside the spectral function which, due to the translational invariance of the 
system, reduces to Na times one typical term of the sum, say the term with j = 0. 
Also, since the spectral function continues to be a function of the momentum k, 
the above mentioned translational invariance introduces a Kronecker (5-function: 
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6, 



I OkAP, 



(3.7) 



B\k,u;) = NaJ2\ifi\cL\GS)\Huj-AEi 
fl 

The next step is to express the electronic operator in terms of rotated electronic 
operators. Even though a closed form expression relating the former in terms 
of the latter is unknown, there are some things that can be done to shed some 
light on the procedure. By using the {U/t) 3> 1 expansion of V{U/t) presented 
in section ()2.2.2|) . we have for example to first order, that 



_ _y(i) 



+ 



+ ... 



+ 



1 

U 



T-u, Ci, 



+ ... (3.8) 



which after introducing the explicit expressions for Tu and T^u and evaluating 
the commutators yields 
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where cr = —a. Note that as expected. 



(3.9) 



(3.10) 



which expresses the fact that the electron - rotated electron unitary transforma- 
tion becomes the identity transformation in this limit. This shows an example 
of how to replace the electronic operator with the rotated electronic operator. 
However, we will not be depending on the large- (f//t) expansion of the electron 
- rotated electron unitary transformation from now on, but instead use physical 
reasoning when introducing the rotated electrons into the problem. Our basic 
consideration is based upon the results of Refs. jSH], |13nj and |131j . In these 
references it is shown that for the one electron spectral weight, the substitution 
c-^ accounts for over 99% of the total spectral weight, as measured by 

the sum rules (I3.6|l . This does not mean however, that we let (f//t) = oo by 
the erroneous assumption that we only keep the first term of the expansion of 
V{U /t), Eq. ()3.9|) . On the contrary, all quantities will be evaluated for the actual 
value of {U/t) that the original problem refers to. For example, the phase shifts 
Qauil) are strongly dependent on (U/t) via the functions ^au^a'u'ig, q'), and there 
will be no limiting procedure in evaluating these functions, which means that all 
exponents will inherit this [U/t) dependence as well. 

The substitution c'^ is just a statement of the fact that the subsequent 

terms only contribute marginally to the total spectral weight, with the i = 1 term 
contributing not more than 1% to the sum rules. However, one could argue that 
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since we are measuring deviations from the exact sum rule there could still be 
some singular behavior that is left unaccounted for. This could be the case with a 
very strong (i.e. narrow) singularity with a dominant contribution as compared to 
other singular features but with a small contribution to the total sum rule. That 
this is not the case is easily concluded from the type of terms in the expansion 
of Eq. ()3.9|) : the higher order terms are generated by particle-hole processes of 
rotated electrons and do not bring about any new types of excitations that could 
lead to some sort of critical behavior that the first term does not bring about 
(an example of a new type of excitation would be, for example, a net creation or 
annihilation of 2 rotated electrons in the {t/U) term). Eq. ()3.7|1 is then rewritten 
as 

B\k,u) = NaY,\{h\~^o.\GS)\^H^ - ^Ei) 5uAP, (3.11) 
h 

where the final states are described by occupations of rotated electrons, which we 
wish to reformulate in terms of pseudofermions. 

3.1.3 Classification of the eigenstates of the model 

The purpose of this section is to "sketch" the decomposition of the state sum- 
mation appearing in the expression for the spectral function, into summations 
over subspaces defined by pseudofermion deviational numbers and occupational 
configurations. The mathematical details necessary for an exact computation of 
the spectral function will then be presented in subsequent sections. For reasons 
apparent in section ()3.1.5j) . the following theory will not include any finite num- 
bers of — ^ Yang holons or HL spinous, to the contrary of the theory developed in 
Refs. |13nj and |131j . This will simplify the definitions of the relevant processes 
and subspaces introduced in this section. 

In section ()2.3.2p the PS subspace was introduced. This subspace is spanned 
by the ground state and all states generated from it by a finite number of electronic 
processes, i.e. by finite deviations AA^ and A(A^| Now, due to Eqs. ()2.98|1 

and ()2.100|) . these quantities are uniquely expressible in terms of deviations of 
pseudofermions, —\ Yang holons, and —\ HL spinous, respectively: 

CO 

A AT = AA^co + 2L^_i + 2 i^N^y 

oo 

A(A^I - Ni) = ANco - 2L^ _i - 2AA^,i - 2 ^ uN^u (3.12) 

i/=2 

which means that we can equivalently say that the PS is spanned by the ground 
state and all excited energy eigenstates with finite deviations of cO and si pseud- 
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ofermions, with finite (or zero) numbers of cv {u > 1) and su {u > 2) pseud- 
ofermions and witli finite (or zero) numbers of — | Yang liolons and — | HL 
spinous, respectively. One sfiould note tliat 1) tlie ground state, labeled \GS), is 
void of cz/ (z/ > 1) and su {u > 2) pseudofermions as well as of — | Yang ho Ions 
and —\ HL spinous, respectively, 2) apart from pure creation and annihilation 
of pseudofermions, the generation of the excited states which span the PS also 
involves a finite number of particle-hole processes in the av = cO, si bands, 3) the 
au 7^ cO, si branches have no ground state pseudofermion occupancy and thus 
do not have any Fermi points. However in these cases, the limiting canonical 
momentum values for the au effective Brillouin zone play the role of the Fermi 
points. To summarize from section ()2.2.5|1 : 
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(3.14) 

In the following, let the index l = ± denote the left Fermi point and/or 
shake-up discrete momentum shift of pseudofermions towards smaller canonical 
momentum values {l = — ) and the right Fermi point and/or shake-up discrete 
momentum shift of pseudofermions towards larger canonical momentum values 
[l = +). Let us now count the number of pseudofermions created and annihilated 
at the Fermi points {au = cO, si) and at the limiting canonical momentum values 
of the effective Brillouin zone {au ^ cO, si) on the one hand, and the number 
of pseudofermions created and annihilated away from these points on the other. 
The principal reason for this division is due to the fact that in the continuum 
momentum limit, the Fermi seas become compact since the momentum spacing 
(27r/L) 0. Hence a non zero phase shift Q^uiQ) inside the Fermi sea can not 
be detected due to the uniform occupation of pseudofermions. 

However, the situation is different with the outmost canonical momentum val- 
ues Lqpau that become shifted to a value for which there are no occupancies on 
the positive (t = -|-) and the negative [l = — ) side of that canonical momentum 
value. This means that the value of Q%j,{iqFau) defines the new canonical Fermi 
points, which are unique for each value of (U/t), n and the transition in con- 
sideration. Thus, creating or annihilating pseudofermions at their Fermi points 
should contribute much more to the dynamics than the corresponding actions on 
pseudofermions inside the Fermi sea. 
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Let the number of av — cO, si pseudofermions created or annihilated at the 
positive {i — +) and negative (t = — ) Fermi points respectively, be denoted by 
^■^avL- For av 7^ cO, si the definition is the same but with "Fermi point" replaced 
by "limiting effective Brillouin zone canonical momentum value" (however, we will 
stick to "Fermi point" even when av 7^ cO, si, keeping in mind that we actually 
refer to the limiting effective Brillouin zone canonical momentum values). If we 
add the extra contribution from the shake-up effect, we obtain a number AN^^ ^ = 
^^aft + ^Q^avl'^'^ of 0;^ pseudofermions at the positive {u — +) and negative 
(t = — ) Fermi points respectively. This "half particle addition" reflects the shift 
from integers or half-odd integers to half-odd integers or integers, respectively, 
of the quantum numbers introduced by the Takahashi string hypothesis. Thus, 
the total number deviation of av pseudofermions at the Fermi points is AN^^ = 
AN^^j^ + AN^^_. Similarily, we define the au pseudofermion current deviation 
AJ^^ logically as the difference between the number of au pseudofermions created 
or annihilated at the right and the left Fermi points respectively, i.e. 2AJ^^^ — 
^^au,+ - ^^au,~- We thus obtain the following 

AA^-^ = AN''' j_ + AN'' 
2AJ^ = AN'' ^ - AN'' 

The corresponding number deviation of au pseudofermions created or annihilated 
away from the right and the left Fermi points respectively, is denoted AN^f . 

These difi^crcnt types of dcviational numbers correspond to different types of 
ground state — > final state processes. We will classify these processes as A, B and 
C respectively, according to: 

• A: Creation or annihilation of av pseudofermions away from the Fermi 
points {qv = cO, si) or the limiting canonical momentum values for the 
effective Brillouin zone {ai> 7^ cO,sl). This is a finite energy and finite 
momentum process, affecting the number AN^f . 

• B: Creation or annihilation of av pseudofermions at the Fermi points {av = 
cO, si) or the limiting canonical momentum values for the effective Brillouin 
zone {au ^ cO, si). For the au — cO, si branches, this is a zero energy and 
finite momentum process affecting the number AN^^. It transforms the 
cO and the si densely packed ground state configurations into excited state 
densely packed configurations. 

• C: Small momentum and low energy particle-hole processes near the av — 
cO, si left and right Fermi points respectively, relative to the densely packed 
configurations obtained through processes B. For these processes, we will 
assume that the cO and the si pseudofermions disperse linearly. 



ANi 



AN^ 

^^Jl + (3.15) 
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At a later stage, we will use this classification of the pseudofermion processes 
when expressing the (rotated) electronic operators in terms of pseudofermionic 
operators. In Ref. |131j there is a similar classification of the different types of 
processes, however also including the numbers L^_i. For reasons apparent in 
section we do not need to include these numbers here. Also, we will not 

include particle-hole processes which are not in the vicinity of the Fermi points, 
due to the very small effect these excitations have on the one electron spectral 
weight pIK] . 

As the ground state is well defined in terms of occupational numbers of au 
pseudofermions, the excited energy eigenstates are well defined in terms of de- 
viational numbers of these pseudofermions. Indeed, each excited state is char- 
acterized by a number A^°^ + ANau of pseudofermions, as well as by a number 
L„ _i of Yang holons {a = c) and HL spinous {a = s). For each combination of 
these numbers, there exists then a subspace of many states all with the same de- 
viational numbers, but with different canonical momentum dependent occupancy 
configurations. This is particularly evident when considering the au = cO, si 
particle-hole processes C, which for each set of numbers {^Nau} au=co si corre- 
spond to many different eigenstates of the model with each of these having a 
different configuration of cO and si pseudofermions, as given by ANau{(l) where 

The subspace of states with fixed deviational numbers AN^i, for all au branches 
is called "cO pseudofermion, holon and spinon ensemble subspace", abbreviated 
" CPHS ensemble subspace" , in Refs. |13Uj and |131j , if we also include fixed (pos- 
sibly nonzero) numbers L^_i of Yang holons {a = c) and HL spinous {a = s). 

Note that there can be many states characterized by the same values of the 
deviational numbers AN^i,, but with different values of the deviational numbers 
AN^^^ (for au = cO, si and l = ±). Moreover, we note that for each state with a 
fixed number of the following deviations and numbers: ANau, AN^^^ and L^_i 
(for au = cO, si and t = ±), there exists a subspace with a total number of A^^^ 
particle-hole pairs for au = cO, si, due to the C processes. A typical element of 
this subspace contains one specific particle-hole configuration of number of 
cO pseudofermion pairs and N^i number of si pseudofermion pairs, with energy 
and momentum as specified by the numbers mau,i. of Eq. ()3.1fi|) . Therefore, we 
will see that the summation over the particle-hole towers of states will reduce to 
the summation over the integer numbers mau,i- 

The total number of subspaces here considered is less than what is being 
considered in Refs. [130J and jl31j . However, since we only study one electron 
spectral functions, to the contrary of the studies of these references, we can afford 
to simplify the subspace descriptions and minimize the total number of subspaces 
needed. 
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3.1.4 Energy, momentum and number deviations 



As before, we will not bother with finite occupancies of _i, due to the find- 
ings of section ()3.1.5p . which will simplify the expressions for the energy and 
momentum deviations. 

The energy deviations AEau = A_Eq,,^(A) + AEau{C) for au = cO, si give 
the total energy deviation from the ground state due to the A and C processes, 
respectively, and similarily AEau = AEc,u{^) for au ^ cO, si branches. We 
remind ourselves that the B processes are zero energy processes. After defining the 
Fourier transform of the electronic operator inside the spectral function, we will 
obtain another delta function for the momenta, of the form 6k,iAP, where AP = 
APa^. Here we have AP^^ = APa^(A) + AP„^(B) + AP„^(C) for ai^ = cO, si 
and AP„^ = APc,^(A)+AP^^(B) for au ^ cO, si (however note that AP,^(B) = 0, 
for > 2). The above energy spectra are given in terms of the energy bands in Eq. 
()2.7H|) together with the pseudofermion number deviations under consideration. 
We thus obtain a connection between the variables of the spectral function and 
the energies and momentum occupancies of our pseudofermions: 
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(az/ = cO, si) 



where the expressions for AP^o (B) and AP^i (B) stem from the effect accounted 
for in section ()2.3.6|) . namely that the au ^ cO, si pseudofermions at the hmiting 
canonical momentum values are felt by the au = cO, si pseudofermions as z/they 
were scattering centers at the cO and si Fermi points, respectively (for an exact 
derivation, see Ref. |13Uj ). Assuming a given momentum and energy reached by 
the A and B processes, mau,L measures the number of momentum steps in units 
of (27r/L), to the left {l = — ) or to the right {l = +), of such a momentum value. 
This simulates the particle-hole C processes in the vicinity of the au = cO, si 
Fermi points. There are many particle-hole processes contributing to one number 
(this description of the particle-hole excitations coincides with that of Ref. 



The number deviations are expressed as ^-functions at the corresponding 
canonical momentum values. However, in order to achieve the right dimension 
of our expressions, we need to remember that the quantities that have to equal 
each other in order for the (5-function to contribute, are the quantum numbers of 
the Takahashi equations, say / and /'. When going to the continuous system, we 
then have that: 



where L = aNa- 

In this fashion, we have for the various number deviations introduced in the 
sections above, that 



IHO!). 




iVa> 1 



(3.17) 



AiV„, = AiV^y + A<, 
AiV„.(g) = AiVr (g) + A<(g) + ^N^^q) 





i=l 




i=± i,=l 



where Np'I = N'^_^ + N'^_. It is these deviations that enter expressions like Eq. 
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(I2.61|l and Eq. ()2.67|1 in the evaluation of the energy deviations and phase shifts, 
respectively. 

3.1.5 Yang holons and HL spinons 

Having defined the relevant pseudofermionic processes and subspaces, we need 
to clarify how the rotated electronic creation/annihilation operator inside the 
spectral function relates to the pseudofermions. The subsequent analysis is further 
explained in Refs. [HS] |l<'^()j even though here we can profit from the simpler 

case of having only one electronic creation or annihilation operator. This fact will 
simplify our studies and we will not need to derive selection rules connecting the 
operators to the total number of L^_i Yang holons {a = c) and HL spinons 
(a = s) that can be created (these selection rules become trivial in the one 
electron spectral function case). 

To start with, let us note that the operator c^-^ stands for 4 different operators: 
Cj^, Cji^ and Cj^ respectively, where j = 0, 1, 2, . . . , Na—^-- However, considering 
our specific ground state and that we do not consider the UHB, we introduce the 
projection operators 1 — {c^j^Cj^j) = 1 — fija^ according to: 

RHB Cj^ Cj„ (1 - hj^) 

LHB cl^5li^-n,,) (3.19) 

where a = —a and j = in the spectral expressions of this chapter and of chapter 

We have many times claimed that we do not need to consider final states with 
finite occupancies of — | Yang holons and — | HL spinons. In the following, we 
will motivate this claim. 

Due to not considering the UHB, we note instantly that we can never create 
any L^_i Yang holons nor any cu {u > 1) pseudofermions, L^_i = N^u = 
always. Moreover, the total number of — | HL spinons can never exceed one. 
This is easily seen by the following two considerations: First, it is impossible to 
form a — I HL spinon either when annihilating a |-spin rotated electron (the only 
possibility is that the |-spin rotated electron came from a si pseudofermion, thus 
leaving an unpaired +| HL spinon) or when creating a |-spin rotated electron (of 
the simple reason that there is no combination of quantum objects that would 
allow a formation of a — | HL spinon in this case), respectively. Second, we can 
create either zero or one — ^ HL spinon(s) when annihilating a j-spin rotated elec- 
tron or when creating a |-spin rotated electron, respectively. In the first case, we 
may (or may not) annihilate the j-spin rotated electron from a si pseudofermion, 
leading to a single unpaired — | HL spinon left in the system. On the other hand. 
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if the spin degrees of freedom of the |-spin rotated electron was a +| HL spinon, 
the resulting system will continue to be void of — | HL spinous. In the second 
case, the spin part of the created |-spin rotated electron can either couple with 
a si pseudofermion (thus decreasing the number of +| HL spinous by one), or it 
can remain uncoupled in the system, giving rise to one — | HL spinon. 

We must remind ourselves here that the Yang holons and HL spinous are 
quantum objects that are invariant under the electron - rotated electron unitary 
operator V{\J/t). Moreover, their energy and momentum values remain constant 
during any ground state final state transition. As pointed out in section ()2.3.6p . 
the Yang holons and HL spinous are neither scatterers nor scattering centers. This 
can easily be seen since they do not suffer any phase shifts under an arbitrary 
ground state — *• final state transition. It follows that these quantum objects do 
not affect the dynamics of the model. 

In the case of the HL spinous, there are some straightforward estimates one 
can do comparing the spectral weight between final states with _i = 1 to final 
states with L^_i = 0, respectively. Let us first form the candidate final state: 



\f) = Si\LWS) (3.20) 

where S"^ brings the LWS state up one notch on the LWS — » HWS ladder. By 
direct evaluation of the norm, we obtain 

(/I/) = {LWS\SJi\LWS) = {LWS\ [s^,Si\ \LWS) - 

-{LWS\SiSa\LWS) = -2{LWS\S'JLWS) = L„ (3.21) 

by using the commutation relations of the SU(2) algebras and by noting that 
Sa\LWS) = 0. This means that a proper normalized state which is not a LWS 
but with one L^ _i occupancy, is given by 

I/) = -^Si\LWS) (3.22) 



Now, by using the explicit form for Ss, Eq. ()2.4j) . together with the usual fermionic 
anticommutation relations of the rotated electrons, we see after some algebra that 



5,4, - 4,5, = -a 4. + = (3.23) 

defining the operator Qs- Remember that according to our notation convention, 
the symbol / stands for creation/annihilation for the operator c^-^, as well as the 
numerical values ±1, as well as the signs ±. 
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For a final state with _i = 1, we now have for a typical matrix overlap that 
\{f\cl,\GS)\ = ^\{LWS\Sscl\GS)\ = 



1 



{LWS\cl^Ss\GS) - a{LWS\gs\GS)\ (3.24) 



where the first term of the last line is always zero since the ground state is a LWS 
(note that we follow previous considerations of letting j = whenever inside a 
matrix overlap expression). Since we will study systems where the magnetization 
ma — > 0, it does not matter whether we choose to study matrix elements with 
cr =t= + or a =|= -. 

This means that by choosing a =| for the RHB and a =t for the LHB, we 
obtain final states which carry no spectral weight if _ i = 1. 

We hence conclude that by choosing the following rotated electronic operators 
for the RHB and the LHB, respectively: 



RHB Cj^i (1 — rij i 



LHB d], (1 - n,, J (3.25) 



we obtain final states completely void of any occupancies of the — | Yang holons 
and — I HL spinous, respectively. Thus, from now on, these quantum objects will 
not enter in the following analysis of the one electron spectral function. 



3.1.6 Restricted subspace approximation 

By restricting ourselves to excited state subspaces such that L^_i = Ncu = 
(z/ > 1) the deviational expressions for AA^ and ^{N^ — N^) given in Eqs. ()3.12|) 
can be simplified. Moreover, due to the studies of Ref. [OSl, final states with finite 
occupancies of Ng,^ (z/ > 3) contribute very marginally to the spectral function. 
In this reference, the following was found for the thermodynamic limit at zero 
magnetization, and for values of the filling < na < 1 and arbitrary values of 
{U/t): Final states with Ngi number of si pseudofermions and with Ng^ = 
for z/ > 2, generate approximatively 94% of the total LHB spectral weight and 
98% of the total RHB spectral weight, respectively. Moreover, final states with 
Nsi number of si pseudofermions and A^s2 > number of s2 pseudofermions, 
but with = for z/ > 3, generate at most approximatively 6% of the total 
LHB spectral weight, and at most approximatively 2% of the total RHB spectral 
weight, respectively. Thus, with this motivation, Eqs. ()3.12|) become: 

AN = AN,o 

A{N^ - Ni) = AN^o - 2ANsi - 4AAr,2 (3.26) 
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However, these are not the only relevant deviations to study. To the con- 
trary of systems with free fermions, or even systems of spinless fermions in the 
{U /t) ^ oo Hubbard model, the total number of discrete momentum values A*"*^ 
is in general not constant, but deviates from their respective ground state values 
according to the specific transition in consideration. We have: 

N* =N -\- 

N* =N*''^ + AN* (3.27) 
AN* =AN -\- AN'^ 

where contra-intuitively ANav 7^ ~^^au general (however, for cO we will 
indeed always have that AA^*q = 0). From section ()2.2.5|) we have that N'^^ is 
the corresponding ground state number, given for the different branches by 

N*f = N^ (3.28) 
N*J = N^-N^ 



as well as the number of holes for each branch 



N^, = N,o-2N,,-2N,2 (3.29) 
N^, = N,o-2N,,-AN,2 



from which it is easily deduced that 

AiV:o = 



Vo ~ " 

an:, = AiVco - AiV.i - 2AiV,2 
AN:^ = AN,o - 2 AiV,i - 3 AiV,2 



(3.30) 



Vs2 

Note that the last deviation also can be expressed as 

AN:^ = ANs2 + A{N^ - N^) (3.31) 

which expresses the fact that in the zero magnetization limit, the entire band 
shrinks as ma vanishes. In the zero magnetization limit, the entire s2 band is 
nonexistent for the ground state since A^*2 — (and similarly for su bands with 
z/ > 3). However, if we create one s2 pseudofermion, there appears a canonical 
momentum band with a single discrete value, in order to accommodate for this 
quantum object. Thus in this case, we will have one static s2 pseudofermion, 
with zero energy and with zero canonical momentum. 
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The last deviational numbers we have to consider regards the eventual con- 
tribution from the shake-up effect, as given by Eqs. ()2.125|1 . In our case, the 
shake-up effect will contribute for the specific branch, if the following correspond- 
ing deviations are odd: 

cO : AA^.i + ANs2 

si : AN,i + ANco (3.32) 

s2 : AN,2 + AN,o 

Due to this consideration and the fact that states with many s2 pseud- 
ofermions are highly unlikely (i.e. produce negligible spectral weight), we will 
confine our final states to having maximum one s2 pseudofermion. Due to the 
findings of Ref . jH^ , we will label the final states as " Basic" or " Exotic" in ac- 
cordance with how much the corresponding final state contributes to the sum 
rule. In the following, if no s2 pseudofermions are created, the corresponding s2 
deviations will not be accounted for. Also, the approximative percentage with 
which the transition contribute to the total sum rule is given after the — >■ symbol 
in the heading of the each table. 
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RHB Basic 



AN = -1 A(A^^ - Ni) = +1 

AN^ = -1 ANsi = -1 

AN^^ = +1 

AiV.i = -l =^ iQcol^TT 

A7V,i + ATV^o = -2 ^ |Q°il=0 



RHB ExoUc 2% 

AA^ = -1 A(iV| - iV|) = +1 

AN^ = -1 AAT^i = -3 AN,2 = +1 

AAT^^i = +3 



ANsi + ANs2 


= -2 


=^ m 


= 


AN,i + AN,o 


= -4 


^ m 


= 


ANs2 + AiVeo 


= 


^ m 


= 



LHB Basic 94% 

AA^ = +1 A{Ni; - N^) = +1 

AiVeo = +1 AN,i = 
AA^,'\ = +1 

AiV,i=0 =^ \Q%\=Q 

A7V^i + A7V,o = +l =^ |Q°il=7r 



LHB Exotic 6% 

AN = +1 A(Ar^ - Ni) = +1 

A7V,o = +1 A7V,i = -2 AN,2 = +1 

AA^^'\ = +3 



AA^.i + AA^,2 = 


- -1 


^ IQ°ol 


— TT 


AA^.i + AA^co = 


= -1 




= TT 


AAr^2 + AAT^o = 


= +2 


=^ m 


= 



In similar fashion, for example, can the most relevant UHB transitions be 
classified, by keeping a non zero number Nd in the expressions for the deviations. 
Note that AA^^i 7^ — AA^^^^ for both of the LHB transitions. Here, this means that 
the transition induces one extra si pseudofermionic hole in the si band. 
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In the following, if not stated otherwise, we will only consider the "Basic" 
transitions, which only involves finite pseudofermion deviations in the cO and the 
si bands. In other words, all other branches will be assumed to be completely 
void of pseudofermions. 



3.2 State dependent dynamics 

3.2.1 Scattering phase shifts: particle-hole processes (C) 

Using the deviational numbers defined in section ()3.1.6|) together with Eq. ()3.18p . 
we can calculate the scattering phase shifts Q%{q) and Qfi{q) by use of the 
quantity defined in Eq. ()2.fi7|l divided by L, which in the continuous momentum 
representation becomes: 

Qtiq) = LJ2 ^au,a'u'iq,q')AN^'.'iq') (3.33) 

We can now specify the scattering phase shifts for each transition and for each 
process. Schematically, we would then have Qauil) = Qtuily-^) + Qtuil^^) + 
Qauily C) for each transition. However, as we will see by the following analysis, we 
will always have that Q%^{q, C) = independently of av and the transition under 
consideration, due to a pairwise cancellation of the phase shift of each "particle" 
and "hole" pair. To show this, we must recall that the C processes are defined 
in the vicinity of the Fermi points only, where the linearization of the dispersion 
relations remains a valid approximation. This means that the momentum values 
for the particle (gp) and the hole {qh) are only a distance 2t{J/L apart, where J 
is a finite number. We thus have for these processes a canonical momentum shift 
given by 

dq' $,,,„,,.(g,g')AiV^t(9') = 

f"'"' if, I 

= 2tt dq ^av,a'u'{,q, q) \5{.q' - q-p) - Hq - Qh) 



27r( ^cv,cc'v'{q,qr) - ^au,a'u'{q,qh) j = (3.34) 

2nj d^ au.a'u' iq,q') 



2tt[ ^au,a'u'{q,qh) + 



2n 



L dq' 

2txJ d^au.a'u'{q,q') 



dq' 



— ^au,a'u'{q,qh, 

g'=ih 

= 0{1/L^) 
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We remember that in our pseudofermion theory, discrete momentum contri- 
butions of order {1/Ly do not have any physical relevance for j > 2. Thus, such 
a canonical momentum shift is to be treated as equalling zero exactly. It is by 
this consideration that we claim that all the particle-hole processes belonging to 
the same tower of states all share the same phase shift. 



3.2.2 Relative spectral weights: Tower of states 

The property of the scattering phase shifts derived in the previous section allows 
us to treat the spectral function in two "steps": the first step due to processes 
A and B, respectively, producing a spectral weight at (/cqjI^o); given by a lowest 
peak weight Aaif^. The second step, due to processes C, produces a spectral 
weight with energy given by ujco + i^si, where Uau = (27r/L) '^^^^Vau^au,L and 
with momentum given by kco + kgi, where k^u = (2vr/L) Luiau,!., respectively. 

Thus, such A, B and C processes contribute to the weight at the point (/cq + kco + 

The final weight is a convolution of the weights for au = cO and au = si, 
respectively, where mau,i. is the number of momenta steps, measured in units of 
(27r/L), to the left {l = — ) or the right {t = +) of the {kQ,Uo) point. The 
superscript (0,0) refers to {niau- = 0,mai,,^ = 0). This procedure is inspired by 
that of Refs. [111-110], where the same problem is studied in the (f//t) ^ oo limit. 
The spinless fermions used to describe the model in these references correspond to 
the cO pseudofermions in the arbitrary [U/t) model. The "spinous" on the other 
hand are carried over from the ID Heisenberg model, which cannot be done in 
the arbitrary [U/t) case. Indeed, the si pseudofermion should not be compared 
to the notion of "spinous" that is used in Refs. [Zll-[in]- In these references, the 
spinous have a spin projected value of ±|, whilst in our pseudofermion theory for 
arbitrary (U/t), the corresponding si pseudofermion is a two-spinon object with 
zero spin projection. 

As distinguished by the numbers (jnau-,''TT'au,+), we then obtain a tower 
of states where mau,+ + fnau- is proportional to the energy of the particular 
particle-hole process and where mau,+ — fTT-au,- is proportional to the momentum 
of the same process. The number of particle-hole processes contributing to the 
{mau-,rnau,+) point grows exponentially as we build the tower. Due to conven- 
tion, we go to successively more negative energies whenever in the RHB and to 
successively more positive energies whenever in the LHB, as we build the tower 
of particle-hole states emerging from momentum k and energy u. 

As in, for example, Ref. [10], we now define the relative spectral weight 
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according to 

)AT (3.35) 



In the cited reference only one one spinless fermion and one spinon is created 
or annihilated, during the one electron addition or removal process. Thus the 
explicit form for aai,{Tnai,^-,mai,^+) is simpler than what we will need here, where 
in general a number of ANau pseudofermions is created or annihilated for each au 
branch. This means that wherever in Ref . jHO] , the excited energy eigenstate has 
+ / spinless fermions or spinous (for / = ±1), we in our case have N^^ + AN^u 
av pseudofermions. In the following analysis, we will explicitly calculate aau{l, 0) 
and aa!^(0, 1) which then easily combines to form a^ui)-, 1). For reasons of clarity, 
we will consider excited energy eigenstates resulting from processes B only, such 
that AN^J = 0, i.e. such that there are no au pseudofermions away from 
the densely packed excited Fermi sea. This means that we will consider AN^u = 
^■^au — ^-^au + ^ ^-^au -■ Furthermore, we will consider only one pseudofermion 
taking part in the most basic particle-hole excitation process: the one in which 
the hole momentum is at the shifted (excited state) t Fermi point and the particle 
momentum is just one step away from this point, i.e. % = % + i{2T[ j L). Adapted 
to our notation the generalized expression for the relative weights, for one particle- 
hole pair, is then 



n sin^(^) nsin^(^) 
%) = '^\o.^.. ^^^^^T (3-36) 



sm 



' (^) n sin^ (i^) n sin^ (^) 



i=i i=i 



Thus, this relative weight is only valid for one particle-hole pair. Moreover, the 
notation aaui^h, %) is allowed here due to the simplicity of the particle-hole pro- 
cess in consideration. For any other particle-hole process, we would have to write 
O'avi'n^av- 1 fTiav,+) siuce the canonical momentum values of each particle-hole pair 
uniquely define two integers {niau-^'mau,^) whilst the converse is not true. As it 
stands, this expression is only valid in the finite system, the continuous momen- 
tum limit will be taken at a later stage. The equivalent expressions for higher 
numbers of particle-hole pairs can be found for {U/t) ^ cxd in Ref. [20] and for 
arbitrary [U/t] in Ref. [TTn] . 

Since we are considering densely packed excited state pseudofermion occupa- 
tion configurations, most of the factors above will cancel with each other, which 
can be easily seen by explicit calculation of a^uisih, %) using the following expres- 
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sions for the relevant canonical momenta: 

- _^ I QtiQh) ^ _ 2vr 27r . ^ . 

qp = qh + i^j- ^ ip = ih + i^ (3.37) 

However, this fact is only due to some quite obvious approximations, based on 
the difference 

5^_,-..|(/^_/.+ «M_l«tM) (3.38) 

We note that in the second term above, the difference of the scattering phase 
shifts, is always bounded and small due to the boundedness of the ^au,a'u'{q,q') 
functions, as compared to the difference \Ij — Ii\ = 1, 2, . . . For example, consider- 
ing the "Basic" excitations of section ()3.1.6p . we have that (Qaui^j) ~ Qauiqi)) /^^r 
attains a typical maximum value between 1 and 2 when qj and are on oppo- 
site sides of the cO or si Fermi sea ^98j, and is always smaller than this value 
for any other pair of Qj and g^. Moreover, since the scattering phase shifts are 
in general continuous functions of the momentum, when the two canonical mo- 
menta are close to each other, the difference of their scattering phase shifts be- 
comes negligible. This means that the difference between the two scattering phase 
shifts can always be neglected. Moreover, since the particle and the hole canon- 
ical momenta both are in the vicinity of the same Fermi point, we have that 
QtiAlh) ~ QtuiQp) ~ Qtui'^lFau) ■ By following this scheme, we obtain for the 
relative weights that: 

K. .in' (i (I, - 4) - 2%^) „„2 _ ;^)) 

"°"'*'*'"Msi„^(i(/,-«-2Lif«l) l\ Bin' (id, -4)) 

(3.39) 

which by explicit investigation of the factors allows for the above mentioned 
cancellation. After this cancellation, what remains is 



sin^ (i^^ + 9lJm^^ (vrn + ^ ^ 
" sin^(vrn + ^ + -^-fr-))sin^(S) " ^''"^^ 
which in the large system limit goes as 

aUciH, ^ (.A<, + 9U;^^^ ' . (3.41) 
1 + I (a<_, - ^^"ff-^ ) cotvrn + O (l/L^) 
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Even though strictly speaking this formula for the relative weight is only valid for 
the lowest particle hole excitations aau{^,0) and 0^2/(0,1), we notice that there 
is a leading order term which does not depend on the system size. This leading 
order term will remain untouched as we go to higher particle-hole processes, as 
relative weights for successively higher particle-hole excitations only contribute 
to the (l/L) term, which is verifiable by suitable modification of Eq. ()3.39|) . 



3.2.3 Relative spectral weights: closed form expressions 

Generally, we have particle-hole processes at both sides of the Fermi sea at the 
same time. This means that in the expressions for the relative weights of the 
preceding section, we will have to consider many particle-hole pairs {qhiApi)i 
(g/j2,^p2), ••• , {'ihis[pi^ 1 'ipisfp'^) ■ The expression for the relative weights for more 
than one particle-hole pair are given in Refs. [HD] {U /t — »• oo) and |131j (arbitrary 
U/t). For example, to continue our study for mav,± = 1 of the previous section, 
it is easily shown that the relative spectral weight of two particle-hole processes 
can be written as (adapted to our notation): 

a„.(l, 1) = a..(0, 1) a„.(l, 0) ^-^2 ^-^2 (\-. r^ (3-42) 

which by introduction of the quantum numbers in the same way as before becomes 



aaz.(l, 1) = aauiS^, l)aai.(l, 0) ■ 



1 + C(1/L' 



(3.43) 



where the exact cancellation of the {1/ L) terms is a result of the mav,L = 1 special 
case. Exactly the same procedure can be repeated for successively higher particle- 
hole excitations, i.e. for numbers mau,± = 1,2,..., with successive contributions 
to the first order correction term. In general as mau,± increases we will have more 
and more factors that do not cancel in the expression for the relative weight. The 
number of factors that do not cancel grows exponentially with mai,,±, with an 
additional exponential increase of the factors that mixes the canonical momenta 
for the holes and the particles of the left and the right Fermi points, respectively. 
As shown above, this effect is not present in our mau,± = 1 example. The expres- 
sions involving many particle-hole pairs are very involved and an exact derivation 
of the various cases for increasing mau,± will not be given here since it is more 
confusing than enlightening. However the method used is exactly the same as in 
the example shown above. Due to the number of non-cancelling factors present, 
we find that we have to relax the 0{1/L) correction term slightly, according to 



0{\nL/L) 



(3.44) 
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which nevertheless allows for a complete description of the spectral weights of the 
processes C in the thermodynamic limit. 

In section ()3.1.3p it was mentioned that the dynamics of the system should, in 
the continuum limit, depend heavily on the scattering phase shifts at the Fermi 
points. Consider now a ground state —>■ final state transition. The pseudofermions 
that are inside the Fermi sea become shifted from one canonical momentum value 
to another that is also inside the Fermi sea, and that thus was also occupied in 
the original ground state. Thus, the corresponding phase shift becomes like a 
grain of salt in a bowl of water that technically is visible during the transition 
itself, but that in the final state becomes "dissolved" in the filled densely packed 
excited energy eigenstate. Note however that this is only true for the continu- 
ous momentum limit. The shifted Fermi points on the other hand, define the 
boundaries for the new Fermi sea. Because of this, the amount that the canonical 
momenta at the Fermi points have shifted, AqFau,L, is an important quantity for 
the description of the dynamics of the model. We can see for example in Eq. 
that 



which means that by defining 



(3.45) 



we have indeed 



aa.(0,l) = a„,,+(l) = 2A+, (3.47) 
a„,(l, 1) = a„,,_(l) a„,,+ (l) = (2A-J (2A+J 

for L — oo. The quantities 2A''^j^ will be very important in our subsequent study 
of this problem. As already mentioned, this analysis can be carried to higher 
orders, which give for example, 

a„.(2,0) = a„.,_(2) = — "^^^^"'^^^^ 



2 

2A+ ('2A+ + 1) 
aUO, 2) = a,.,+ (2) = \ "'^ ^ ' (3.48) 



and so forth, reaching the general result 



, . , 2A^,(2A^, + l)...(2A^, + m..,-l) 
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Algebraically, this can be rewritten as 




2A' -1 



(3.50) 



where the approximation is almost exact except for the first few positions in 

the tower of particle- hole states, provided that l/Y^^ ^ 0. For = there is 

no change in the canonical momentum value of the l Fermi point and hence there 
is no dynamical change between the ground state configuration and the excited 
energy eigenstate configurations of au pseudofermions at that Fermi point. 

Eq. ()3.50|) shows the power law type behavior of the spectral weight generated 
by the particle-hole processes, with exponent equal to {'2A''^^ — 1). This exponent 
obviously changes for each position in the {k, u) plane that the processes A and 
B brings the excitation to. Furthermore, the total spectral weight at a certain 
position in the tower of states will then become a summation of contributions 
from different particle-hole processes originating from neighboring points in the 
{k,uj) plane due to processes A and B. Therefore, the final exponent will be 
different from the exponent given here, as we will take into account contributions 
from many different overlapping towers of states. 

3.2.4 Scattering phase shifts: A and B processes 

In the previous section, we derived a closed form expression for the relative spec- 
tral weight for the particle-hole processes following a given A and B process. 
The spectral weight generated by the particle-hole processes is crucially depen- 
dent on the phase shift that the processes A and B produces. Since ANauig) = 
AN^jf{q) + AN^j^{q) for the A and B processes, we similarity have by definition 

Qtvi^QFav) = Qti^^\LqFau) + Qti^\l^qFau) ■ Wc have 



2n 



E W'^' (^^^-' ^^a'u'iq') + ^N^M) (3.51) 



where by the use of Eq. ()3.18j) we see that 




av,a'v' i^qpau, q') AN^,^,/,^^g^^g,)6 {\q'\ - qF,a'u') 




(3.52) 
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To this end, we will define some quantities i^ua'u' where j = 0, 1: 

^au,cx'u' = ^av.a'u' + ^au,a'u' {^QFavi QPa'u') + (" 1 V (l^QPau, —QPa'u') (3.53) 

which implies that ^au,a'u' {t'lFau, i-'^lFa'v') can be expressed as 

/ 1 

^au,a'u' {t'lFau, iFa'u') = 2^au,a'u' + 2^au,a'u' ~ ''^L,L'^au,a'u' (3.54) 

The introduction of these quantities simplify the expressions of the scattering 
phase shifts as we will see below. Together with the identity AN^,^, ^, = l'AJ^,^, + 
AN^,^,/2 we can substitute Eq. (ITMj) into Eq. (ITH^ to obtain 



27r ^ 2 

which obviously implies that 



SO that we finally obtain, by performing the a'v' = cO, si summation, 



(3.56) 



2 A' -ff^ AJ^+f^ A/^ + /^° ^^^^ I ,cO ^^ii I \LqFau 

(3.57) 

where generally 



^^iNF)f ^ 

'^^^^^^^ = J2-gn{AN^,^,) 



{l-qFau, q'l) + ^av,a'v'{LqFav, + • • • + 



+ ^au,a'u'{l-qFau, 



a' u' ' 



(3.58) 



is the scattering phase shift at the l Fermi point due to the created or annihilated 
a'u' pseudofermions with canonical momenta g^, • • q'l^iyNF >, for aV = cO, si. 
Note that for the "Basic" transitions of section ()3.1.6|) . we will always have that 
\AN^,I^,\ = 0,1. The quantities 2A''^j^ ultimately control the behavior of the 
spectral function in the {k,uj) plane. We will see later that different linear combi- 
nations of 2A^ and 2A^^ form the exponents of the one electron spectral function 
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for the excitations for the "Basic" transitions, which lead to a power-law behavior 
of the spectral function. This mixing of the various 2A^^'s is a result of the fact 
that several different processes contribute to the same vicinity of a typical point 
in the {k, u) plane, as is evident in the case of overlapping towers of states. 



3.2.5 Excited state characterization 

In section ()3.1.fi|l . we found that nearly all of the total spectral weight was dom- 
inated by two types of transitions for both the RHB and the LHB. These transi- 
tions were called "Basic" and "Exotic", respectively, due to their pseudofermionic 
content (as shown by the occupancy number deviations) and their contribution 
to the total sum rule. As mentioned in that section, and as implicitly assumed in 
the subsequent sections, we only consider the "Basic" transitions. These transi- 
tions are such that the only bands with finite occupancies are the cO and the si 
bands. They contribute to about 98% (RHB) and 94% (LHB) of the total sum 
rule, respectively. In this section we will express the momentum k and energy u, 
which are the variables of the spectral function, in terms of the pseudofermionic 
or the pseudofermionic hole canonical momenta. The procedure is actually quite 
straightforward, considering that we have already calculated all the relevant quan- 
tities needed. All that remains is to specify exactly how the specific transitions 
here considered affect these quantities. From section p.l.ip we have that the 
variables of the spectral function obey 



- Egs = 5Z 

au=cO,sl 



Pfi - Pgs 



ai'=cO,sl 



AK,(A)+ AE„,(C) 
AP„, (A) + AP,, (B) + AP,, (C) 



(3.59) 



Since we will only consider the "Basic" transitions we can simplify the expressions 
of Eq. ()3.1fij) according to: 



IAE = 1 [sgn (AATj^^) e,o(gco) + sgn (AN^/) e,i(g,i)] (3.60) 
ZAP = I [sgn (AiVj^^) g,o + sgn (AN^/) q,, + AkpAj[, + 2kF^Ajf,] 

where we have used the fact that we will never have more than one cO or si 
pseudofermion created or annihilated. Qco and g^i are the corresponding momen- 
tum values for the created or annihilated scattering centers, which should not 
be confused with the Fermi momentum qpau, nor with the limiting momentum 
value for the au effective Brillouin zone Judging from Eq. ()3.57p . we see 
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that we will have qualitatively different expressions for the quantities 2A^^ when 
pseudofermions or pseudofermion holes are created at any of the Fermi points on 

the one hand, and when they are not, on the other. This leads us to consider 
four different cases (where in the following "P" stands for those Pseudofermions 
created or annihilated away from the Fermi points): 

1. 2P contribution: Neither the cO nor the si pseudofermion or pseudofermion 
hole are created at any of the Fermi points. This contribution will lead to 
the overall "background" of the weight distribution of the spectral function, 
as both pseudofermions or pseudofermion holes are dispersive, leading to 
contributions over nearly the whole range of allowed k and lu values. 

2. s-branch (IP): The cO pseudofermion or pseudofermion hole is created at 
the left or the right cO Fermi point and the si pseudofermion hole is created 
away from any of the si Fermi points. This will lead to a line in the {k, u) 
plane, following the dispersion of the si pseudofermion hole. 

3. c-branch (IP): The si pseudofermion hole is created at the left or the right 
si Fermi point and the cO pseudofermion hole is created away from any of 
the cO Fermi points. This will lead to a line in the {k, co) plane, following 
the dispersion of the cO pseudofermion or pseudofermion hole. 

4. Fermi contribution (OP): Both pseudofermions or pseudofermion holes are 
created at their left or right Fermi points, respectively. This contribution 
leads to a spectral weight distribution in the vicinity of certain points in 
the (kjUj) plane. 

This distinction is merely due to the number of different ways that we can 
calculate spectral weight in (k, uj) domains topologically different from each other, 
as well as due to the observation that the numbers AN^^ and AJ^^ are different 
in each of the cases above. Before we move on, however, there are two special 
cases which will be important in the following. They are classified as 

1. 2P contribution - Border Lines: These line shapes share the definition of the 
general 2P contribution described above, but with the additional require- 
ment that the velocity of the non Fermi cO pseudofermion is equal to the 
velocity of the non Fermi si pseudofermion hole. This extra requirement 
confines the 2P spectral weight to certain fines in the (fc, a;) plane. We will 
see that the general expressions for the 2P contribution are singular as the 
two velocities become equal, giving rise to a divergent feature of the overall 
spectral function. 
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2. s-branch and c-branch - Luttinger contribution (IP): This "partial line 
shape" shares the definitions of the s-branch and the c-branch described 
above, with the difference that the dispersive si (s-branch) or cO (c-branch) 
pseudofermion has a finite but small energy. These excitations, which are 
close to the Fermi energy level, belong to the subspace of excitations usually 
described by the Luttinger liquid theory. We will see that excitations be- 
longing to the branch lines but in the region of very small energies will have 
different final expressions of the full spectral function, as compared to the 
expressions of the finite energy s-branch and the finite energy c-branch, re- 
spectively. The values of the critical exponents at low energies are obtained 
in this regime, and have been shown to reproduce known results obtained 
by conformal field theory 1^21 11211 1123 • 

To make the following analysis easier, we note that 

2AJl = A<;f_, - A<;f_ + (3.61) 

by using AiVf,^, = AN^^i^^ + iQlJ2n and where = QlJii = -1, 0, 1. 

Lastly, we will express the RHB and the LHB annihilation and creation rotated 
electronic operators in terms of pseudofermionic operators. The latter operators 
are defined in Eq. ()2.1()7j) . and denoted f^au and fg^au in the canonical momen- 
tum representation, and fj^^^^ and fj^au in the effective au lattice representation 
(where j is denotes the effective lattice site coordinate). The rotated electronic 
operators occurring in the expression for the spectral function are local opera- 
tors and can be uniquely expressed in terms of local pseudofermionic operators. 
Without further due, the conversion between the two representations results in 
the following leading order expressions for the au = cO, si pseudofermions: 

RHB Cj- 1 (1 - %,-f) ^ fj,cofj',si 

LHB (1 - J ^ 4o (3.62) 

where the cO pseudofermion is annihilated (RHB) or created (LHB) at position 
Xj = CLcoi and the si pseudofermion hole is created (RHB) at position Xji = a^ij', 
such that Xj ~ Xj' |13Uj . In the latter case, we also produce a si pseudofermion 
hole, however, this hole results from the emergence of one extra canonical mo- 
mentum value in the si momentum space, and not due to the destruction of 
a si pseudofermion, as confirmed by the analysis in section ()3.1.6p . The LHB 
transition is further described in section ()3.2.7p . 
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In Eq. ()3.62|1 . we use the approximation sign for the following reason: There 
should be a pre-factor on the right side of the equalities above, which in Ref. 
is denoted 1/Cj. For example, the strict operator equality for the dominant con- 
tribution of the LHB transition, is given by Eq. (47) of that reference. However, 
for our cases, this pre-factor reduces to a simple phase factor, since the quantities 
denoted Gc and Gj in this reference are both equal to one, as can be seen in Eq. 
(57) of that reference. Moreover, the argument of this phase factor is given by 
the discussion and the equalities on page 17 of that reference. Since this phase 
factor will not be important for the subsequent analysis, it is omitted here. 

A note on the effective av lattice: The lattice site index j is not arbitrary. 
Rather, a rotated electron being created or annihilated on a lattice site position j, 
corresponds in the cases considered here to creation or annihilation of av = cO, si 
pseudofermions on strictly defined au effective lattice site coordinates. These 
coordinates can be found in the discussion on page 12 of Ref. We find 

that the cO effective lattice site is the same as the rotated electronic lattice site, 
whilst the si effective lattice site j' equals the closest integer number to jn^ 
jn/2 as m — » 0. Note that Xj ~ Xj>, i.e. that the two pseudofermions are 
(approximatively) created or annihilated at the same spatial coordinate. 
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3.2.6 The RHB "Basic" transition 



The RHB " Basic" transition involves creation of one cO pseudofermion hole and 
one si pseudofermion hole. It is characterized by \Q%\ — n and IQ^il — and 
thus Aco = sgn((5co) and Asi = respectively. In the case of the "s- branch" and 
the "c-branch", we will need to specify at which Fermi point the cO and the si 
pseudofermion hole, respectively, is created. For this, we will define the quantity 

1. RHB 2P contribution: —2kp < Qco < ^kp — kp^ < g^i < kp^ 



ATVjVF ^ _i AiVfo = AJfo = Aeo/2 
i^N^/ = -1 /^Nf-y = AJjJ = 

gS^) {i2kp) /2n = Aeo (-1 + Co,co) /2 
QtP (''kpi) /2n = A,o Ci,co/2 

g*J^^) {L2kp) /2n = - (<l>,oM^2kp, q,o) + <^coA^2kp, q,,) 



Qsi i^^Fi) /27r = - ^siM^^Fi, qco) + ^si,si{^kpi, qsi) 





= eco(gco) + esi(gsi) 



k = qco + qsi - Ko^kp 
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2. RHB s-branch: Qco — ico^kp — kpi < qsi < kpi 
AiV,r = AN^ = -1 AJfo = (Aeo - .co)/2 



AN^/ = -1 ANf^ = AJfi = 



qS^) (62fc^) /27r = .(1 - ecVo)/2 + (Aco - ^co)(ec'o,co 
gff) (.^^J /2n = -.e°i,co/2 + (Aco - ^co)ei,co/2 
g*^^^) (i2A;^^) /27r = -$co,.i(i2A;f , 



27r 

OA' - I "^°l.cO (AcO - <'co)C\,cO g^^(^-^) (-^fcp^ ) 
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3. RHB c-branch: —2kF < ?co < ^kp qsi — hikpi 
ATVjVF ^ _i ^ ^jF ^ 

ANf/ = ANf^ = -1 AJfi = -t,i/2 

qS^) {i2k,.) /27r = -.e°o,.i/2 + Aeo(ecVco - l)/2 - ^.iec'o,.i/2 
QtP (ikj.^) /2n = .(1 - e°i,.i)/2 + AeoC\,co/2 + - ei,.i)/2 
q5^^) (t2A;p) /27r = -$eo,co(i2A;p, ^^o) 



27r 



" I 2 + 2 2 + 27r 

^ = eco(<Zco) k = QcO — \co2kp + tsi/Cpi 
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4. RHB Fermi contribution: Qco — ico^kp Qsi — tsikpi 

AiV,r = AiVe^ = -1 AJfo = (Aeo - .co)/2 
AN^/ = ANf^ = -1 AJfi = -i,i/2 

qS^) (.2^;^) /27r = .(1 - ^,00 - a,si)/2 + (Aco - ^co)(ec'o,co - l)/2 - ^.iec'o,.i/2 
QtP (ikp^) /2n = .(1 - e°i,eo - e°Mi)/2 + (Aco - ^00)^1,00/2 + ^.i(l - eMi)/2 
g^""^) {L2kp) /2n = 
^r^^^ i^kp,) /2n = 

oAt / ''(^cO.cO + CcO.sl) (AcO — ''co)CcO,cO , '•sl^cO.sl 

2A^o = I ■ ^ — + 



_ I ^(^°i,cO + C°i,sl) _ (AcO - <-co)Csl,c0 _^ ''slCl.sl 



O; = /C = isl^n ~ ('^cO - tco)2/CF 
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3.2.7 The LHB "Basic" transition 



The LHB "Basic" transition has a pecuharity that is absent in the RHB case, 
namely that even though no si pseudofermions are annihilated from the system, 
one si pseudofermion hole is being created. The created f-spin rotated electron 
opens up the canonical momentum space for the si pseudofermions, allowing 
one more discrete canonical momentum value in the si band, without changing 
the number of si pseudofermions present. Thus we have that AA^^i = and 
AA^^]^ = +1, which induces a shake-up in the si band, i.e. \Q%\ = and = vr 
and thus Aco = and A^i = sgn (Q'^i) respectively. The canonical momentum 
value of this si pseudofermion hole controls the dynamics of the si band, and we 
will thus obtain different dynamical descriptions of the spectral weight depending 
on whether this hole is created at any of the Fermi points or away from the Fermi 
points. Actually, the si pscudofermionic current will depend on where the si 
pseudofermion hole is created, which is also why the number A^i will not appear 
in any of the expressions below. The fact that IQ^il = tt does not mean that we 
first shake-up the si band and then create a si pseudofermion hole in it, rather, it 
is the emergence of one extra canonical momentum value, i.e. the increase in the 
total number of allowed canonical momentum points, which induces the shake-up 
effect. 

Thus, in spite of the creation of one si pseudofermion hole, the expression 
of the rotated electron creation operator in terms of pscudofermionic operators 
does not include a si pscudofermionic annihilation operator. However, the si 
pseudofermion hole that emerges, gives rise to a finite current if it appears on any 
of the si Fermi points (which for ma — > coincide with the effective Brillouin 
zone limits). 

The way to see this is that if the hole emerges inside the si Fermi sea, it 
"pushes" all the si pseudofermions on its left side a half step towards the left 
and all the si pseudofermions on its right side a half step towards the right. This 
means that we will have a "half-particle addition" at each Fermi point, implying 
that AN^i ~ i + I ~ -'- ^ven though AJj^ = 0. If the hole emerges at the c 
Fermi point then it will only have si pseudofermions on the right (t = — ) or on 
the left (i = -I-) side of it, and the above mentioned "pushing" will result in a 
global relocation of all si pseudofermions present in the band. In this case we 
have that AA^^ — (since the emergence of the hole at one of the Fermi points is 
cancelled by the appearance of another "pushed" si pseudofermion at the other 
Fermi point) and that 2AJ^ = —tgi = — sgn(gsi), where qsi is the momentum of 
the emerging si pseudofermion hole. With this peculiar effect in mind, we obtain 
the following characterization of the LHB " Basic" transition: 
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LHB 2P contribution: 2kF < \qco\ < — kpi < Qsi 



AN^/ = +1 AN^, = AJ[, = 



AN^/ = -1 ANf^ = 1 AJfi = 



qS^) {i2kj.) /27r = .e°o,.i/2 
gff) (Lk^,) /27r = - l)/2 

g^-^^ (i2A;^^) /27r = (^^ooA^'^kF, Qco) - ^coA^'^^f, Qsi) 
Qti^^^ (i^kFi) /27r = ^$a,co(''fcFi, Qd)) - ^siMi^^fi, Qsi) 

... AcVi , QT'\^^kF) \ 

" 1^ 2 + 27r j 

^^'^ " 1^ 2 + 27r 

= eco(9co) - k = qco-qsi 
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LHB s-branch: Qcq — idd^kp — kpi < qsi < kpi 

ANf/ = -1 ANf^ = +1 AJfi = 

qT i^'^kp) /27r = .(e°o,co + i%,s. - l)/2 + ^00(^^00 - 
gff ) /27r = .(C°i,,o + - l)/2 + icil.J'^ 

q2^^) (t2A;f ) /27r = -$eo,.i(i2A;f , 
Qfi^'^^ (t/cFi) /27r = -*^i,5i(tA;F;, g^i) 

_ I (Co.cO + CcO,sl) _^ ^coCcO,cO _^ QcO ^ {l^kp) 



27r 

^^si - I o + - + 27r 
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3. LHB c-branch: 2kF < \qco\ < Qsi — hikpi 

AiV,r = +1 AArf„ = AJ^ = 
AN^/ = ANf^ = AJfi = -i,i/2 

qS^) (62fc^) /27r = -i.iec'o,.i/2 
gff) (.^^J /27r = (1 - /2 

Qti^^'' {i^kpi) /27r = ^si,co{'^kFi, Qco) 
" 2 + 27r J 

~[ 2 ^ 27r 
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4. LHB Fermi contribution: Qco — ico^kp qsi — hikFi 
ATVjVF ^ AN[, = +1 AJfo = .,o/2 
ANf/ = ANf^ = AJfi = -i,i/2 

qS^) (.2^^) /27r = .(a,co - l)/2 + ^co(ec'o,co - l)/2 - ^.i^c^o, 
QtP (ikj.^) /2n = .e,co/2 + ^00^1,00/2 + ^.1(1 - eMi)/2 
qS^""^ (^2A;^) /27r = 
Q!f (^^Fi) /27r = 

^^^^ " 1, 2 + 2 2 
cu — k — icaikp — isikpi 
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Chapter 4 



The One Electron Spectral 
Function 



4.1 Basic Derivation 



4.1.1 Matrix elements and pseudofermion operators 

In this chapter, we express the state generators and the operators of the matrix 
elements of the spectral function in terms of pseudofermionic creation and anni- 
hilation operators. We remind ourselves that the operator expressions presented 
here are special cases of the ones calculated in Ref. In the following, we 

will only focus on the necessary ingredients of the one electron spectral function. 
Thus, for example, Fourier transforms of a product of an arbitrary number of 
pseudofermionic operators will not be formally performed here, but the proce- 
dure will be described in general terms and the final results presented. Moreover, 
in this section, we will temporarily revert to the momentum representation (which 
is also used in Ref. and thus write c|.^ for the rotated electron creation 

{I = +) or annihilation (/ = — ) operator. 

In section fj2.2.5p . we defined qpau as the positive ground state au = cO,sl 
Fermi momentum in the thermodynamic limit, qpco = "^kp and qpsi = kpi- These 
quantities are dubbed q^^^^ in Ref. |131j . First off, we will define the following 
quantities: 



(lFau,i — 



qFav,i — 



L-qPau + ^qFau,L 




+ 



L 



(4.1) 
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which are nothing but the deviations of the Fermi point momenta and canonical 
momenta, respectively, in the excited state configuration relative to that of the 
ground state. We notice that since the ground state consists of two densely packed 
minimum energy Fermi seas (i.e. with all the canonical momenta symmetrically 
distributed around zero), one for cO and one for si, it can be quite easily expressed 
in terms of cO and si pseudofermionic creation operators 

\GS)= n n/.uio) (4-2) 

au=cO,sl 1=1 

where |0) is the pseudofermion vacuum state. 

Now, a typical matrix element occurring in the definition of the spectral func- 
tion of Eq. p.ip . is of the form {fi\c[JGS) . Here 1/;) denotes an excited energy 
eigenstate and \GS) is the ground state of Eq. ()4.2|) . Due to the findings of 
chapter (jH}, we can rewrite the rotated electron creation {I = +) or annihilation 
(/ = — ) operator, in terms of operators that creates and/or annihilates au = cO, si 
pseudofermions according to the A and the B processes, respectively. Thus, let 
us write = c^fco-^Bo-) where c^^^ is associated with the A processes, and c^^ 
with the B processes. The latter operator creates a number X]a!/=co si ^-^au 
pseudofermions at the av = cO, si Fermi points, producing a densely packed 
state denoted \GS) = c'^^\GS). This state can be expressed solely in terms of 
pseudofermion creation operators as: 

i^)= n n fLM (4.3) 

au=cO,sl i=l 

where qi = -qpau + ^QPau- and Qno^+anp, = Qfou^ + ^qFau,+- The matrix 
element {fi\c''i^JGS) can then be written as {fi\c''^i^JGS). 

Consider now the state {fi\ = {fi\c''^ka- Upon acting onto (//|, the op- 
erator cj^ removes the finite energy pseudofermions or pseudofermion holes 
created under the \GS) \fi) transition, due to the A processes, and thus 
{fi\ = {fi\(^A,ka is also a densely packed state. The state \fi) can be generated 
from the pseudofermion vacuum, by acting with pseudofermion creation operators 
carrying canonical momentum values Qi = qi + Q%^{qi)/L: 

1/')= n n (4.4) 

ai'=cO,sl i=l 

Here, Q^j^{qi)/L is the scattering part of the canonical momentum shift of 
the excited energy eigenstate 1/;), relative to the corresponding ground state 



118 



discrete momentum g^. This means that after taking into account the A and the 
B processes, the matrix element (/i|c^ fco-|G'5') can be written as {fi\GS). 

These processes take the excitation to a certain energy and momentum in the 
(fc, u) plane, upon which the spectral weight is calculated. On top of this we 
build a tower of particle-hole states, generated by the application of the operator 
^l-hau oiito the state \fi). This operator reads 



^p-h,au W fqi,avfq[,au (4-5) 
i=l 

for a total number of N"^ particle-hole pairs in the vicinity of any of the two av 
Fermi points. 

In order to proceed, we focus now on the matrix element {fi\Fp_h^oLu\GS) . This 
quantity originates from the matrix element {fi\Fp_h,avc\„\GS) , that involves the 
final state with N'p^ particle-hole pairs, namely F^_^^J\fi). This state obeys the 
following equality: 

(/d^p-h,..5L|G'^) = {h\Fp^n,.AM\GS) (4.6) 

Importantly, we note that the excited states |/;) and F^_f^^J\fi) have the same 
canonical momentum shift Qau{,<ii)/ L = Q^^^/ L + Q^^{qi)/ L, for every g^, due to 
the findings of section 1)3.2.11) . 

Since Fp_h^ai/ involve pseudofermion operators with canonical momenta in the 
vicinity of the av = cO,sl Fermi points, whilst c^^^ can be expressed in terms 
of pseudofermion operators with canonical momenta away from the au = cO, si 
Fermi points, we have that 



(4.7) 



and therefore, {fi\Fp_h^auc''^ka\^^) — {fi\Fp-h,au\GS). Since the momentum val- 
ues Qi of Eq. ()4.3j) and the canonical momentum values qi of Eq. ()4.4j) are slightly 
different, the general matrix elements of the form {fi\Fp-h,au\GS) are finite, even 
for ^ 1. 

In conclusion, after taking into account the occupancy configuration trans- 
formations produced by the A and B processes, the typical matrix elements 
{fi\Fp-h,auciJGS) can be rewritten as {fi\Fp_h^au\GS). The states \GS) and \ fi) 
correspond to the densely packed occupancy configurations given in Eqs. fl4.3p 
and ()4.4|1 . respectively, and the hermitian conjugate of the operator Fp_h,ai/ is 
given in Eq. ()4.5|) . This overlap involves an excited state whose particle-hole 
occupancy configuration includes pseudofermion particle-hole pairs in the 
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au = cO, si bands. In the subsequent section, we will then sum over all pos- 
sible particle-hole occupancy configurations corresponding to the same tower of 
states. We recall that the value of the phase shift Q^ui'^lFau) is constant for each 
particle-hole tower of states. 

The only remainder of the finite energy pseudofermion (or pseudofermion 
hole), is then a Fourier canonical momentum summation, introduced by the 
Fourier transform of the corresponding pseudofermion operator defined on the 
effective au lattice. This Fourier transform is formally treated in Ref. |lH()j . 
The finite energy canonical momentum will in the following be denoted qau, for 
au = cO, si. In section ()4.2|) . this summation will be crucial to the derivation 
of closed form expressions for the full spectral function. We note that the en- 
ergy lAE and the momentum lAP are functions of these canonical momentum 
variables. 



4.1.2 Convolutions 

Having dealt with the contribution from the operators creating or annihilating 
pseudofermions away from any of the Fermi points, it still remains to treat the 
problem involving the matrix overlap with the operator Fp_f^^j^, Eq. ()4.6p . In 
the following, we define the subspace (tau as equivalent to the subspace of Ref. 
|131j called " J-CPHS-az/-(C)" occurring for example in Eq. (45) of that refer- 
ence. Thus, the states of this subspace are described by pseudofermion occupancy 
configurations generated by the C particle-hole processes, with fixed values of the 
deviation numbers {ANau} and {AA^^^^^} (for au = cO, si and l = ±). Let a 
typical element of this subspace describe a particle-hole occupancy configuration 
with energy Uau = ^Eau{C) and with momentum ka^ = APau{C), according to 
Eq. ()3.16|) . We define the total particle-hole energy and momentum as the au 
sum of the energy and momentum of the individual branches: Au = Uco + t^si 
and Aoj/v = kco + ksi, respectively. Note that the last definition also defines the 
particle-hole velocity v. 

In order to evaluate the relevant matrix elements, we have already seen in 
section 1)4.1.11) that we can describe a general excited energy eigenstate in terms of 
pseudofermion creation and annihilation operators in canonical momentum space, 
acting upon a vacuum state. As a consequence, the spectral weight associated 
with the amplitude -Mq^^ = \ {fi\Fp-h,ai^\GS)\ can be written as a convolution 
between the cO and the si branches. This involves the definition of an auxiliary 
function that reads 
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= ^ f E E /""^ r^°(^)^-/-° _ _ ^^^^ . 

V ^ / e-^,, e-^^ ^0 J-sgn{v)Au;/v,o 

■ 6{Auj -u'- uJsi)5{Auj/v -k'- ksi) M'q^^^ ■ M'q^^ 

where Auj and Auj/v are the particle- hole energy and momentum, respectively, 
The expression of the amplitude -Mq^^ in terms of pseudofermion creation and 
annihilation operators acting on the pseudofermion vacuum reads: 

QFai^ 

+Ag 

In Eq. fl4.8|) . the factor {1/Na) in the second line stems from the introduction 
of the discrete k' sum. The jacobian that arises when making the momentum a 
continuous variable cancels this factor and the pre-factor of (27r/L)^ arises due to 
turning the Kronecker (5-functions into Dirac 5-functions, according to Eq. ()3.17p . 

The introduction of sgn(f ) is necessary in order to keep the integration limits 
of the k' integration in proper order: for / = — we have that the particle-hole 
excitations grow in the negative uj direction whilst for / = + we have the opposite 
situation. In other words, in the tower of particle-hole states, positive velocities 
for / = + is equivalent to negative velocities for / = — and vice versa. This ensures 
that the lower integration limit is always smaller than the upper integration limit. 

The energy Aa; runs by definition over a small range of particle-hole energies 
from zero to a small negative number (RHB) or from zero to a small positive 
number (LHB), whilst the momentum runs over a small symmetrical interval 
around zero. 

The " Luttinger contribution" specified in section ()3.2.5p is not described by the 
function given by Eq. ()4.8p . In the region of small excitation energy, the branch 
hne group velocity equals one of the velocities of the particle-hole excitations. 
In this case, the pseudofermion or pseudofermion hole is created at the Fermi 
points by the B processes. We can then view the build-up of the a-branch line, 
as moving this pseudofermion along its dispersion starting at one of its Fermi 
points, and in this way tracing out the branch line in the (fc, uo) plane. However, 
for the first low lying energies that the av pseudofermion assumes after having 
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left the Fermi point, the dispersion is in its hnear regime. This means that for 
this intermediate region, the C processes replace the A processes in generating 
the branch line. This special case will be treated separately in section ()4.2.5p . In 
the following, we will always assume that we are outside of the region belonging 
to the Luttinger regime. 

We note that the final expression of Eq. ()4.8|1 is equivalent to Eq. (45) of Ref. 
|131j . by identifying the spectral function Bq^^ occurring in that reference with 
the expressions found here, according to: 



4.1.3 The energy cutoff Q, 

In this section, we will focus on the highest possible excitation energy value that a 
tower of particle-hole states can achieve. This energy is then ultimately a measure 
of how good an approximation the linearization regime will be, since, by forcing 
this value to be very small, we obtain a low tower of states albeit with good 
accuracy. Letting this value grow, we indeed obtain a higher tower of particle- 
hole states, allowing us to take into account more and more particle-hole processes 
and thereby aspire to account for almost the entire spectral weight of the problem 
(as measured by the sum rules). However, if we let this energy to be too large, 
we will start to consider non physical processes as the linear approximation of the 
dispersion relations become less and less accurate. 

This quantity hereby introduced will be called the energy cutoff, denoted by Vt. 
There is in principle an exact value for fi, which depends on the width of the linear 
regime of the pseudofermion band and on the four 2A^^'s, given in Eq. ()H.57j) . 
However, since the dependence of VL on these quantities is an open problem, we 
will replace the exact value by an effective value. In the following, we will assume 
a constant cutoff for the entire (/c, u) plane for all types of contributions. We can 
then deduce an average value of such an effective f2 by imposing the sum rules. 

As a prologue to the introduction of the cutoff into our expressions, we should 
first mention that the pseudofermion velocities 





(4.10) 
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obey the following relations (for m 0): 

U 

Vsl ^ VcO < n < 1, — < oo 

t 

Vsi^O ^ ^ oo (4.12) 

U 

Vsl ^VcO J ~^ 

which are easily deduced by suitable expansions of Eq. ()2.76|) . 

From Eq. ()4.8|1 . we see that the domain of Au is defined by the value of Q. 
Indeed, we have that < lAuj < Qhj the definition of the cutoff. For the energy 
u due to the A and C processes, we have that u = Au + IAE. In other words, the 
criteria induced by the cutoff is that the spectral function treated in Eq. ()4.8|) . 
should be multiphed by two ^-functions: 9 {l[uj — lAE]) 9 {Q — l[u — lAE]), in 
order for the energy Aa; of that equation to stay in the valid regime for the 
particle-hole tower of states. 

For the momentum part, we will use the fact that between the two pseud- 
ofermion branches available for particle- hole processes, it is the si branch that 
always has the lowest Fermi velocity between the two, as stated in Eq. fl4.12j) . In 
other words, the velocity v of Eq. ()4.8j) must be such that |f | > Vsi. This can be 
implemented by introducing 9{\v\ — Vsi) = 9{l [u — lAE] — Vsi\k — IAP\) 

Thus, we arrive to the following expression for the full spectral function: 



B\k,uj)= ^ ^ ^ 9{l[uj -IAE])9{VL-I[uj -lAE]) ■ 

■ 9(l[uj- lAE] - Vsi\k - IAP\^& ( fc ~ ZAP ' ~ ^^^) ^^"^^^ 
B'{v,Au) = ^-^ du' / dk' BQ^,{k\uj')BQ^,{Auj/v-k\Auj-u') 



where the primed summation in the first line stands for nai/=co si 
case of the 2P contribution, V- for the au = cO or au = si IP contribution 
or just 1 (no summation) in the case of the OP contribution. Here, the variable 
qav stands for the canonical momentum of the created or annihilated finite energy 
av = d), si pseudofermion. 

Note that when creating this finite energy av pseudofermion, the range of the 
qau canonical momentum summation goes from —(f^y to —qpau and from qpau to 
I'au- When annihilating the finite energy av pseudofermion, this range goes from 
—QPau to qpau- 
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4.1.4 The lowest peak weight ^(^,0) 



In order to calculate the lowest peak weight, from which the particle- hole excita- 



tions part, we need to consider matrix overlaps with F 



t 

p—h,au 



1. This means that 

we do not allow any particle-hole excitations for the lowest peak weight, defining 
A^^'^^ as the spectral weight associated with the A and B processes. To evaluate 
the matrix overlap we will use the following general result, which is easily proved 
by induction: 



M 



M 



Koin4n4io)i 



/4 ' flM 



4m' •^9Af 



(4.14) 



In other words, the evaluation of the spectral function Bq^^ {k', u') with -Fp_^ 
1 reduces to evaluate a (A^^^^ + AN^,^) x (A*'^^, + AN^,^) determinant. The anti- 
commutator is given by Eq. ()2.114|) . but that expression simplifies here since all 
of the creation operators refer to the ground state with Qaui^) = for all q in 
the entire Fermi sea. This means that Eq. ()2.114|1 becomes: 



N* 



^iaav{q' -q) /2 ^~iQav{q) /2 _ 



sm 



}av{q) 



sm 



aav[q-q') 



(4.15) 



which can be substituted into the expression for the determinant. In the evalu- 
ation, we will use the following exact result for the so called "Cauchy determi- 
nants" : 



det ( ^ ) = (-1)*^(^^-)/^ . (4.16) 

After some algebra, we find that by introducing the expression for the anti- 
commutator into the determinant, we can express each entry as the difference 
l/[cot(g'/2) — cot(g/2)], which allows for a direct comparison with the Cauchy 
formula. We thus arrive to: 
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^(0,0)^ -Q ^(0,0) (4^^7) 



au=cO,sl 

i-1 



4(0,0) ^ t TT 



sm 



aui i=l 



Qauiqi)\ j=l 



n sin^ (^) 



where one should not mistake canonical momenta and momenta with the same 
index i as being equal: it is just a way to enumerate the momentum values, hence 
qi 7^ Qi in general. This formula is exact, however, it poses some problems from a 
numerical point of view. This is due to the fact that in this expression, we have 
as many scattering phase shifts as pseudofermions. This makes the problem of 
finding the dependence of on the system size and the filling a very tedious 
problem. In other words, we are not able to deduce a closed form expression 
for the productorials. The same expression is used in Ref. [MI for the limit 
{U/t) — >■ 00, but for a constant scattering phase shift Q. This simplifies the 
problem tremendously, since we have a similar cancellation of the factors in the 
above expression, as in section ()3.2.2|) for the relative weights. In that section, 
the calculations depended on the fact that for each tower of states, the scattering 
phase shift is constant. Here, on the other hand, we cannot depend on such a 
result. However, inspired by the solution of the problem in the (U/t) 00 limit, 
there are some asymptotic behaviors that can be deduced. 

Since the scattering phase shifts away from the left and the right Fermi points 
do not contribute to the leading order terms of the dynamical quantities, we will 
use a trial approximation of letting all the scattering phase shifts in the above 
expression be equal to Qauit'QFau), for t = ±. This would allow us to evaluate 
the lowest peak weight for each t, leading to a weight that we can call A^auT We 
then form the product 

4r^-n^K (4-18) 

which should be a reasonable approximation in the thermodynamic limit. Note 
that this approximation must coincide with the expressions of the known limits, 
in this case the limit (U/t) 00. Inspired by the results of Ref. [ZSj-IHni, we 
here propose 

4(0,0) ^ V fiQaujl^ClFau) - Tt) 

where the even function /(x) is the same function as that of Ref. jHO], and Sau is 
an unknown quantity, depending on the density, the magnetization and {U/t). By 
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comparison with the corresponding expressions of Ref. [HHI, we have ScoSsi — * 1, 
and at the same time Sco — > sin irna, in the hmit (U/t) —>■ oo. Thus, we must have 
Ssi —>■ ^/ sinvma in that hmit (for m —>■ which is always our case). Note that 
(1) the square root of the function f{x) will guarantee the correct behavior of 
the full lowest weight at (U/t) oo and that (2) we will in the following assume 
the quantity to be equal to its value in this limit. One should note that the 
argument of the function f{x) is here shifted by vr due to the fact that in Ref. 
[HDj, the ground state phase shift is taken from the ID Heisenberg model for the 
spinous, and is equal to tt exactly, whilst in our case due to the normal ordered 
formulation of the problem we have by definition that Qai^{iqFau) = always 
for the ground state. Another example of this is also related to the scattering 
phase shift, namely the exponent in the denominator. We note that in our case, 
we can create or annihilate any number of pseudofermions at the Fermi points, 
whilst in Ref. jHOI, there is always exactly one spinless fermion and one spinon 
being created or annihilated, respectively. This is manifested by the absence of the 
numbers ANau,L in that reference, which ultimately leads to a different expression 
for the exponents. Thus, the generalization from Ref. [80: is two-folded: from 
infinite repulsion to arbitrary repulsion and from one quantum object creation / 
annihilation to many quantum object creation / annihilation. 

The form of Eq. ()4.19|) follows the form of the equivalent quantity of Ref. [101 • 
This form for the lowest peak weight is due to the equivalence of Eq. ()4.17j) to 
the corresponding expression of that reference. Due to this likeness, the thermo- 
dynamic limit expression of the lowest peak weight must by necessity be of the 
form presented in Eq. ()4.19j) . 

4.1.5 Merging the lowest weight and the relative weights 

Having dealt with the relative weights in sections ()3.2.2j) and ()3.2.3j) . and the 
lowest peak weight in the preceding section, we can finally put some real meaning 
to Eq. ()3.35|) . by noting that the weight ^2^^"'^.-'™'^'^.+) Qf ^i^g^^ equation is nothing 
but the spectral weight we would obtain using the derived expressions for the 
spectral function, with -Fp_;j given by Eq. ()4.5|) . This spectral weight is given by 
BQ^^{k', u') of Eq. ()4.10|) . which means that we can replace the explicit "absolute 
value of the squared overlap" with the following summation over the allowed 
positions in the {k, u) plane, of the particle-hole excitations 



(4.20) 




[mau,+ + m 
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which leads to 



where we have used the following identity for the ^-functions: 

Sia + x + y)5{h + x-y) = \5 - x\ 5 - y\ (4.22) 



and in the thermodynamic limit aau{rnaij,+, mau~) = aau,-{^au-)ciau,+ {^o 
Here and in the following we use the following expressions for the energy and 
momentum in terms of the numbers rriau/- 

= + > (4.23) 

From Eq. ()3.50|) we have that 

aauA^au,:) = 9{m^u,.) [T (2A^J]-^ (m,,,,)'^^""^ (4.24) 

where the ^-function seems unnecessary due to Eq. ()4.2Hj] . however it will be 
handy in the continuous limit as we shall see later. Substituting Eq. ()4.23|) into 
Eq. ()4.24|) . and substituting the resulting expression for aau,L into Eq. 1)4.211) . 
together with the expression of Eq. ()4.19|1 for the lowest peak weight, we arrive 
to the following expression: 



where we have defined the convenient sum 

2A,, = 2A+, + 2A;, (4.26) 
To recapitulate, we have now that 

B\v,Auj) = -^^ / du' / dk'- (4.27) 

■B^Qjk', uj')Bl^JAuj/v - k', Aco - co') 

Introduce now the following change of integration variables: 

x = uj'/Auj I sgn(f) , ,,,, (AuY , , , , 

' } =^ ^^dJdk' = -dxdy (4.28) 

y = sgn{v)vco k' / Aio j 27r ^ttVco 
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which means that < x < 1 and —l<y<l and also that 



lAu [x + tsgn{v)y] 



au = cO 



l[uj' + LVayk'] 



lAuj 



(r 






VcoJ _ 



(4.29) 



au = si 



Thus the quantity to integrate becomes: 



2tcv, 



cO 



n 

au=cO,sl 



Cl-2Ac,^ 



n 



r(2Ai„) 



L = ± 

X + Lsgn{v)y 



Y\_0\^x + Lsgn{v)yj 6\^1 — x + LVsiSgn{v) 

1 — X + ifsisgn(t>^ 
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2A'„-1 



(4.30) 



2A' -1 



where the first hne of this expression is constant in the integration. 
Let us define the function jF(l/ti) as 



Hl/v) = ^ fdx f dy n 



e[x + Lsgn{v)y^ e(l-x + Lv^isgn^v 



1 

If I 



y_ 



X + Lsgn{v)y 



1 — X + 6f sisgn(v) 



r(2A^i) 

1 y i\2A^i-i 

\v\ Vco\) 

This function plays the role of the function called Fq in Ref. |131j . and differs 
from this function only due to the various factors of fco and Vsi occurring in the 
definition of Fq. In this thesis report, we have chosen to collect all of these factors, 
as can be seen in Eq. fj4.32|) . We get, after some algebra simplifying the constant 
factors, that 



(4.31) 



BHv.Aoj) 



V 



2A,o „l-2A,i /IX \ Co-2 



(4.32) 



=± 



where we have defined Co = Yliau'^^oL'^- 

Note that this expression is equivalent with that of Eq. (61) in Ref. |131j . 
even though the definitions of the different quantities differ from each other. It 
is a simple algebraic task to extract all the different square roots and powers of 
VcQ and Vsi from Fq of Eq. (61) and (62) of that reference, to show that they are 
equal to the ones presented here. The difference in powers of these two velocities 
stem from choosing the cO branch to be the convoluting branch. This introduces 
a factor (l/fco) in the jacobian of Eq. ()4.28|) . 
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4.2 Closed form expressions 



4.2.1 The final step: canonical momentum integrations 

The expression for the spectral function derived in section ()4.1.5|1 . Eq. ()4.32|1 . 
is the base for all the expressions given in the remainder of this thesis report, 
except for the expressions for the Luttinger contribution, which will be presented 
separately. The function jF(l/t>), where v is the velocity of the particle-hole 
excitations, is crucial for the evaluation of the state dependent pre-factors of 
the spectral function. The particle-hole energy Aa;, and the exponent C^q — 2 oi 
Eq. ()4.32|) . will become crucial for the power law type behavior of the spectral 
function. The double integral occurring in the definition of T{l/v) can not be 
expressed in a closed form and has to be treated numerically. 

What remains now is to perform the finite energy pseudofermion summations, 
with summation variables gco and qsi in case of the 2P contribution, g^i in case of 
the s-branch, g^o in case of the c-branch, and with no summation at all in case of 
the OP Fermi point contribution. The canonical momentum qai, was introduced 
at the end of section 1)4.1.11) . and entered the derivation of the spectral function 
in Eq. ()4.13|) . Moreover, since we only consider one type of transition, called 
the "Basic" transition defined in sections ()3.2.(jj) and ()3.2.7|1 . we can drop the 
summations over the deviation numbers {ANay} and {AA^^j,^^}. 

We have that 



where a = c or s for the c-branch line and the s-branch line, respectively, and the 
two special cases tu) and B'-^^^^{k,uj) will be treated separately. 

The domain defining function T^ph{v, u) illustrates that we are always integrat- 
ing over such canonical momentum values, which all contribute to that certain 
point in the {k, u) plane, to which the A, B and the C processes bring the exci- 
tation to. Note that the integration domain, as controlled by T>ph, is very small 
as compared to the whole (/c, a;) plane. It covers a region in this plane, which 
comprises all particle-hole excitations that can reach the momentum k and the 
energy u. Note that the spectral weight at this point has contributions from 
particle-hole processes originating from many different surrounding points in the 
{k, u) plane. In other words, there are many A and B processes whose excitations 




Bop{k,u;) 
Vphiv,uj) 



Vph{v,oj) B\v,uj- lAE) 

e{i[uj - iAE])e{n - i[iu - iAE])e{\v\ - v^i) 



(4.33) 
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brings us to a finite momentum and a finite energy in the vicinity of {k,uj), and 
that by the C processes actually reaches the point {k,uj). This demonstrates the 
overlap of the different particle-hole tower of states, originating from the vicinity 
of {k,u!), and contributing to the spectral weight at [k^oj). 

4.2.2 2P contribution 

We will perform the 2P contribution integrations explicitly, to obtain the final 
expression for B\p{k^ oj). The calculations here are depending on the fact that the 
2P contribution covers a compact subspace of the (fc, uj) plane. Thus a typical 
point in this compact subspace is both reachable by a A and B process, and 
by a similar process for a slightly different set of canonical momenta values and 
then adding particle-hole excitations, due to process C. The limiting lines of this 
complete 2D subspace of the (fc, uj) plane are the lines which limit the effective 
Brillouin zone, and the so called border lines described in the subsequent section. 

Consider an energy uj and a momentum fc, reached by the A, B and the C 
processes. Due to the compactness of the subspace of the (A;, uo) plane, we can 
define a canonical momentum value for each av branch, such that 

u = l/\E{ql) 

k = lAPiC) (4.34) 

where lAE{q) and lAP{q) are given in Eq. ()3.60|) . Note the adopted shorthand 
notation: technically we should have written "lAE{q^Q,q^i)" (and similarily for 
lAP ) in the equation above. 

Let now g^^, + qau denote a typical canonical momentum value, that due to the 
A and B processes have energy and momentum lAE{q^j^+qau) and lAP{q^^+qau), 
respectively. We are looking to integrate over canonical momentum values (jau, 
such that the differences [u — lAE{q^^ + qau)\ and \k — lAP{q^^ + qciv)\ are 
inside the allowed region for the particle-hole towers of states, i.e. such that 
< /[cj - lAE] < Vt, and |f | > Vsi- 

This means that 

UJ - lAE{q^^ + qau) = -l^ sgn (AA^^^) Vay{qay)qav 

au 

k - lAP{ql + g,,) = -/ ^ sgn (A</) (4.35) 

av 

by expanding in powers of qav and retaining the linear terms only. Defining the 
energy variable uj' = uo — lAE and the reciprocal velocity variable z = 1/v = 
{k — lAP)/{uj — lAE), we can solve the relationships of Eq. ()4.35|) for q^u, to 
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obtain 



1 - Va'u'{qa'u')z 



Vau{qau) — Va'u'iQa'u') 
1 - Va„{qau)z 



Qa'v' — 



/sgn {AN%) uj' 



(4.36) 



At this stage, we remind ourselves that we are integrating the variable Qai,, 
as we are scanning for all possible A and B processes contributing to the energy 
and momentum given by u and k above. However, the domain of integration is 
governed by the quantity Vph defined in Eq. ()4.33|) . which is expressed in terms of 
uj' and z. The jacobian of this transformation is equal to uj' /\vau{(^au)~'^a'u'{(^a'u')\■ 
We thus obtain 



so that the final expression for the full spectral function in the case of the 2P 
contribution becomes: 



where we have u = ^(M^lsi), k = K^lcOjj^) as well as Co = Co(gco 

as defined in sections ()3.2.6|) (RHB) and ()3.2.7|) (LHB) (in these sections, the 

canonical momentum is denoted qav)- 

We see that this expression becomes singular as |fco(^co) ~ "^silo^ii)! ~^ 0- 
deed, the transformation from the q^y variables to oj' and z is not well defined 
as the two velocities approach each other, which is clearly seen from the jaco- 
bian of the transformation: it shoots off to infinity. The mathematical condition 
'^co(q'co) ~ traces out the border lines in the {k^u) plane. 



B2p{k,uj) 





(4.38) 
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4.2.3 Border lines 



The border lines are truly the "borders" of the spectral weight due to the A 
and the B processes. For example, consider the RHB "Basic" transition. The 
point in the (A;, a;) plane, reachable by the 2P contribution, corresponding to 
pseudofermion canonical momenta both equal to zero, belongs trivially to the 
border line, since i'co(O) = Vsi{0) = 0. But this point also has the smallest energy 
of all one electron removal excitations, i.e. eco(O) + 6^1(0) < eco(q'co) + ^si{(lsi) 
for all Qco 7^ and g^i 7^ 0. This means that there are no finite energy and 
finite momentum processes available to put spectral weight below this point. The 
same reasoning can be applied to the one electron addition case as well. Note 
however, that we can still reach areas outside the border lines by particle-hole 
excitations originating from a point sufficiently close to the border line itself. This 
effect forces us to consider two separate cases: one in which we consider areas 
of the {k, uj) plane available to either A and the B processes on the one hand 
or C processes on the other, and one in which the {k, uj) plane is reachable by 
the C processes only. We will thus have two different contributions, here dubbed 
^B^irderi^^^) former case and -Bso^derl^' ^) latter. 

Let us focus on -B^order(^' '^)- '^^^ parametric equations of the border lines 
are given by the energy cubl and the momentum ksL, according to 



au 

These are then the energy and the momenta of the border lines themselves. As 
before, we now fix a point in the vicinity of this line, with energy and momentum 
UJ and k, such that this point is reachable from a small but finite region of the 
border line by particle-hole processes. We introduce the small variation q^^, + q^v 
{av — cO, si), where q^^^ is a canonical momenta bringing the excitation to the 
fine. Since this fine is truly a "border" of the spectral weight, this means that in 
this first step, we are scanning for particle-hole processes slightly below {I = +) 
or slightly above {I — — ) the line itself, which motivates the introduction of qau- 
In these regions, we have points which are reachable by the A and B processes 
as well, whilst on the other side of these lines, the {k, uj) plane is completely 
void of any spectral weight from the A and B processes, so that the only spectral 
weight stems from particle-hole excitations originating from regions slightly below 
{I — +) or shghtly above {I — — ) the fine itself. 

We have now the following slight displacement in the energy and the momen- 
tum: 
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^E{q^^ + q^y) = Iujbl + VauiQau) sgn{AN^f)qau 

av 

AP(e + 9a.) = IkBL + sgn(AiV^/)g«, (4.40) 
The particle-hole velocity v is thus 

uj - lAE{q^^ + qc,y) uj - ujbl - lvau{qau) Y.au sgn(AiV^^)go 



lAP{qS, + q^,) I sgn(AiV„r )g„. 



which is equivalent to 

5:sgn(AA'r)«„. = -^^^^^ (4.42) 

where az/ = cO, si as usual and Vaui^au) could stand for either pseudofermion 
velocity, since they are equal to each other. Now, this equation permits us to ex- 
tract the relationship between the two different deviations (jau, that the border line 
demands. This constraint comes from the condition that the two pseudofermion 
velocities must be equal, and hence, a small deviation in the canonical momentum 
of one of the branches induces a small change in the canonical momentum of the 
other branch, in order to keep the velocities equal. We can simplify matters a lot, 
by restricting ourselves to infinite particle-hole velocity, which would make Eq. 
()4.42j) equal zero. Like this, we are only considering points in a straight vertical 
line originating from the point {kBLi^BL}- This means that we only consider such 
particle-hole excitations which are straight over (/ = — ) or straight under {I = +) 
the border line (which is something we can do without loss of generality since we 
will integrate over all relevant momentum and energy values). We obtain: 

^ sgn(AiV^^)g„^ = q^^ = -gau.a'u'qa'u' 

av 

9au,a'u' = sgn(AiV^/)sgn(AiV„^,^) = ±1 (4.43) 
This means, for example. 



AE{q^^ + g„,) = sgn(AiVj^^)e,o(g,°o + Qco) + sgn(AiVi{^)e,i(g,1 - gco,siqco) 



■lUBL+^-f 



sgn(AiVj^^)a,o(ge"o) + ^gn{AN^/)aM^, 



NF\ 



(4.44) 
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where the energy is expanded up to the second order in qco, and we have defined 
the pseudofermion accelerations aay{q) = dvau{,<i) / dq. Now, we need a way to 
treat the difference in the pseudofermion velocities. We obtain: 



= sgn(AiV„y sgn(AiVj,^)a,v'(?~.v') (4-45) 

a'u' 

and thus 



n — ^n-nM T^T N F aA^lau + (jau) Va'v'{(laiu' 9au,a'v'qau) ( A AP,\ 

g.. - sgn( AyV„, ) ^^^^^ sgn( AiVj,^)a„v' iCw) ^ ^ ^ 

which is the quantity that, when squared, can be introduced into Eq. ()4.44|) . to 
produce 



I {2l[uj - lobl]^ 

\Vav{qL + ^o.u) - Va'u' {Cu' ' 9av,a'u'qav) \ VEaJ«"^(^aJI 

(4.47) 

which finally yields the expression for the fully integrated spectral function, for 
the spectral weight in the vicinity of the border line, just above (/ = — ) or below 
(/ = +) this line: 



2e(Q-l[uoBL- ^])v 



-2Aco „,l-2Asi 
cO ^sl 



7rCo\/|aco(gco)l + k^ifeDI 
— j {21[uj-ujsl]) " j ^ dzT{z) (4.48) 



where a; — cj^j;, is a small energy, sufficiently small to be reached from the border 
line by some particle-hole processes, as demonstrated by the ^-function. Note 
that as before, the energy lo and the momentum A;, is connected to the energy 
and the momenta of the dispersive quantum objects through uj = uj{q^Q ,qsi)y 
k = HQco ' Qsi) as well as Co = Co(gco > Qsi) and ubl = ^BLiQco , Qsi)- Finally, the 
factor of 2 arises from the consideration that there are always two values {q^Q , 
(jgi) contributing to the spectral weight at the same point {k,uj). 

The other border line expression is very similar to the one derived above. 
The differences stem from the fact that above (/ = +) or below (/ = — ) the 
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border line, there is no spectral weight due to the A and B processes alone. The 
spectral weight of that region is generated by the C processes, on towers of states 
originating from the region below (/ = +) or above (/ = — ) the border lines. 
The first difference is due to the uj' integration of Eq. ()4.37|1 . In this case, this 
integration can not run from tu' = since we do not have any finite energy and 
finite momentum processes at the (A;, cj) point under consideration. In fact, we 
must move away a minimum distance oi lo — lobl from this point, in order to 
reach the "allowed" region for the A and the B processes. Thus, the integration 
can only start at this energy value. Furthermore Eq. ()4.42|) introduces a factor 
1/[1 — zvav{(i^^\^ since we can no longer scan all particle-hole energy values from 
to i7 at infinite particle-hole velocity, as above. Moreover, we are integrating 
over a two dimensional region that is tilted with a slope proportional to Vaui^au)- 
The velocity of the border line measures at what angle it cuts through the region 
available for the particle-hole excitations. Due to the inclination of the line, there 
is thus a non negligible region for which the particle-hole processes can not enter 
the domain of the A and B processes. However, the size of the region which is 
available depends also on the value of the particle-hole energy (and yet for some 
particle-hole velocities, we will never reach the region allowed for the A and B 
processes). In other words, the interplay between these quantities influences both 
the region for which we will have a finite spectral weight as well as the limits of 
integration in the variables depending on the particle-hole excitations, in our case 
in the variable z = 1/v. These considerations are further explained in Fig. 1)4.11) . 
We will skip the mathematical details of this analysis and merely present the 
result. In the expression below, the first ^-function refers to a border line cutting 
through the particle-hole region in such a way that the "base" of this region is 
completely inside the allowed region for the A and B processes. This is why 
the accompanied integration can run over the entire particle-hole velocity range, 
from —Vsi to Vsi- The second ^-function corresponds to a border line cutting 
through this "base" and the accompanied integration limits are thus modified to 
only integrate in the allowed region for the A and B processes. Note that the 
integrand is always the same, with the term that was zero in the ordinary 2P case 
now replaced with a term proportional to l[u — ujbl]^°- The quantities with index 
au remain unspecified since these quantities are equal for both pseudofermion 
branches. 
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Figure 4.1: Schematic figures of the geometrical considerations needed for 
the integrations in the z = 1/v variable, for spectral weight contribution to 
Bs'^orderik-oo), foT l = +. In the left fi gure, we see a border line with a "smaller" 
velocity \vau{(jau)\- The border line divides the depicted region into two parts: 
below the line, where the A and the B processes allocate spectral weight, and 
above the line which is only reachable via particle-hole excitations from energy 
and momentum points below the line. The line with dotted endpoints originates 
from a point under the line, and through the C processes reaches the point {k, lu). 
Note that the arrowheaded horizontal fine covers the entire particle-hole range, 
i.e. —1/vsi < z < 1/vsi, in contrast to the right figure where the "larger" value 
of |vai/(?Qi/)| causes a cut of the base of the triangle. On the right side of this 
line, we have only spectral weight due to the C processes. The integration range 
is fimited by the velocity of the border fine, but also by the value of a; — AE. 
Indeed, if we consider the left figure with a larger value of cu, the border line 
will cut through the base of the left triangle as well. Note that for the spectral 
function B^^^^{k.uj), the {k^uj) point is below the border line. 
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(4.49) 



4.2.4 The a-branch lines 

The a-branch lines [a = c, s) are lines which contribute significantly to the over- 
all shape of the full spectral function. These lines are formed by the au pseud- 
ofermion or pseudofermion hole assuming values along the entire range of its 
dispersion, whilst the other pseudofermion or pseudofermion hole is created or 
annihilated at one of its Fermi points. The only reservation we will have in this 
section is that when the dispersive pseudofermion is sufficiently close to either 
one of its Fermi points, the mathematical treatment of the problem will become 
different, and will be dealt with in the subsequent section. 

We saw in sections (I3.2.6|) and ()3.2.7p . that the scattering phase shifts have 
different expressions in the 2P compared to the branch line cases. Indeed, 

in the former case we have two scattering centers, which disperse in the many 
body system. Here, however, we have only one dispersive scattering center, as the 
other pseudofermion or pseudofermion hole is confined to one of its Fermi points. 
Thus, the displacement of the integrating pseudofermions in sections ()4.2.2|1 and 
()4.2.3|) . is in this section confined to a line. This means that we fix a point in the 
{k, u) plane in the vicinity of the branch line, and then integrate over canonical 
momentum values on the line only, in order to account for the spectral weight 
due to the processes that brings us from the line to the point in the {k, u) plane 
under consideration. In other words, if the integrating pseudofermion would leave 
the branch line, we would be considering spectral weights described by another 
set of the quantities 2A''^j^, than that of the branch line itself. 

We now have that the small canonical momentum denoted qau, only varies 
inside a small domain on the branch line such that the {k,uj) point can be reached 
by particle-hole processes with energy less than or equal to Q. Since (jau will vary 
on a tilted line, with the amount of inclination proportional to the velocity of 



137 



the branch hne, we have to consider a situation topologically equivalent to the 
one considered for the spectral function B^,^^^^{k, to), i.e. the situation where the 
slope of the branch line might be such that for some values of the particle-hole 
velocity v, a non negligible region of the tower of states can not be reached from 
the branch line, assuming that we only allow particle-hole energies between and 
fl. This topological effect is depicted in Fig. ()4.2|) and further discussed in Ref. 
|131j . The effective impact that these considerations have, is the introduction of 
the integration limits and the 6'-functions in the expression given below. 

Motivated by the discussion above, we now introduce a small canonical mo- 
mentum value Qau, which varies on the branch line around the canonical mo- 
mentum value g^j,. Thus the energy and momentum of this slightly displaced 
canonical momentum will be 



AE{q, 



NF\ 
av ) 



- sgn(AA^, 
:sgn(AiV,y 



Qau '^au \Qau ) 



(4.50) 



To integrate over the particle-hole contributions, we will need the particle-hole 
velocity v, which is readily found to be: 



u - lAE{q 



r,0 

av 



qav] 



NF' 



k - lAPiq^, + ~q^,) 
We thus obtain 



^a.(e)-sgn(AiV^; 



NF\ 
au ) 



l[uj — luJo\ 



where oja = sg'n.{AN^f)e{q^ 
The jacobian becomes: 



is the energy of the a-branch line, where a 



1[UJ — luJa\ 



-dv 



l[uj — luJa] 



JauyHau 



[I - zvau{q^y)y 



rdz 



(4.51) 



(4.52) 



c, s. 



(4.53) 



where the neglected sign is later taken care of when defining the integration limits. 
This jacobian is presented here because it will change the behavior of the spectral 
function in the energy. The particle-hole energy occurring in the argument of the 
spectral function & can now be expressed as l[uj — lua] ■ v/ {y — Vau{q^y))- By 
expressing the energy in this form, we show that we are scanning the branch line 
in such an interval, where it can reach the (/c, uj) point under consideration. Note 
for example that directly under this particle-hole point, i.e. for v = ±oo, we have 
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Figure 4.2: Schematic figure for a typical brancli line integration, for / = +. The 
line with dotted endpoints has one end fixed at {k, uj), and one end varying on the 
branch line itself. The latter point must not be further away in energy as what 
is dictated by Q. Note that, depending on the branch line velocity and the value 
of — AE, we will have a similar situation as already discussed for the border 
line case, accounted for in Fig. ()4.H) . This means that the branch line may cut 
through the base of the triangle in this picture, rendering a smaller integration 
interval of z than what is depicted here. 



that this energy expression is equal to /[cj — lua]- We are then, by introducing a 
small quantity Qau, and allowing it to be both positive and negative, scanning the 
branch line in an interval that covers both sides of this energy point. Thus, by 
defining the energy in this way, we are automatically accounting for the region of 
the branch line of interest. We have that 



& (v, 1[UJ - lUo] ■V/{V- Vo^uiQau)) 



4.['t^Y"f,. m_ (4.54) 



So that the full spectral function in the vicinity of the branch lines becomes: 
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(4.55) 



where the factor sgn(g^j,) is introduced together with the integration hmits, in 
order to always produce a positive number from the z integral. Note that this 
expression becomes singular whenever we approach the branch line for states 
such that Co ~ 1 < 0. This expression will be responsible for the characteristic 
line shapes of the spectral function, following the dispersion of the dispersive 
pseudofermion or pseudofermion hole. However, when it enters the linear region 
of its dispersion, this expression ceases to be valid. 



4.2.5 The Luttinger contribution 

The "Luttinger contribution" is a special case of the a-branch line, defined in 
section ()3.2.5|) . where the dispersive pseudofermion is very close to one of its 
Fermi points so that the dispersion relation is in its linear region. This case 
needs to be treated separately from the general a-branch case since the formulas 
applied in that case are not valid as the dispersive pseudofermion enters the linear 
region. We remind ourselves that in this region, the dispersive pseudofermion is 
in the same region as some of the particle-hole excitations. In other words, the 
created pseudofermion or pseudofermion hole, and the particle-hole excitations, 
share the same velocity. In this way, the "Luttinger contribution" case arises from 
a "velocity resonance effect". 

For this reason, we have to take a step back in our analysis, all the way to Eq. 
fl4.27|l . In this equation, we must change the domain of the k' integration, due to 
the limited range of momenta available for the linear regime. We will introduce 
a small quantity denoted Aq which measures the width of momentum over which 
we will integrate, for each value of uj'. This more careful procedure of integrating 
the particle-hole processes accounts for the linear regime, as the integration runs 
successively along the dispersive line, according to the integration limits of the k' 
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integral. We will thus let 

I duj' / dk' — > duj' / dk' (4.56) 

We then define Ag = 2Txy/L where ?/ is a number between and 1, and will be 
further specified later. Since the Fermi velocities differ from each other for the 
different branches and for the two different Fermi points of the same branch, we 
will change the notation of the spectral function B\ to In the following we 
let cO = si and si = cO: 



1 /-AlJ l-iU}' /Vav 

B''' = - dJ dk' B^Q^^ik', uj')B'q^{^uj/v - k', Au - u') 

^ Jo JiiLu'/Va^-Aq) 



Aq 



/ du' Bq^^{luo' /vav — LAq/2 , lo')Bq_{Alo /v — LLO' /Vau + LAq/2 , Auj — Lu') 
Jo 

This intermediate step can now be continued by using Eq. ()4.25|1 . and by 
changing integration variable from uj' to x = uj'/Auj. There will be some constant 
factors in the following expressions, which we for now bundle up into one overall 
constant, denoted by C. In this way, we will arrive to the following expression 
after some straightforward algebra: 

C^(/Acu)'^°-='(47rt;,,)2-2^-(47rt;^)2-2^^ [' dx ■ (4.57) 
271" Jo 
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iAq] 
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, [1 LX iAq]Y^^-W 

1- X + LVa!7< h —r- > ] [ X + L Vo 

[V Vau 2 J y V 

We note now that for the au branch, the argument of the ^-function is always 
larger than 0, which makes it convenient to perform the l' product explicitly for 
this branch, which leads to 

^ / j LX LAq\\ ( f LX .Ag^^'^-^-' 

[[0 \X + LVau < — \ j [X + LVau 



i'=± 



By extracting all the different exponents of the Fermi velocities, of the energy 
Au and of the momentum Aq, as well as of numerical factors and factors of vr. 
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we find that by defining the following function: 



TcvA^jv) = Do 

■ Y[e(i-x + 



a'u'=cO,slL'=± -I 



dx {2x) 



(4.59) 



i'=± 



LUVt, 



-X 



1 -X + 



X 



V, 



we can reach the following expression for the spectral function: 



(4.60) 



where we have defined (^au^t = Co ~ ^^au ^^e factor (1/Na) comes from the 
definition of the quantity y occurring in the definition of Aq: y = nIt^^"^^""' . This 
y is chosen so that y ^ Q when '2A''^^ and so that y ^ 1 when 2A^j^ 1, 
where it is assumed that for the Luttinger case < 2A^j^ < 1. 

The remaining procedure is now exactly equivalent to that of the a-branch 
case, with the same considerations as already dealt with. This is the consequence 
of the Luttinger contribution being a "special case" of the a-branch line: what 
remains is to disperse our pseudofermion, albeit confined to the linear regime 
of the dispersion relation, and fix a point in the (/c,a;) plane which is in the 
vicinity of this dispersive line. Where in the a-branch case we dealt with the 
energy difference cu — cua, here we deal with a similar energy difference u — ujLutt- 
The exponent of this energy difference in the former case was Co ~ 1) here this 
exponent is C,av,L — 1- The arguments of the ^-functions will consist of Fermi 
velocities, and not the general momentum dependent velocity |''''qi/(9)|, simulating 
the confinement to the vicinity of the Fermi points. Moreover, the jacobian of 
the integration will be identical to the a-branch case, but with a velocity equal 
to iVav instead of a velocity Vav{q^jj). Thus, with an analysis identical to the 
one of the a-branch line, we obtain the full spectral function for the Luttinger 
contribution: 
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(4.61) 



4.2.6 Fermi point contribution 

The last case to consider is the OP case, where both pseudofermions are confined 
to one of their l Fermi points. We will have 4 different points in the (fc, a;) plane 
of this kind. The spectral function expression will then be valid in the vicinity 
of these 4 points only, as the particle-hole contributions allow the spectral weight 
from this case to extend a maximum energy of VL from the Fermi level. There 
is actually not much work to be done to account for the spectral function for 
this case, since there is no finite energy pseudofermion to integrate over: both 
pseudofermions or pseudofermion holes are confined to one of their Fermi points. 
The resulting expression for the spectral function in this case can thus easily be 
read from Eq. ()4.32j) : 



u) = e{n - lu;)e{liu) ■ J^{l/v) (4.62) 

where the introduced 6'-functions restrict us to particle-hole energies between 
and Q. We note that this spectral function contributes in the vicinity of the 
specific zero energy momentum points specified in sections ()3.2.6|) (RHB) and 
(j3.2.7|l (LHB), respectively. However, of all these points it is only in the vicinity 
of {kp, 0), that the spectral function has a singular behavior. Indeed, for the other 
points, the exponent Co~2 is positive and hence the spectral weight vanishes as the 
zero energy level is approached. Moreover, in the vicinity of [kp, 0), this spectral 
function has the smallest value for the exponent of all exponents derived. 
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Chapter 5 

One Electron Spectral Weight 



5.1 General {U/t) and n dependence 

The expressions for the spectral weight distributions derived in chapter ^ depend 
on the value of the ratio of the effective Coloumb interaction strength U and the 
transfer integral t, as well as on the electronic density of the system n (note that 
we always have the magnetization m — >■ 0). This follows from the {U/t) and 
n dependence of most quantities involved, as for example the phase shifts, the 
dispersion relations eauil) and their corresponding group velocities VauiQ) (these 
quantities are defined by Eqs. ()2.6H) . ()2.7H) . ()2.73p and ()2.76p . respectively). The 
dispersion relations determine the shape of the branch lines in the {k, u) plane. 
Since one of the pseudofermions or pseudofermion holes is created at one of its 
Fermi points, we can associate these characteristic lines with distinctive charge 
type (cO) or spin type (si) excitations (where the "type" of the excitation stands 
for the charge or spin content of the dispersive pseudofermion or pseudofermion 
hole). These spectral features are only singular, however, for negative exponents 
which produce divergent expressions as we approach the branch lines. In the 
following, whenever referring to electrons, we will use units such that the lattice 
constant a = 1. 

We remind ourselves that all of the 2A''^^ quantities are larger than zero and 
thus Co > as well, where (o is defined in the text under Eq. ()4.32p . This 
means that the general 2P contribution does not exhibit any singular behavior, 
with the sole exception of the border line case. As we approach these lines, the 
exponent becomes negative, equal to (—1/2), and thus we would expect a singular 
"rim" along the line described by the condition fco(9co) — '"sii^siJ^ where k and 
u depend on g^Q and through the relationships given in sections ()3.2.6|) and 
()3.2.7|) . It is interesting to note that the "velocity resonance effect" of having 
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the two elementary excitations propagating at the same group velocity, produces 
significant spectral features far away from the zero energy Fermi level. 

From now on, we will change back to the original notation by letting q^u ~^ Qau 
denote the canonical momentum of the created or annihilated pseudofermion, 
independently if this is a finite energy pseudofermion or not. Coming back to 
the branch lines, these are controlled by an exponent = Co ~ 1) which may or 
may not be greater than zero (Cc denotes the c-branch line exponent and the 
s-branch line exponent). If we are aspiring to compare our theoretical results with 
experiments, then these singular features must be visible in a {k, u) photo emission 
or photo absorption scan of the spectral weight of the material in question. 

In the notation of sections ()3.2.(i|l and ()H.2.7|1 . we will find in the subsequent 
sections that the lines exhibiting singular behavior, correspond to the numbers 
Aco = f'co (RHB s-branch line), Aco = Lsi (RHB c-branch line), Lco = + (LHB 
s-branch line), and tgi = + (LHB c-branch line) respectively. For these cases, the 
momentum dependence of the branch line exponent is plotted for various values 
of (U/t). 

One can envision the a-branch line contribution as a dispersive au pseud- 
ofermion or pseudofermion hole, "moving" along the line dictated by its dispersion 
relation. As it disperses, the particle-hole towers of states gives rise to spectral 
features in the vicinity of the branch line. However, as this pseudofermion or 
pseudofermion hole reaches one of the end points of its dispersive line shape, it 
enters an intermediate regime where the valid expression for the spectral weight 
is not that of the a-branch line, but rather that of the Luttinger contribution. 
This special case can be likened to the border line case, since it too arises from a 
"velocity resonance effect": when the dispersive pseudofermion or pseudofermion 
hole enters the regime where its dispersion relation becomes linear, it has a veloc- 
ity equal to one of the velocities of the particle-hole excitations. In this case, the 
spectral features are described by another exponent than that of the a-branch 
line case. 

The exponents obtained for the Luttinger contribution is equivalent to the ex- 
ponent obtained by low elementary excitation energy methods, such as conformal 
field theory [T^ . 

As the above mentioned au pseudofermion or pseudofermion hole approaches 
the very end of its dispersive line shape, it enters yet another regime within the 
Luttinger-liquid behavior, here dubbed the OP regime. In this regime, notable 
singular spectral features can be found, as it corresponds to the most divergent 
exponent. We will expect some diverging peaks at the zero energy Fermi points, 
i.e. for k = kp. However, for other integer multiples of kp (at the zero energy 
level), the exponents are positive and does not give rise to any singular features. 
Lastly, since kp = 7m/2, we have that the distinctive Fermi point peaks, as well 
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as the shape of the branch hnes, change proportionally to n. As we approach half 
filling, n — > 1, more and more spectral weight is transferred from the LHB to the 
UHB. Indeed, for n = 1, the LHB is completely empty. We will briefly discuss 
the filling dependence of the exponents in section ()5.4|1 . 

Thus, the two ingredients one needs in order to deduce the general spectral 
weight behavior are the dispersion relations and the values of the exponents, 
respectively. However, the greatest numerical challenge is to compute the pre- 
factors of the spectral function expressions, which were found to be proportional 
to an integration in the z = 1/v variable of the function T{z). The results for 
the theoretical spectral features reported in the remainder of this thesis work 
were obtained by employing the programming language Fortran as well as the 
Mathematica software. 

Numerical considerations: 

Since we assume that the magnetization is vanishing, the si pseudofermion 
band is almost completely filled. However, at strict zero magnetization, we have 
that 2A^]^ = by definition, since the si pseudofermions become non dynamical. 
This leads to an ill defined expression for the function J^{z), Eq. 1)4.311) . Also, 
many quantities show discontinuities in this limit, i.e. that limm^o/(^) 7^ /(O)- 
See for example Figs. ()2.5j) - ()2.7|1 . where g') has a sudden jump at the 

boundary q = ±fc^| = i/c^. Therefore, we do allow a very small yet finite magne- 
tization in our calculations. In this way, we avoid problems with having coinciding 
values for the Fermi momenta and the limiting momenta for the Brillouin zone. 
Typically in our calculations, we have kp-^ — kp^ < 0.001. 

To arrive to a suitable value for our cutoff Q, we assume for the following 
discussion that we are in the vicinity of a spectral feature described by a negative 
exponent. The particle-hole tower of states will then produce a decaying tail, as 
the particle-hole energy increases. The value of the cutoff Q must be chosen in 
such a way as to properly account for this decaying tail. If chosen too small, there 
will be an unphysical (abrupt) end to the tower of states, producing a step-like 
feature at the cutoff. However, if chosen too large, we will take into account 
unphysical processes as we approach the cutoff energy. Moreover, the cutoff has 
to be chosen so that the sum rules are satisfied. Under the approximation that the 
cutoff is only weakly state dependent, we find that an average value of f2 ^ 0.2t 
produces a spectral weight that fulfills these criteria. 
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5.2 The RHB spectral weight 



In this section we present the full one electron spectral functions for the RHB, 
as obtained by use of Eqs. (IOH|l . (lOHll . (gSH), and (jO^ . re- 

spectively. Furthermore, we present the momentum and {U/t) dependence of 
the a-branch line exponents, given by Co ~ 1 = ~1 + Yliav'^^oLu^ where 2^av = 
^^^_i_2A^^ is defined by Eq. (I4.26|l . and 2A^j, is given by the expressions found 
in section ()H.2.fi|) . We also plot the regions in the {k,uo) plane where the contri- 
butions to the one electron spectral function generates a finite spectral weight. 
These regions were obtained by the defining equations for k and u;, respectively, 
presented in section (K12.(ij) . 

From section ()3.2.6|) . we see that the sign of the shake-up phase shift, equal to 
Aco = ±, can be combined with the sign of the Fermi point of that pseudofermion 
hole which is confined to such a point, when listing all possible branch lines. We 
remember that the spectral function is an even function of its momentum variable 
and that hence for simplicity we are only interested in positive momentum values. 
We now have two distinct c-branch lines, both involving a si pseudofermion hole 
being created at its positive Fermi point. 

In Figs. ()5.1|) . ()5.2|) and ()5.3|) . the lines originating from excitations with mo- 
mentum values outside the first Brillouin zone (i.e. such that k > n) are folded 
back into this zone. Moreover, for < < fcj? we have two c-branch lines joined 
at A; = 0. The one lowest in |ci;| can equivalently be described in terms of the 
other one, but with negative momentum values. This c-branch line segment is 
then folded over into the positive momentum region. The momentum dependent 
branch line exponents are plotted in ()5.3p . for the branches with negative expo- 
nents, i.e. for the singular c-branch line between —kp and Sfc^ and for the s-branch 
line between —kp and kp. All other branch lines have positive exponents and thus 
their weight decreases as we approach the branch lines. For the Fermi point con- 
tributions, the only negative exponent occurs at = kp- For {U/t) = 100 this 
exponent is (p ^ -0.867 whilst for (U/t) = 4.9 it is (p ^ -0.951. 

The spectral weight distribution for the entire {k, u) plane for arbitrary values 
of (U/t), n and magnetization m ^ 0, was obtained by the use of the expressions 
derived in the previous section. Our results for (U/t) = 100 should be very 
similar to other results valid in the large (U/t) limit. In Fig. ()5.4|) . we plot the 
spectral function for (JJ/t) equal to 100 and in Fig. for (JJ/t) equal to 4.9, 

respectively. The former case should then be compared with the corresponding 
Fig. ()5.6|) . originally presented in Ref. [7^] and valid for {U/t) oo only. That 
reference uses properties of the Hubbard model unique for the infinite repulsion 
case, and does not use the pseudofermion representation per se. However, the 
representation used in that reference is related to ours since in the {U/t) oo 
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Figure 5.1: The region of the (fc, uj) plane with a finite spectral weight from the 2P 
contribution (left) and the branch lines and the border line (right), respectively, 
for the one electron removal band (RHB) with {U /t) = 100, n = 0.59 and m 
0. In the right figure, there are two c-branch lines emerging from the point 
{k,u!) = {kF,0). The one extending towards smaller momenta is characterized by 
Aco = '-si = ~1) the other one by Aco = —^si = —1- The former line obeys 
k = Qco + kp and the latter k = q^o + Skp, in the m — limit. Both of these 
lines are folded back into the positive momentum section of the first Brillouin 
zone. Note the almost completely fiat si dispersion. The border line can be 
seen connecting the minimum energy points of the two c-branch lines, having 
VcoiQco) = Vsiiqsi) ~ 0. 
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Figure 5.2: The region of the (fc, a;) plane with a finite spectral weight from the 2P 
contribution (left) and the branch lines and the border line (right), respectively, 
for the one electron removal band (RHB) with {U /t) = 4.9, n = 0.59 and m — > 0. 
The main difference as compared with the (U/t) = 100 case, is that the si s- 
branch lines have now a non negligible energy width. The s-branch line between 
—kp and kp is characterized by Aco = '•co = il (the two choices of the sign 
produces two superposing line shapes), and the s-branch line between Skp and 
5/ci? is characterized by Aco = —Lco = — 1 {^kp is folded back to 2n — 5kp). The 
border line velocity vbl assumes all values in the domain vbl = VcoiQco) = Vsi{qsi) 
and brings the region of finite spectral weight to smaller energies than that of the 
(U/t) = 100 case due to the larger value of Vsi- 
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Figure 5.3: The value of the exponents for the c-branch hues (left) and the s- 
branch lines (right), for various values of (U/t), n = 0.59 and m — 0. For the 
s-branch line exponent, the tick mark at /ci?/2 is inserted to aid the eye. Note 
that the c-branch line segment between —kp and is folded over into the positive 
momentum region. The value of the c-branch exponent for {U/t) = 100 is almost 
constant. For the folded momenta values, the values of the c-branch line exponent 
for intermediate values of (U/t), is smaller than the corresponding exponent for 
(U/t) = 100. For this subbranch, the exponent is smaller for smaller (U/t), 
however for all other values of q, we have the opposite situation. The s-branch line 
exponent is always smaller for smaller (U/t), however with a decreasing difference 
as we approach the Fermi momentum kp- 
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Figure 5.4: The one-electron removal (RHB) full spectral function, for {U/t) = 
100, n = 0.59 and m ^ 0, as viewed from two different angles. The most 
divergent peak is to be found at the zero energy k = kp point. From this point, 
one c-branch and one s-branch emerge. For both of these branches, we have the 
intermediate "Luttinger contribution" which brings the spectral weight down as 
compared to the value at {k,uj) = {kF,0). For the c-branch, the total weight 
does not change significantly for negative velocities. However, as this line passes 
the zero velocity point the weight starts to vanish as we approach k = 3k p. The 
border line contribution can be seen to produce very little weight, however visible 
in the figure. The s-branch weight decreases continuously as we depart from 
k = kp and approach = 0. The entire s-branch is concentrated at excitation 
energies close to zero, due to the very weakly dispersing si pseudofermion holes 
for large values of {U/t). 
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Figure 5.5: The one-electron removal (RHB) full spectral function, for {U/t) = 
4.9, n = 0.59 and m — 0, as viewed from two different angles. Some of the 
features are similar to the ones of the (U/t) = 100 case, for example the strong 
divergence of the spectral weight at the point {k,uj) = {kF,0). The s-branch is 
however much more dispersive, as can be seen on the s-branch line feature between 
k = and k = kp. Note that the border line is no longer flat, as compared to the 
{U/t) = 100 case, mainly due to the larger si group velocity. 



153 




Figure 5.6: The one-electron removal (RHB) full spectral function from Ref. [7^ . 
for (f//t) = oo at quarter filling n = 0.5 and m ^ 0. Note the overall agreement 
with Fig. ()5.4j) . for example the c- and the s-branch line features, and the Fermi 
peak singularity. Moreover, the c-branch line feature fades away as it approaches 



limit, the cO pseudofermion becomes the "spinless fermion" of that reference. The 
spin part of that reference is described by the ID Heisenberg spin hamiltonian. 
Here we use it as a reference for the validity of our results. 

The features of the "large (U/t)" spectral function obtained by using the 
pseudofermion representation are described in the caption of Fig. ()5.4|) . It is 
a verification that the "Basic" transition is a good approximation to the total 
spectral weight of the large (U/t) spectral function, by comparison with Fig. ()5.(jj] . 
Considering other transitions (for example, the previously described "Exotic" 
transition) will only modify the total spectral weight very slightly. For example, 
these other transitions will not bring about new features to the overall spectral 
function, but add small corrections to the already existing features and ultimately 
make so that the exact sum rule will be satisfied. This is the reason of not having 
any tick marks on the z-axis of the figures presented here: considering more 
transitions could, however slightly, shift the total weight. The general shape is 
however directly proportional to the probability of finding the created electronic 
hole at momentum k and energy u. As a final remark, we note that our result 
that the 2P "background" contribution is indeed very small, is confirmed by the 
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studies of Refs. |ZH]-|Hn!. 

The main difference between the full spectral function of the large {U /t) case 
and of the intermediate {U/t) case, is the increase of the velocity Vsi- This 
increase introduces spin related excitation energies significantly departed from 
the zero energy level. Since we always have that Vcq > f^i, the border line covers 
the entire si band, but only a segment (symmetrical around zero) of the cO band. 
Due to the larger value of Vsi for the intermediate [U/t) compared to its 

value for large {U/t), the border line extends to higher energy excitation values, 
than for large (U/t) case. Our results for (U/t) = 4.9 are presented in Fig. 

5.3 The LHB spectral weight 

In this section we present the full one electron spectral functions for the LHB, 
as obtained by use of Eqs. (g^Hl), (jO^ - (gH}, and (gSSl), re- 

spectively. Furthermore, we present the momentum and {U / 1) dependence of 
the a-branch line exponents, given hj (q — 1 = —1 + J^au'^^^u, where 2Aau = 
T.,=±'^Ku is defined by Eq. (jOHl), and 2A^^ is defined in section (jTTTjl . We 
also plot the regions in the {k, u) plane where the contributions to the one electron 
spectral function generates a finite spectral weight. These regions were obtained 
by the defining equations for k and u, respectively, presented in section ()H.2.7j] . 
With an eye to the applications of the theory, presented in chapter (0), we have 
chosen the intermediate value of (U/t) to be equal to 5.61 and not 4.9 as in the 
RHB case. 

The one electron addition spectral function is described by creation of one cO 
pseudofermion and the appearance of one extra si pseudofermion hole. The LHB 
"Basic" transition is described in section ()3.2.7p . from which we find that there 
is a smaller number of branch lines than in the RHB case. 

This results in part from the subtle effect that the si pseudofermion hole is not 
created at the expense of a si pseudofermion. The number A^^i remains constant 
under this transition. The contrast to the RHB case described in the previous 
section can be described by the fact that the si pseudofermion current is zero for 
all values of the si canonical hole momenta Qsi different from qpsi = ^kpi ~ kp- 
For positive values of the momentum k, we hence have one c-branch and one 
s-branch in total. 

The general domains of finite LHB spectral weight is given in Fig. ()5.7p 
for {U/t) = 100 and in Fig. ((^3) for (U/t) = 5.61. Note that the cO total 
bandwidth eco(±vr)— eco(O) = 4t is independent of {U/t). Since a cO pseudofermion 
is created in the "Basic" LHB transition, having canonical momentum values 
between — vr and —2kp and between 2kp and vr, respectively, we see that the 
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Figure 5.7: The region of the (fc, (-j) plane with a finite spectral weight from the 2P 
contribution (left) and the branch lines and the border line (right), respectively, 
for the one electron addition band (LHB) with {U /t) = 100, n = 0.59 and m — 0. 
The si band is nearly dispersionless, in analogy with the large {U/t) case for the 
RHB. There is no finite spectral weight for < A; < fci?, as both branch lines 
originate at k = kp and extend into regions with larger k. The border line is 
nearly fiat, due to having Vgi ^ 0. 
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Figure 5.8: The region of the {k, uj) plane with a finite spectral weight from the 2P 
contribution (left) and the branch lines and the border line (right), respectively, 
for the one electron addition band (LHB) with (U/t) = 5.61, n = 0.59 and m — 0. 
The basic topology is the same as for the (U/t) = 100 case, but with a larger 
value for v^i, which infiuences the extension of the s-branch line and the border 
line onto higher values of u, as compared to the (U/t) = 100 case. 
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maximum excitation energy of the cO creation is governed by eco(±7i"). This value 
can be read off of Fig ()2.2j) . and increases for decreasing (U/t). This explains the 
difference in the scale of the energy axis in Figs. ()5.7p and ()5.8|) . The s-branch is 
described by momentum values k = 2kF — Qsi and runs through values between 
kp and 3k p, respectively, and not between and kp as in the RHB case described 
in the previous section. 

The values of the branch line exponents are plotted in Fig. ()5.9p . for various 
values of {U/t). Similarily to that case, we are always restricting ourselves to 
momentum values k such that < /c < vr, which explains the apparent double 
valuedness of the c-branch exponent. The two c-branch lines, corresponding to 
the two possible values of Lsi = ±, maps onto one single continuous line by only 
allowing positive values of k inside the Brillouin zone, in the following way: The 
c-branch line feature with Lsi = + ranges from k = — vr — kp to k = —Skp (for 
negative cO pseudofermion canonical momentum values) and from k = kp to k = 
IT — kp (for positive cO pseudofermion canonical momentum values), respectively. 
Similarily, the c-branch line feature with Lsi = — ranges from k = — vr + kp to 
k = —kp (for negative cO pseudofermion canonical momentum values) and from 
k = 3kp to k = n + kp (for positive cO pseudofermion canonical momentum 
values), respectively. 

For example, in the (U/t) = 5.61 curve of Fig. ()5.9j) . we see the value of the 
c-branch line exponent for the Lgi = + subbranch, for k ranging between kp and 
TT — kp. The other subbranch (i^i = — ), k ranging between 3kp and n + kp, starts 
at 3kp with a value of the exponent roughly around —0.23 and then decreases as 
k increases beyond k = it. The value of the exponent at k = tt + kp (backfolded 
to momentum n — kp) is roughly equal to the value of the exponent belonging to 
the subbranch with Lsi = + at the same momentum value. 

We note that the s-branch exponents are monotonously increasing with k, 
becoming larger than zero for a large segment of the total branch line. We are 
thus expecting that the spectral weight of the s-branch line will vanish completely 
as we travel along the line from k = kp towards k = 3k p. This characteristic 
behavior of the s-branch line is indeed verified in Figs. ()5.10|) and ()5.1H) . where 
the fading away of this line feature is evident. 

Another characteristic feature of the LHB is the importance of the border lines, 
which in general carry much more spectral weight than their RHB counterparts. 
This effect has been attributed to a van-Hove singularity [ZH]-|Si]j in the strong 
coupling limit. 

The exponents of the Fermi point singularities are (p ~ —0.889 for (U/t) = 
100 and (p ~ —0.965 for (U/t) = 5.61. The strong negative exponent in both 
cases motivates the high peak at the point {k, uj) = {kp, 0) for both values of {U/t). 
For the c-branch line, however, we have that Cc grows with decreasing {U/t), 



157 



<c(q) ^,(q) 




Figure 5.9: The value of the exponents for the c-branch hne (left) and the s- 
branch line (right) for the LHB, for various values of (U/t), n = 0.59 and m — * 0. 
The value of the c-branch exponent for (U/t) = 100 is almost constant. Note that 
the exponent for the c-branch line segment for momentum values g > vr is folded 
into the first Brillouin zone (the exponent with canonical momentum q = ir + kp 
is folded back into the Brillouin zone at the small break of the continuous line 
given by Cc(q'), visible between 2kF and 3kp). The c-branch line exponents start 
from q = kp with roughly the same slope, but with an increasing value of the 
curvature for decreasing values of (U/t), and thus the values of Cc(q') increases for 
decreasing {U/t). The s-branch line exponent is negative for momentum values 
close to kp, but grows monotonously with q and becomes positive after some 
specific value of q and remains positive for the remainder of the branch line. We 
would thus expect the s-branch line feature of the full spectral function to vanish 
as this value of q is reached and passed. Note that Cs('?) depends almost hnearly 
on q for {U/t) = 10, and for smaller values of this ratio the dependence is mainly 
concave, whilst for larger values of this ratio, it is mainly convex. 
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Figure 5.10: The one-electron addition (LHB) full spectral function, for {U/t) = 
100, n = 0.59 and m — 0, as viewed from two different angles. The strongest 
divergency of the spectral function occurs at the zero energy Fermi point k = 
kp. The c-branch originating from this point has a smoothly decreasing spectral 
weight and is almost vanishing at the other end (towards the zero energy point 
at = "ikp). The origin of the large rim of the border line is discussed in the 
text of this section. However, this line seems here a little bit rugged which is due 
to the numerically sensitive calculations as well as to limitations in our " constant 
cutoff' and LHB "Basic" transition approximations. 

producing a weaker spectral weight as compared to larger values of this ratio. 

The spectral function of Fig. ()5.10|1 should be compared with that of Ref. [TH] . 
given in Fig. ()5.12|) . In this reference, there is no division of the different types of 
contributions leading to the total spectral function, but all types of final states fall 
into the same mathematical treatment, in contrast to the pseudofermion method. 
This means that there is no clear division between the contributions of the branch 
lines and the border lines, for example. However, there is also a positive effect of 
this: there are no troublesome crossover regions in which it is not clear exactly 
which type of contribution should be valid. For the c-branch line, in the vicinity 
of the zero velocity point, we have that this contribution superposes on the border 
line contribution. Moreover, due to the flatness of the si band, it is not clear 
how important the Luttinger contribution will be for the border line, since it can 
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Figure 5.11: The one-electron addition (LHB) full spectral function, for {U /t) = 
5.61, n = 0.59 and m — 0, as viewed from two different angles. The non zero 
dispersion of the si band introduces a curved border line, as well as a curved 
s-branch line. The latter is vanishing due to the positiveness of the exponent (s, 
as confirmed by Fig. ()5.9|) . The exponent for the Fermi point contribution is 
even smaller here than in the {U/t) = 100 case, explaining the strong peak at the 
point {k,Lj) = {kF,0). 
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Figure 5.12: The one-electron addition (LHB) full spectral function from Ref. 
|79j . for (f//t) = oo at quarter filling n = 0.5 and m 0. This figure should be 
compared to Fig. ()5.10|) . We see that the characteristic spectral features shown 
here is also present in Fig. ()5.10|) . For example, the pronounced border line 
singularity, the cO pseudofermion branch line, si pseudpfermion branch line, and 
the strong divergency at the point {k,u!) = (kpyO), are all features accounted for 
in both figures. 



161 



be argued that the Luttinger hquid region increases as the dispersion relations 
becomes flatter, and thus "more linear". These small uncertainties cause the 
rugged appearance of the border line of the pseudofermion method, which is a 
numerical feature totally absent in the method of Ref . [ZH] • The reasoning behind 
this ruggedness in the iU/t) = 100 case implies that for lower values of (U/t), 
this numerical effect would be less pronounced, as the si dispersion becomes less 
flat and hence the Luttinger contribution more confined. This is verified in Fig. 
(15.1111 . where virtually all of this ruggedness is non existent. 

5.4 Density dependent exponents and disper- 
sions 

In this section we will briefly discuss the filling dependence of the a-branch line 
exponents, without plotting the full spectral function. As before, these exponents 
were calculated by use of the defining equations in sections ()3.2.fi|) and ()H.2.7j) . 
respectively. When discussing these exponents, we will exclusively focus on the 
same branch lines as in the previous sections of this chapter. Also, the energy 
dispersions presented here were obtained from Eq. ()2.7H1 . 

All of the figures presented here are valid for (U/t) = 10 and m — 0. In Figs. 
fl5.1H|l - ()5.15|l . we present the cO, si, and cl pseudofermion dispersion relations for 
n = 0.35, 0.59, and 0.85 respectively. 

For the RHB, the dispersing quantum objects are cO and si pseudofermion 
holes. Thus, as n decreases, so does the value of the Fermi momentum, and 
hence the cO and si pseudofermion holes will have a smaller canonical momentum 
range in which to disperse. This means that the characteristic peak at the Fermi 
point {k,uj) = {kF,0) approaches the zero momentum and zero energy corner in 
the {k,u!) plane, as n — > 0. The minimum value of the au = cO, si dispersion 
relations, namely €au{0), approaches zero in this limit. Thus, the a-branch lines 
will also shrink towards {k,uj) = {0,0). 

For n ^ 1, the situation is reversed: the cO and si pseudofermion holes will 
have an increasingly larger canonical momentum range in which to disperse. As 
before, the position of the peak at (/c, a;) = (fci?, 0) will move as n varies. The a- 
branch lines will extend over increasingly larger portions of the {k, u) plane, as 
the domains of the pseudofermion energies and canonical momenta increases with 
increasing n. 

For the LHB, the situation is somewhat different since in this case the dis- 
persive quantum objects correspond to the cO pseudofermion (and not the pseud- 
ofermion hole), and the si pseudofermion hole, respectively. For the cO pseud- 
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Figure 5.13: The cO pseudofermion dispersion relations for various values of the 
electronic density n, {U /t) = 10 and m 0. According to the definition of the 
Fermi momentum, qpco = 2A;f = Trn, we see that the Fermi momenta approaches 
the limiting values for the effective Brillouin zone as 1, and goes to zero as 
n goes to zero. The bandwidth is constant and equals 4t independently of the 
electron density. 
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Figure 5.14: The si pseudofermion dispersion relations for various values of the 
electronic density n, {U/t) = 10 and m ^ 0. The band shrinks with decreasing n, 
with a decreasing bandwidth, and with Fermi points ^qpsi = ^kp = ±7™/2 
as n ^ 
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Figure 5.15: The cl pseudofermion dispersion relations for various values of the 
electronic density n, {U/t) = 10 and m — * 0. Even though we do not consider 
final states with finite occupancies of cl pseudofermions, the dispersion is included 
here for completeness. Both the energy bandwidth and the values of the canonical 
momenta at the effective Brillouin zone boundaries decrease as n — ^ 1. This is 
basically an effect of the diminishing number of doubly occupied, and empty, 
rotated electron sites, respectively, in this limit. 
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ofermion, we have that for decreasing values of n, the domain in which the cO 
pseudofermion can disperse increases. For very low densities, the added cO pseud- 
ofermion will have almost the entire cO band "for itself, and will thus yield a 
more extended c-branch line as for higher densities. We remember that the cO 
pseudofermion domain of dispersion in this case corresponds to canonical momen- 
tum values in the (— vr, — vrn) and the (vrn, vr) domains, respectively. The s-branch 
line, however, follows the same general behavior as in the RHB case, with its 
branch line feature shrinking towards the point (0,0). 

Once again, the situation is reversed in the opposite limit, n ^ 1. Here, the 
added cO pseudofermion will only be able to disperse along the "wings" near the 
canonical momentum values q = ivr (corresponding to the shrinking regions for 
which eco{q) > 0). According to the sum rule of Eq. ()3.6|) . the total spectral 
weight of the LHB vanishes in this limit, transferring its weight to the UHB. 

Note that in all of these cases, the region in the [k, lj) plane with contributions 
from the border line, is governed by the value of f^i- Irrespectively if we study 
the RHB or the LHB, all canonical momentum values in the si band such that 
^si{(l) < have a corresponding canonical momentum value in the cO band, such 
that the pseudofermion group velocities are equal to each other. 

The a-branch line exponents are plotted in Figs. ()5.1()|1 and ()5.17|1 . Their 
dependence on the filling n is discussed in the corresponding captions. Generally, 
for the RHB, we have that as n decreases, so does the total spectral weight, in 
accordance with the sum rule. Indeed, the value of the negative c-branch line 
exponent increases with decreasing n. However, the s-branch line exponent does 
not. The conclusion of this is that a study of the branch line exponents alone, 
is not sufficient to characterize the behavior of the branch line spectral feature. 
As n increases we have that the RHB c-branch line feature exhibits a stronger 
divergent behavior as n ^ 1, even though this effect is not "dramatic". 

For the LHB, we have an interesting effect for the s-branch line. Namely, 
as the momentum increases from k = kp, the exponent for the s-branch line 
becomes positive at a certain momentum value larger than kp. However, this 
momentum value approaches kp a.s n decreases. Hence, the weight of this branch 
line must also decrease as n decreases. We do not detect any similar effect for the 
c-branch line exponent. It seems reasonable to assume that the disappearence of 
the spectral weight in the LHB case as n 1 is in the pseudofermion picture 
linked with the disappearence of a dynamical cO branch, rather than to the values 
of the branch line exponents themselves. 

In this chapter, we have studied the behavior of the a-branch line exponent 
only. For a complete understanding of the filling dependence of the total a-branch 
line spectral behavior, we need in addition to study the corresponding behavior 
of the pre-factors. 
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Figure 5.16: The exponents for the c-branch hne (left) and the s-branch hne 
(right) for the RHB, for various values of ra, {U/t) = 10 and m — > 0. For a 
fixed momentum value q = qo, we see that Cc(q'o) is increasing as n decreases 
(i.e. by following a vertical line upwards). Hence in general, the c-branch line 
feature diverges more slowly as the density decreases. Note however, that the 
value of this exponent at the Fermi point kp/ir = n/2, does not alter significantly 
between different densities. An interesting effect occurs at n = 0.85 for the c- 
branch line: the value of the exponent at the Fermi point is larger than the 
corresponding value at momenta between and kp, in contrast to the values of 
the other exponents of the c-branch line. For the s-branch line, no linear-type 
trend can be deduced, as the intermediate density value produces the smallest 
exponent for small momentum values. For any n, is always negative however 
increasing as q increases from towards kp. As with the c-branch line, the 
exponent for the s-branch line does not vary significantly at the corresponding 
Fermi points kp/'n' = n/2. 
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Figure 5.17: The exponents for the c-branch hne (left) and the s-branch hne 
(right) for the LHB, for various values of n, {U/t) = 10 and m — 0. For the c- 
branch line, the exponent is always negative and not varying significantly with n, 
roughly larger than —0.40 and smaller than —0.30. For the s-branch line however, 
as n decreases, the region for which Cs(Q') < shrinks, and thus we expect that 
the s-branch line becomes less and less significant. Note that the length of the 
s-branch line shrinks as n decreases, whilst we have the opposite dependence on 
the filling for the c-branch line. In conclusion, as n decreases, the weight in the 
vicinity of the s-branch line decreases whilst the weight in the vicinity of the 
c-branch line does not. 
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Chapter 6 



Applications - Experimental 
Spectral Weight 

6.1 The organic compound TTF-TCNQ 

We will in this section devote some attention to the organic material "Tetrathi- 
afulvalene Tetracyanoquinodimethane" , abbreviated TTF-TCNQ. For tempera- 
tures above the broken symmetry state (linked to a Peierls transition as further 
described in this section), it is characterized as a metallic "charge transfer salt" 
consisting of linear stacks of planar molecules. In the subsequent sections, some 
general properties of this material will be discussed conveying the reasons for why 
it constitutes a reasonable quasi ID material allowing it to be compared with the 
theoretical results obtained so far. It is not our intention to explain in greater 
detail the rich physical literature that exists regarding TTF-TCNQ, nor to make 
a "from first principles" derivation of its physical properties. This is outside the 
scope of this thesis report. However, the interested reader could use the references 
given in this and in the subsequent section for a deeper study of the properties 
of TTF-TCNQ. 

The charge transfer occurs between the two types of molecules, i.e. between 
the stacks: an approximate 0.59 electrons per molecule is transferred from the 
TTF to the TCNQ molecule, which drives the stacks metallic. In the metal- 
lic phase, the electrical conductivity is about three orders of magnitude larger 
in the intra-strack direction (for both molecules) than in other directions |137j . 
a property attributed to the crystal structure This is manifested 

through TT-type orbitals overlapping in the conduction direction, i.e. overlapping 
with neighboring molecules belonging to the same stack. The high conductivity. 
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as compared to inorganic metals, has been a key motivator for scientific inves- 
tigations since the 1970's. Actually, TTF-TCNQ is one of the most celebrated 
organic conductors which have been widely examined since its discovery |138j 
mainly regarding its electronic conductivity and optical properties, which have 
been thoroughly investigated in for example Refs. For our purposes, 

the material makes a good candidate for a "one dimensional electronic system" 
regarding the one electron removal spectral properties of the metallic phase. The 
crystal structure of TTF-TCNQ, along with its Brillouin zone, is depicted in Fig. 

The non Drude behavior of the metallic regime is described in Ref. |144j . a 
regime reached above the critical temperature Tp = 54 K. This critical temper- 
ature is about half the value of the predicted mean-field weak coupling value, 
a deviation attributed to strong ID fluctuations ^16j. Below this temperature a 
charge density wave builds up in the TCNQ chains, eventually turning the system 
into an insulator at 38 K (the corresponding temperature for the TTF stacks is 
49 K |144j ). The critical temperature is linked with a Peierls transition, where 
an electronic gap opens up due to the molecule lattice displacements |145j - [TT7] . 
This means that even though the theoretical model developed here refers to zero 
temperature, we have to study the material at temperatures higher than the crit- 
ical temperature in order to reach the metallic phase, for which our theory is 
valid. 

The charge transfer between the molecules shifts the intra-molecular density 
of electrons by utcnq = 0.59 for TCNQ and thus by uttf = 2 — 0.59 = 1.41 for 
TTF, i.e. with a hole density of 0.59. Thus, due to the particle-hole symmetry 
of our model, it suffices to study systems with a density of 0.59. In this way, 
"removing an electron" (RHB) for the TCNQ translates into "adding a hole" for 
the TTF (LHB) , and hence the spectral weight of both transitions can be mapped 
onto the same {k,Lj) region, with u <0. We must note however, that the transfer 
integrals for the individual stacks are different from each other. 

The fact that Coloumb interaction plays a key role for the electronic structure 
of TTF-TCNQ is not a new claim |136j |158j , and from physical properties other 
than the photo emission spectrum, one would expect that U ~ At, with a slightly 
higher value for this ratio for TTF, than for TCNQ. By comparing with the photo 
emission spectrum, it is a simple fitting procedure (with the pseudofermion energy 
dispersions as the fitting functions) to deduce the values of {U/t) yielding the 
best match between experimental and theoretical results JH] (note that the cO 
pseudofermion bandwidth is constant and equal to 4t). In the following, we shall 
use the results of Refs. |148j and [EI, where an almost perfect agreement between 
the theoretical and the experimental spectral weight is found for t = 0.40 eV and 
U = 1.96 eV for TCNQ, and t = 0.35 eV and U = 1.96 eV for TTF, respectively. 
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Figure 6.1: The crystal structure of TTF-TCNQ. The angles about the a axis, 
as measured with respect to the c axis, ensures maximum covalent bonding along 
the stack direction, hence making the electronic conductivity strongly anisotropic. 
As a result, the electrical conductivity in the b direction is three orders of mag- 
nitude larger than in the other directions. The Brillouin zone (right) shows the 
high symmetry points in the corresponding reciprocal space directions. Thus, we 
measure electronic momentum along the b* direction, identifying the centre of 
the Brillouin zone F, as the zero momentum point. 
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We note that these values of the Hubbard parameters produce an exact fit with the 
pseudofermion energy bandwidths. We thus obtain effective values {U/t) = 4.9 
for TCNQ and {U/t) = 5.61 for TTF, respectively, which furthermore coincide 
with the findings of Ref. |144j (and will be further demonstrated below). 

One last note should be made on the subject of "dimensional crossover", i.e. 
the phenomenon that the quasi ID material, for some reason and in some regime, 
changes its physical behavior to start acting like, for example, a 2D layer. In Ref. 
|15(ij it is confirmed that the intermolecular transfer integral is on the order of 5 
meV ~ ksTp in the TB direction. This means that at sufficiently low tempera- 
tures, we will start to see hopping between identical molecules belonging to dif- 
ferent stacks. Moreover, below the Peierls transition we will also have important 
contributions to the dynamics of the system from electron-phonon interactions. 
These considerations force us to only study the system at temperatures T > Tp. 

6.2 ARPES experiments on TTF-TCNQ 

ARPES is an abbreviation for "angle resolved photoelectron spectroscopy" and 
basically stands for the experimental procedure equivalent to what a theoretician 
would call " one electron removal" |149j . The basic physical considerations regard- 
ing the ARPES technique will not be accounted for here, other than just stating 
that the energy and direction of the photoemitted electron defines the quantum 
state of the material sample. Thus, varying these parameters, it is possible to 
obtain a full energy and momentum map. The ARPES technique does not de- 
mand a certain environment or sample temperature per se, and the experiments 
can be conducted at zero magnetic fields. Usually, this technique is employed for 
probing the shape of the Fermi surface of the sample, however in our case this 
technique is used to study properties of the bulk material itself. 

To motivate this, we use the reasoning of F. Zwick et al |15()j . In this reference, 
special concern is taken regarding surface effects: surface sensitiveness to radiation 
damage, aging of cleaved surfaces, eventual perturbations on the bulk charge 
balance and lattice periodicity due to surface effects, possibility of an insulating 
surface, and so forth. The "passing" of these tests for TTF-TCNQ allows us to 
conclude that the ARPES measurements actually do give us information on the 
bulk properties of the material. Some of the references that F. Zwick et al use 
in this analysis include Refs. |151j - [TFf] . with the main conclusion that cleaved 
surfaces of TTF-TCNQ are "highly ordered and retain the periodicity of the 
bulk" up to the penetration depth used by the ARPES technique. By keeping 
the material sample at an ambient temperature of 150 K, F. Zwick et al jl5Uj 
measure the spectral weight along the b well as along the perpendicular 
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a axis. The last measurement serves the purpose of excluding this direction 
as a conductive direction, showing no dispersive behavior. The F point of the 
Brillouin zone is the natural zero momentum point, with the spectral weight 
in the b direction symmetrical around this point. Further on, Zwick goes on 
discussing the temperature dependence of the spectra, and how the size of the 
Peierls gap fits with various theoretical models. The interested reader should 
note that the spectral weight distribution of TTF-TCNQ indeed does have a very 
interesting temperature dependence, which is summarized for example in Refs. 
[TH] |150j . However, the conclusions are expected: the experimental results do not 
fit with a Fermi liquid description (an example of this is the absence of the usual 
metallic Fermi edges of the spectral weight; instead a complete supression of the 
quasi particle weight is found), neither do they fit with a strong electron-phonon 
interaction description. Due to the incoherency of the spectral weight and the 
interaction dependent singular lines in the mapped (fc, u) plane, one should use 
ID correlated fermion models in order to explain the experimentally obtained 
spectral features. 

But there is one feature not explicitly touched in that reference, a feature upon 
which our theory triumphs or fails completely: the separation of the charge-type 
cO excitations and the spin-type si excitations, usually referred to as the "spin- 
charge separation" |118j . In contrast to standard Luttinger liquid theory, where 
one studies the spin-charge separation for low lying excitations only, here this 
separation occurs for finite energy excitations. In recent studies of the spectral 
behavior of the ID Hubbard model, such separate charge and spin excitations 
were identified |112j |113j . In our language, however, this separation is mani- 
fested through the different types of contributions which lead to the full spectral 
function, with the branch lines and the border lines as the most obvious features. 

From an experimental point of view, we find that Ref. |137j presents ARPES 
measurements on TTF-TCNQ with details of the spectral features not previously 
reported. In that reference, the ARPES procedure involved a momentum and 
energy resolution of 0.07 and 60 meV, respectively, and with an ambient 
temperature of T = 61 K. Moreover, a comparison with the predictions of density 
functional band theory is made. Unfortunately, the line predictions of this method 
fails completely, producing non physical band dispersions (i.e. predicting spectral 
weight in the {k, u) plane along lines with the wrong bandwidth) and missing some 
experimentally proven spectral features all together. This is illustrated in Fig. 

We have already discussed and shown in section (0), the various line features 
that the pseudofermion description of the Hubbard model produces. In Fig. fl6.3|) 
we show the same grey scale density plot of the experimental spectral weight of 
Ref. |137j . but now fitted with the characteristic branch lines and border line 
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Figure 6.2: A density plot of the obtained ARPES dispersions of TTF-TCNQ 
according to Ref . |137j , together with a density functional band prediction regard- 
ing the singular line features (of the same reference). We see major quantitative 
and qualitative discrepancies, for example the mismatch between the lines a, b 
and c and the experimentally obtained data. Also, we see that the band theory 
is incapable to reproduce the experimental features labeled d and d', which is 
perhaps the most serious failure of this method. 
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Figure 6.3: A density plot of the obtained ARPES dispersions of TTF-TCNQ 
according to Ref . |137j , together with the hne features of the pseudofermion rep- 
resentation of the ID Hubbard modeL The hne between —kp and +kF encirchng 
the symmetry point F, as well as the lines emerging from ±kp and extending to 
higher values of Ic^l, find their exact equivalence in the RHB s-branch line and 
c-branch line theoretical features, respectively, with (U/t) = 4.9, see Fig. ()5.2j) . 
The other line shapes, originating at +kp and fading away as we follow the dis- 
persion towards the symmetry point Z, are exactly matched by the LHB s-branch 
line and c-branch line theoretical features, respectively, with {U/t) = 5.61. Note 
that the latter s-branch vanishes quite rapidly as we depart from the Fermi mo- 
mentum. The dashed line between the symmetry points Z — kp and Z + kp is 
nothing but the {U/t) = 5.61 LHB border line. The LHB features are adopted 
from Fig. ()5.8j) . TTF has a hole concentration of 0.59, and thus TCNQ has an 
electronic density of 0.59. Due to the particle-hole symmetry of our model, the 
spectral weight of both transitions can be mapped onto the same {k,uj) region, 
with u; < 0. 
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features of our pseudofermion theory. These hue features are exactly the same 
as those given in section (jSj) for the relevant values of (U/t), but with the TTF 
spectra folded onto negative energies, illustrating the donor-acceptor relationship 
between TTF (donor, holes) and TCNQ (acceptor, electrons). Note that a dif- 
ferent choice of (U/t) for any of the two materials would produce line features 
not corresponding to the experimental line features. Now, it is straightforward 
to invert the spectral function of Fig. (jS.lljl . from positive to negative energies, 
and add it to the spectral function of Fig. to produce the full theoretical 

spectral weight corresponding to the ARPES experimental result. Figs. ()6.4j) and 

(in3). 

Lastly, we present the density plot corresponding to the spectral function 
provided in Figs. ()6.4|) and ()6.5|) . which is provided in Fig. ()6.6p . These plots 
lets us identify the types of excitations responsible for the total spectral weight 
of the system, hence chartering previously unknown territory. The dominant line 
shapes are thus due to separate charge type and spin type excitations for all 
relevant excitations, not just the low lying ones. Moreover, the relatively large 
spectral weight quite deep inside the band can here be identified with a "velocity 
resonance effect": the border line is the line where the two types of excitations 
propagate with the same group velocity. 
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Figure 6.4: Theoretical spectral weight for TTF-TCNQ, obtained by use of Eqs. 

(jOHll . (lOni) . (gSH), and (I1S21), respectively. This figure is nothing 

but the superposition of the two spectral weight distributions already given in 
Figs. ()5.5|1 and respectively. We expect the singular features of this 

spectral function to be visible in suitable photo emission studies of TTF-TCNQ. 
This figure is used to obtain Fig. ()6.6|) . which indeed confirms the agreement 
between theoretical and experimental spectral weights. 
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Figure 6.5: Theoretical spectral weight for TTF-TCNQ, obtained by use of 
Eqs. fR38|l . (j08|l . (j09|l . (p35|l . (gSl}, and (jUl, respectively. This figure is 
identical to Fig. ()fi.4j) but shows the spectral weight distribution from a different 
angle. 
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Figure 6.6: Theoretical density plot of the spectral weight for TTF-TCNQ. 
Technically, it is this figure which should be compared with the experimentally 
obtained ARPES dispersions. We see that we have an almost complete agreement 
between this figure and Fig. ()6.3|1 . Important features include the two charge type 
excitations both originating at k = kp, one related to the TCNQ stack and one 
to the TTF stack. They both fade away at higher momentum values. In contrast, 
the TCNQ spin type excitation carries much more spectral weight than its TTF 
counterpart, with the opposite relationship for the border lines. 
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Chapter 7 



Conclusions and Discussion 



The pseudofermion dynamical theory reported in this thesis describes the quan- 
tum objects that diagonahze the normal ordered ID Hubbard hamiltonian. In 
the pseudofermion representation, the scatterers and scattering centers are r^-spin 
and spin zero objects, and hence the S-matrix describing the scattering events 
between them is merely a complex number: it is given by the phase shift of the 
ground state — > final state quantum mechanical excitation. Thus, the scattering 
events of these quantum objects are of zero energy forward scattering type |104j . 

The form of the S-matrix, being a one dimensional matrix, is crucial to the de- 
velopment of the dynamical theory. This ultimately results from the "diagonal" 
form of the pseudofermion anticommutation relations. The au pseudofermion 
or pseudofermion hole S-matrix fully controls the one electron matrix elements 
between the ground state and excited energy eigenstates through these anticom- 
mutation relations. Indeed, the anticommutator can solely be expressed in terms 
of the S-matrix, as shown in section ()2.H.5|1 . The form of the pseudofermion 
S-matrix constitutes an important new result of this thesis report and of Ref. 

m- 

The studies of Ref. jHB] showed that the various quantum numbers introduced 
by the Takahashi string hypothesis describe occupancies of pseudoparticles. In 
this reference, the original electrons were related to the "rotated electrons" via 
a unitary transformation V{U/t), as described in chapter Q. The double oc- 
cupancy, the single cr-spin occupancy (cr =1, |) and the no occupancy number 
of the rotated electrons are good quantum numbers for all values of {U/t). The 
separated charge and spin degrees of freedom of the rotated electrons give rise 
to the pseudoparticles. The related pseudofermion description differs from the 
pseudoparticle description by a shift in the discrete momentum values of order 
(1/L), which are associated with the scattering phase shifts due to the ground 
state final state transitions ^Hj . 
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The pseudoparticles have residual energy interactions which makes them un- 
suitable for the development of a djTiamical theory. Indeed, the residual energy 
interaction prevents a pseudoparticle wave-function factorization. However, for 
the pseudofermions, we have no such residual energy interactions, which indeed 
allows for such a factorization This wave- function factorization is valid for 
the normal ordered ID Hubbard model at all values of the energy ratio (U/t), 
filling n and magnetization m. 

In contrast to the usual low-energy Luttinger liquid theory, the theory reported 
here allows us to categorize a separation of the charge type degrees of freedom and 
the spin type degrees of freedom at a finite energy excitation scale. However, in the 
low energy elementary excitation regime, the results of the conformal field theory 
|124j |125j coincide with the dynamical pseudofermion theory, as demonstrated 
inRef. p. 

The pseudofermion dynamical theory presented in this thesis report, was orig- 
inally inspired by the (U/t) ^ 1 methods of Refs. pTj-jSni, where the spectral 
properties of the ID Hubbard model for (U/t) —>■ oo were studied. For arbitrary 
(U/t), the dynamical theory allows us to calculate general closed-form analytical 
expressions of the finite energy one electron spectral weight distributions of a ID 
correlated system (with on-site electronic repulsion). This derivation is done in 
detail in chapters (jHl) and ^ and constitutes important new contributions to the 
understanding of the spectral properties of the ID Hubbard model. This work 
was also presented in Ref. |131j . 

The canonical pseudoparticle-pseudofermion transformation involves a mo- 
mentum shift Q^^{q)/L. This shift is zero for the original ground state, for which 
the au = cO, si pseudofermions have well defined Fermi points iqpau (where l = ± 
denotes the left and the right Fermi point, respectively). The (/c,a;) dependent 
exponents of the theory are then described in terms of the canonical momentum 
shifts of these Fermi points associated with a finite-number electron excitation. 

The dynamical theory applied to the case of one-electron excitations is pre- 
sented in this thesis report for the cases of one-electron removal (RHB) and one- 
electron lower Hubbard band addition (LHB). The closed form expressions are 
explicitly derived in chapter (jH) and plotted in chapter 

The spectral properties of the RHB and the LHB cases can be categorized ac- 
cording to different types of contributions (this categorization is also possible for 
the one electron upper Hubbard band addition but is not presented here). These 
contributions, corresponding to different regions of the {k, uj) plane, each have dif- 
ferent sets of {k, uj) dependent exponents and pre-factors. By this classification, 
we are able to identify practically all features of the spectral weight of the ID 
Hubbard model, in terms of cO and si pseudofermion or pseudofermion hole exci- 
tations. An example is the border line, which is found here and in Ref. |131j . to 
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generate divergent spectral features when the cO pseudofermion or pseudofermion 
hole and the si pseudofermion hole propagate with equal group velocity. 

On a microscopical level, for each of the "contributions", the au = cO, si 
particle-hole towers of states give rise to an orthogonal catastrophe. For any such 
tower of states, we can associate an exponent, which in the case of being negative 
produces a power-law type decay of the spectral weight as we move away from 
the base of the tower. The spectral weight distribution associated with these 
particle-hole processes is controlled by the value of the exponent but also by the 
value of the pre-factor, both of which are constant for one specific tower of states. 
The pseudofermion dynamical theory is capable of explicitly calculating both of 
these quantities, for the entire {k, u) plane. 

Thus on the theoretical side, we are able to predict the position and origin 
of the one electron spectral singular features of a ID correlated metal. However, 
we are also able to connect our theoretical predictions to experimental results. 
The singular behavior of the spectral function, as predicted by the explicitly 
calculated values of the relevant exponents, leads to a spectral weight distribu- 
tion which should be detectable by photo emission and / or photo absorption 
experiments. It turns out that within the approximation of only considering 
the leading order elementary processes to the RHB and the LHB one-electron 
spectral weight, we are able to reproduce, for the whole energy bandwidth, the 
experimental spectral distributions found for the organic compound TTF-TCNQ 
by high- resolution ARPES. These new results are presented in chapter (jH)) and in 
Ref. |148j . The TTF-TCNQ high- resolution ARPES experiments were reported 
in Refs. jTHI PZ] |137j . In conclusion, the dynamical theory presented here allows 
for an understanding of the elementary quantum processes that give rise to the 
spectral features of TTF-TCNQ. 

With the advent of new experimental techniques that allow for a high-resolution 
study of the spectral features of quasi- ID materials, the pseudofermion dynam- 
ical theory has yet many challenges ahead of itself. One of the most exciting 
recent experimental setups is the "optical lattice" in which ultra-cold fermions 
are trapped in a potential well, forming a "real" ID quantum chain. These sys- 
tems can be described by the ID Hubbard model, with the electrons replaced by 
ultra-cold fermionic atoms. Even though some preliminary results already have 
been reported this technique is still at its infancy. However, it does 

promise an unprecedented control over the necessary parameters (such as the on- 
site Coloumb repulsion U and transfer amplitude t) , enabling the high-resolution 
measurements necessary for a complete understanding of these materials. The 
theoretical spectral weight expressions obtained in this thesis report should be 
taken into account when characterizing the experimental spectral features ob- 
tained through this method. Indeed, this experimental setup comprises one of 
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the most exciting future applications of the pseudofermion dynamical theory. 

In conclusion, we see that the pseudofermion dynamical theory is a suitable 
theory for the study of the spectral properties of the ID Hubbard model, yielding 
results in good agreement with the behavior of the one electron spectral function 
in the {U/t) — >■ oo limit (ZSj-iSO], and the low-lying elementary excitation limit 
ins and with experimental results HE] [HI [MI- 
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